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ABSTRACT 



A Mechanical Engineering Department project in which the VAX/VMS system was 
utilized to create an interactive menu driven program to solve basic preliminary ship 
design problems. 

Enhancement of an existing program w'as initiated to improve the user interface by 
adding user-friendly help information. Also, routines were written to calculate 
propulsive power requirements based upon the ship form coefficients selected and com- 
parisons made using the Method of Silverleaf and Dawson and the Admiralty Coefficient 
prediction method. Further computational routines were added to predict range and 
endurance figures for estimated voyage data and selected propulsion plant types, using 
the U.S. Navy Design Data Sheet DDS9400-1 methodology. Finally, the detailed 
printed report generated by the system was updated to include reports of these calcu- 
lations for the users design study. 
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I. INTRODUCTION 



A. BACKGROUND OF GRAPHICAL USER INTERFACES 

Advances in the electronics industries since the 1970's have enabled the development 
of increasingly more capable and complex computer systems. To harness these systems, 
software engineers have produced programs that enable todays users to have a degree 
of computational power that was unfathomable just 10 years ago. As a result, software 
programs have been increasing in complexity every year. This has turned out to be a 
double edged sword; on the one hand the desktop computer power has enabled the de- 
velopment and solution of problems that would require the resources of a mainframe 
computer not too long ago, but it has also meant that the user is confronted with the 
requirement to comprehend all of this capability. This has led to the development of 
user interfaces that seek to simplify and manage the task confronting the user [Ref 1: 
pp. 115-139]. Among these user interfaces have been a host of graphical systems that 
attempt to simplify the program presentation to the user by means of visual icons and 
multiple screen windows for specific tasks. These interfaces are termed Graphical User 
Interface, or GUI for short. 

The very first development of the GUI occurred in the late I970's at XEROX Cor- 
poration's Palo Alto Research Center (PARC) and was embodied in the XEROX Star 
Office System computer [Ref 2]. The Star System was not agressively marketed and 
languished in the PARC laboratories until 19S4 when its basic ideas were consolidated 
into the APPLE Corporation's Lisa, and subsequently the Macintosh system, which was 
ultimately more commercially successful and exists even this day. Since that time several 
other GUI's have been developed and popularized for the express purpose of interfacing 
between man and machine [Ref 1 ; pp. 337- 371]. 

B. 'TOOL BOX' AND PRELIMINARY SHIP DESIGN 

This brings us to the point of developing a program that will enable a user to solve 
basic preliminary ship design problems without necessarily being completely familiar 
with the inner details of the program code. Having at our disposal a Digital Equipment 
Corporation (DEC) VAX/V.MS Network, we have taken a basic ship design program 
previously developed [Ref 3 ] at the Naval Postgraduate School and enhanced its use- 
fulness by the careful application of grapliical presentations. This program as it existed 
previously computed some basic ship static and hydrodynamic stability parameters. 
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predicted the ship's turning circle, and printed a detailed hard copy report of those cal- 
culations if so desired. It began the design process utilizing basic GUI principles but 
did not reach its potential with respect to its ability to communicate with the user. This 
does not at all reflect upon the abilities of its developer, but only indicates the scope of 
the software developement task. In the short time available to the preliminary author, 
he was able to implement many of the VAX Systems User Interface Services (UTS) 
routines in the initial version. However, time and other constraints kept him from a full 
implementation of the VAX systems capabilities for graphical presentation. This is the 
point where this author has engaged the project and made some important enhance- 
ments of the basic code and also some additions to its computational abilities. 

Figure 1 on page 3 (from Ref 4, p. 15), shows a generally accepted version of the 
Iterative Design Spiral, typically used in Preliminary Ship Design to refine a design as 
conditions and data become more concrete. TOOL BOX , as it existed, sought to define 
the early parameters required for the design, and attempted to show quickly that the 
design would be stable and maneuverable, essentially completing the boxes labeled (a) 
and (b) in Figure 1. From this point we will progress further along the iteration pattern 
and provide the calculations indicated in box (c) of Figure 1, and go back and try to 
enhance some of what was accomplished before. The modular concept of TOOL BOX 
supports this type of application where the design calculations will be performed time 
and again, as the design data is finalized. Thus each module will provide the specific 
calculation required by one of the Design Spirals major elements, and the modules can 
be accessed in any order required. 

Additionally, we will add some commonly utilized GUT features to enhance the 
programs ease of use and information presentation, and attempt to make full use of the 
VAX'V.MS Operating Systems capabilities for Graphical Programming in High Level 
Languages. 

C. WHY WINDOWING SCHEMES? 

One common feature shared by most, if not all GUI's, is the ability to display mul- 
tiple windows on screen simultaneously, with separate windows for separate tasks. This 
scheme has the advantage of concentrating the user's attention on a certain window that 
is considered the active window. Windowing GUI's are present in all forms of computing 
environments today, from Microsoft Windows for MS-DOS based systems, to the Apple 
Macintosh, to a variety of systems for UNIX and its derivatives, and many others. 
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Figure 1. Iterative Design Spiral (typical) 



However windowing systems do allow for the presentation of information to the user in 
a manner not obtainable by command-line systems. 

For example, one of the most powerful ship design programs available today is 
ASSET from Boeing Company, based in Seattle, Washington. This program employs 
strictly a command-line interface where the user inputs various commands to the system 
and awaits a prompt to continue. Nevertheless, for all the power available in ASSET, 
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it suflers from this command-line interface in the worst fashion. A user either has to 
know the commands available and all the options for them, or they must have the ref- 
erence manual handy at all times. Since manuals are generally not available for all users 
at once, it is considerably harder to master all of its features without constant use and 
practice. 

My primary goal in this project was to enhance the methods used by TOOL BOX 
in which information was presented to the user by utilizing some basic GUI presentation 
methods such as 'pop-up' dialog boxes that would convey some important or useful in- 
formation to the user at the appropriate time, then disappear when not required. This 
form of windowing has the advantages of getting the user's attention quickly since the 
display screen changes in a most abrupt fashion. Since the user will now be concen- 
trating on the small piece of information conveyed by the dialog box, and any options 
will be displayed for them at this time, no requirement for reference manuals or mem- 
orization is needed. This technique is employed by all the GUI's available today, so the 
method is consistent with current software design practice. GUI traits that we will in- 
corporate will be timed display of choices or information, visual prompting for needed 
input, multiple display of windows, and so forth. Also, we shall improve the 
calculational power embodied in TOOL BOXhy adding other routines based around the 
same central core code and utilizing the same interface. 

D. STANDARDS AND CONVENTIONS 
1. Attribute Blocks 

Attributes are used to describe the basic default format representation of text 
and graphics. These attributes are stored in blocks that can be modified by the pro- 
grammer. Whenever any modification (for example the screen font) is desired, it is 
necessary to copy the basic attribute block 0 to another numbered block and then 
modify the attribute concerning screen fonts. Otherwise a previously modified block can 
be used as the basic block and modified again. Attribute block 0 can never be changed, 
as it contains all the system defaults. All previous modifications remain unless changed 
by the current modification. The system programmers references [Ref 5 Chapter 9, 
section 2.2 ] contain more on what modifications can be performed and how to accom- 
plish the desired changes. The previous author has utilized blocks I thru 10 for his 
preliminary code; in order to ensure that we don't mistakenly modify a block that is in 
use in some other undocumented portion of the program, we will use blocks 20 thru 50 
for our additions to the code. The documentation for the previous section of code is 
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incomplete, thus we will skip blocks 1 1 thru 19 as a safety feature. The DEC VAX/ VMS 
system allows up to 255 different attribute blocks, so we should have no problem utiliz- 
ing separate blocks for this part and any subsequent additions. Optimization and min- 
imization of code and resources is not an objective at this point in the development. 

2. Code Documentation 

This is the standard we will use to add documentation and comments into the 
code section to explain each step to a reader. All too often a source code that has not 
been documented internally by the author will be reviewed by another programmer and 
the lack of clear concise documentation will lead to lengthy frustration. Therefore, our 
method will be as follows; 

• Comments in the source code shall begin with C*****, and every line that is not 
intuitively obvious will have a comment preceeding it. 

• Subroutines will be preceeded by a description of the purpose and utilization of the 
routine. If the routine is called by more than one section, it will be noted. 

• New sections or changes added by this author to existing code will be annotated 
as changed, so credit may fall to the appropriate author. 

• Complete copies of the TOOL BOX source code, including parts not covered in this 
document, may be obtained from Professor F.A.Papoulias, Code ME, Pa, Naval 
Postgraduate School, Monterey CA 93943-5000 or from Professor J.F. Hallock, 
Code ME/HI, Naval Postgraduate School, Monterey CA 93943-5000. 
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II. DIALOG BOXES 



A. FUNCTION OF THE DIALOG BOX 

The use of 'pop-up' dialog boxes is fundamental in the operation of most GUI's in 
that, at some point during the execution of the program, some element of information 
will need to be conveyed to the user. TOOL BOX is designed in this manner to overcome 
the limitations of comntand-line based ship design analysis programs thus ever}’ effort 
will be made, mostly via context sensitive dialog boxes, to present information to the 
user as it is required without requiring the reference manual to be handy. Basically the 
dialog boxes function as follows: 

1. The point of the program that the author deems important calls the dialog box 
routine. 

2. The dialog box is presented, displaying the information or choices available to the 
user. This involves VAX LTS routines to create a display, to create a window, to 
modify the text as the author desires for presentation, and to execute its function. 

3. If some action is required by the user, the program waits for that action. 

4. If only information is to be displayed, then a suitable time is given for it to be 
viewed by the user. 

5. Upon completion the dialog box is removed from view by the system. This step 
involves VAX LTS routines to time the display to completion, delete the display 
(and associated keyboards, windows, and viewports) and remove itself from the 
screen. 

The best resource for the actual routines used and the parameters required is the 
specific system manual [Ref 5 Chapter 18]. This is an indispensable source for the 
VAX/VMS programmer. 

B. TYPES OF DIALOG BOXES 

Two distinct types of dialog boxes were required to be used. The similarities are 
apparent, yet each functions differently enough to warrant separate discussion. 

1. Visual Display Information Box 

This type of dialog box is used to display information to the user at the appro- 
priate time. No conversation between the user and the system is required, therefore the 
routine can be fairly simple. An example of this type of display is the subroutine 
SHOW_SAVE used during any attempt by the user to save a data file. Its purpose is 
simple: it reminds the user not to save the file under the same name and extension as 
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was used previously in another part of the program, since this will cause the new data 
file to overwrite the old one. This subroutine code is simple and is reproduced in Ap- 
pendix A. Essentially the subroutine is called during an attempt to save a data file by 
the system. It pops up and writes to a new display box and window the message about 
using a different filename and extension. It does not need any action by the user, as it 
is for informational purposes only. Since no action is required, it must remove itself 
from display after an appropriate time has passed, yet give the user enough time to read 
and comprehend its message. For this time delay we have used a system library call to 
the VAX function LIBSSPAWN( 'WAIT 00:00:15'). This function 'spawns' a new 
process that calls the intrinsic VAX function WAIT for a delay set as 15 seconds. The 
time format is HH;.\1M;SS in hours, minutes, and seconds. Upon completion of this 
delay, the display window is deleted and the system proceeds at the point the delay 
process was spawned. This method is rather elegant in that only system functions are 
used for the delay. The UTS routines required for this action are 
UIS$CREATE_DISPLAY, UIS$CREATE_WINDOW, UIS$SET_FONT, 

UISSSOUND_BELL, UISSTEXT, and UISSDELETE_DISPLAY. The font and bell 
routines are optional; however, they add to the display by making it both audible and 
visually dilferent than the existing display when the window is created. Also, the bell 
volume is fully variable by the programmer. A value of '0' for the second argument will 
set the volume to its most quiet position, and an '8' will set it to the loudest position. 
We have used a value of '4' which creates a medium loud bell, loud enough to draw at- 
tention but not to cause the user to be startled. If no volume is specified, then the sys- 
tem defaults the bell volume to the workstation volume setting. The window created by 
this process is shown as Figure 2 on page 8 below. 

The presentation of this display lasts approximately 15 seconds depending upon 
the number of windows and processes currently open on the terminal, and so will be 
slightly longer if multiple concurrent process are active. Upon removal of itself, the 
display returns to the window that was active when the dialog box routine was called. 

2. Conversational Dialog Bo.x 

This type of dialog box is used to display information to the user and also to 
give the user a choice for proceeding. An example of this type of procedure is the sub- 
routine NOWHERE that is called by the mouse selection of a Reserve Module from the 
Main Menu. This routine uses essentially the same UIS routines as the previous 
method, but also incorporates a FORTRAN PAUSE statement that interrupts the 
progress of the program. Using the PAUSE statement here gives the user two choices: 
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Figure 2, Visual Display Information Dialog Box 

1. Enter 'CONTINUE' or 'C' and the program will resume at the point of inter- 
ruption. 

2. Enter 'EXIT' and the program will terminate at that point, deleting all displays and 
windows associated with it, and return the user to the active process that was used 
to call the program TOOL BOX. 

This selection of choices was chosen since at this point the user is in the Main 
Menu Window and the direction to proceed w^ill be either, (1) to continue if that selection 
is entered, or (2) to end the program and exit if the user has selected these future mod- 
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ules just out of curiosity. The user need not see what has been typed at the terminal 
when a choice is made; however, if the calling window is visible at the bottom of the 
screen as in Figure 3 on page 10, then the typed text reply will appear here. This is only 
necessary if the user happens to be a poor speller, since the program will not proceed 
without the entry of a properly spelled command. Thus the VAX shorthand command 
for ' CONTINUE' as 'C' is best utilized. Otherwise, the user enters a choice by typing 
on the terminal keyboard followed by a carriage return, and the program proceeds along 
that selection. The result of this process is showm in Figure 3 on page 10 after presen- 
tation of the dialog box. 

No timing of the display is required since it will remain there until a proper 
command is entered. The routine UISSSOUND_BELL is used again to signal audibly 
the user and gain his attention; the bell volume again is set to a medium volume at '4'. 
The code for this routine is reproduced in Appendix A. 

3. Other Enhancements 

Although we have not used the technique here, is is possible to have a conver- 
sational dialog box that shows the users reply in the window that requests it. For this 
routine, we would use essentially the same format as in the subroutine NOWHERE, 
except we would attach and enable a keyboard to the display by utilizing the LTS rou- 
tines U1SSCREATE_KB and U1SSENABLE_KB. Then a small routine would have to 
be added to accept and display keyboard input into the window. This is not dilTicult; 
and it might be helpful by allowing the user to see what has been typed as their reply, 
thereby alleviating some errors caused by typing mistakes. This technique has not been 
used here since there are only two allowable responses when using the FORTIUAX 
PAUSE statement; CONTINUE and EXIT. Any other response will be ignored without 
exception and no action will be taken by the program until a proper response has been 
entered. Upon completion of the routine, the attached keyboard would have to be dis- 
abled and deleted by LTS routines and the text reply passed to the calling program as 
input for some action. This sequence is possible and some use may be found for it later 
in the development. 



9 







JE . 


HCIPWINOOU ' 




- - - THIS FEATURE NOT AVAILABLE --- 
TYPE •• CONTINUE” TO RESUME 


1 




1 STATIC 5Ta| 


OR • 


EXIT” TO END PROGRAM 


1 


1 






1 




MANEUVERIN( 






1 



POWER PREDICTION 



endurance 




FILE UTILITIES 



A program for preliminary 
analysis of 

SHIP CHARACTERISTICS 

PROFESSOR F. A. PAPOULIAS 
AND 

GERALD K. MCGOWAN 
DAMES R PLOSAY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY. CALIFORNIA 
1989/90 



GENERATE RERQRT 



EXIT THE PROGRAM 



I $RUN TOOLBDX 
rORTRAN PAUSE 



Figure 3. Conversational Dialog Box Display 



10 



III. PROPULSIVE POWER REQUIREMENTS 



A. BACKGROUND 

One of the more useful pieces of information to have during prehminary ship design 
is an estimate of the vessel's propulsion plant power requirements. This is especially 
handy since the size of the plant required to drive the vessel will also give an initial es- 
timate of volume required to enclose it, manpower required to operate it, and the cost 
to purchase and maintain it. So a single 'ball-park' figure will go a long way towards 
helping to advance the preliminary’ design. 

One of the foremost sources of information and research conducted in this calcu- 
lation was Rear Admiral David W. Taylor, USN(ret,dec), for w’hom the David Taylor 
Naval Ship Research and Development Centers (DTRC) in Carderock and Annapolis, 
Maryland were named. His studies for the then US Na\7 Construction Corps at the 
beginning of this century still serve as required reading for students of ship design and 
naval architecture [Ref 6]. Some of this research has been duplicated and enhanced by 
more recent scholars utilizing more powerful calculational abilities than Adm. Taylor 
ever envisioned. These studies [Refs. 7, 8 , 9, 10] have resulted in some variations of 
Adm. Taylor's basic equations and will be examined here for two methods of power 
prediction that we will utilize. 

B. METHODOLOGY 

Many methods of power prediction exist in the literature and a short amount of re- 
search will produce most of them. The list of methods begins with the 'Quick and Dirty' 
method of Admiralty Coefficients, which assumes that an already established hull form 
(parent ship) with its basic parameters is known and that it is simply being modified for 
the new design. This method has been slightly modified [Ref 8 : pp. 308- 310] and is 
used here in that form. There are several graphical methods for power prediction that 
require charted information on the specific hull form or a family of hull forms close to 
the design in question. These methods will not be used since those charts would have 
to be available in digitized form for all ship classes of interest. There is also an entire 
class of prediction calculations based upon the total estimated hull form resistance to 
movement through the water. This method will not be used here because it entails much 
more detail than we desire, or have at this stage of design. Then there is a method ori- 
ginally proposed by Silverleaf and Dawson [Ref 9 : pp. 167- 196] and subsequently 



modified by Stian Erichsen [Ref. 10 : pp. 83- 115] in a University of Michigan College 
of Engineering departmental paper in 1971. This method will also be presented here and 
the user can choose the desired result with care. 

Some notes about the methods presented here must be discussed first, in order that 
no potential conflicts arise later on: 

• These methods were developed for monohull surface craft only; no capability to 
predict power requirements for catamaran, trimaran, SWATH, etc., types of vessels 
is intended. 

• The resulting power requirement figures do not include added resistance due to 
waves, fouling, or other environmental factors. 

1. Method of Admiralty Coefficients 

The first method to be used is based upon the Admiralty Coefficient since a 
great deal of information about that Coefficient is available in the literature. This 
method assumes that the ship's resistance is all frictional, and thus the power required 
for propulsion varies proportionately as the cube of the ships speed. The Admiralty 
coeflicient for a vessel when no parent vessel is used is given by : 



A,~i.lOUL +^) 



(I) 



where L is length of vessel in meters between forward and aft perpendiculars and V is 
speed of vessel in meters per second. Otherwise, if a parent vessel is available, the 
Admiralty coefficient is given as 



Ac- 



^2/3 j/3 



( 2 ) 



where P, is parent ship power in kilowatts needed to make a speed of V meters per sec- 
ond and A is parent ship displacement in metric tonnes. 

Then, using that coefficient to characterize the vessel, we take 



^2/3 J.3 






Ac 



( 3 ) 



where P„ is power in kilowatts of the new ship needed to make the same speed V and A 
is the new ships displacement mass in tonnes. 
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An improved version of this formula that incorporates almost all of the essential 
ship parameters except that it neglects the influence of the block coefficient Cg is given 
by 

F^(33.0-0.017L) 

P = (4) 

15,000- 

where n is propulsion shaft revolutions per second needed to make a speed of V meters 
per second and all other parameters are as described previously above in metric units. 
This form is easily equated to our more common form of English Horsepower (hp) units 
by a simple conv'ersion factor, and is the equation we will use for our first method of 
power calculation in the program. Also, some obvious limitations to this formula ap- 
pear readily enough and should be noted ; 

• L cannot be larger than 1941.177 meters or the numerator becomes a negative 
value. 

• The product of shaft speed n in revolutions per second and the square root of L 
cannot exceed 136.364 or the denominator becomes a negative value. 

• A realistic maximum shaft speed of 4 revolutions per second ( 240 rpm ) limits ship 
length even further to 1162.196 meters. This is not too much of an overall limita- 
tion, but it is present nonetheless. 

These limitations will not affect us for the most part but should be noted in any 
event as they could end up being important. It should be noted further that this formula 
(4) tends to fit special cases and types of ships. When used outside of these allowable 
ranges, the results tend to produce erroneous expectations. According to Harvald, 
[Ref 7 : p. 290], "It is not often mentioned within which area the formulas can be ap- 
plied". Thus we have a method that can be quite accurate for power prediction, if we 
happen to fall within the allowable range, and which otherwise could generate a false 
prediction. For the first method however, it will be satisfactory for an initial estimate. 
We can then couple these results with the second method. If the results are close, we 
can assume that our prediction is reasonable. If they depart, then we will have to use 
them with caution. 

2. Method of Silverleaf and Dawson 

The second method that we will use for power prediction is given by Silverleaf 
and Dawson [Ref 9] as modified by S. Erichsen [Ref 10). It uses almost all of the ship 
specific parameters of the design and is usually highly reliable for ships within the fol- 
lowing ranges of speed - length ratio VJfL and beam - draft ratio BIT : 
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( 5 ) 



0 . 4 <(^^)< 1.2 



and 



2.0 ^(y)^ 4.5 

with Q and length - beam ratio LjB for single screw ships 

0.50 < Cg < 0.86 

3.33 <(-^)< 9.50 

D 



or for twin screw ships 

0.54 < Q < 0.80 

3 . 80 ^( 4 -)^ 11-50 

D 



( 6 ) 

( 7 ) 

( 8 ) 

( 9 ) 

( 10 ) 



Where the parameters are specified as follows for these and the following equations : 
is ship speed in knots 

Vg is the ships Boundary Speed, which is that speed, for a given hull form, below 
which the resistance coefficient does not vary greatly, and above which it begins to 
increase rapidly, [ Ref 9 : p. 168 ]. 

L is ship length between perpendiculars in feet 

B is ship beam in feet 

T is ship draught in feet 

A is displacement mass in tons 

Q is the Block Coefficient given by 



r 35A 
® LBT 



for ships in standard seawater. 

The Silverleaf and Dawson formula then is given by 

_ {\+xWl^vlK^rt^yK^ 

427.b/o(f/4oo-^LS/>) 



( 11 ) 



( 12 ) 
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where Pj is delivered horsepower required at the propeller. Equation (12) will be used 
as the second method of power calculation in the program development. The terms in 
this formula are all specified for design parameters that may or may not be same as the 
model from which the formula was derived. Thus some corrections need to be applied 
to each parameter below. For a design that does not exactly match the model design, 
then : 



Utilizing ship speed at the Boundary’ Speed Vg expressed in knots 

Us = (1.7-1.4Q)Vr 



(13) 



Correcting for ship Length : 



(1 + .r) = 



0.85 



for L > 1000.0 



0.85 + 0.001 85 ' 6 fori < 1000.0 



(14) 



Correcting for the Beam - Draft ratio 



Kbt — 



0.982 



0.96 + 5.4=" 10“^10°-^’T) 



for BjT<2A 



for BIT >2.4 



(15) 



Correcting for ship speed F* not equal to the Boundar>^ speed Vg : 



K^ = 



2.1S-1.25{-f-) + 5.5{^f 



V, 



for U/Ug< 1 



21.2-43.2(-j7^) + 23.0(-f^)^ forU/Kg>l 



(16) 



Correcting for a ship hull that benefits from a design with a bulbous bow 



Kd — 



1.00 



0.95 



for ships without a bulbous bow 



for ships with a bulbous bow 



(17) 
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Correcting for Hydrodynamic Efficiency of the hull form at L = 400.0 feet ; 



^^400 ~ 



2.60-0.2917- 



V 






2.38 -0.2917- 



. 1/6 



for a single screw 
for twin screws 



( 18 ) 



with a correction for ship length L not equal to 400.0 feet 



^LBP — 



0.9196 + 2.3rHO~^L - 7.5*10~*L^ for L 400.0 



1.0 



otherwise 



(19) 



Correcting for open water efficiency at propeller speed not equal to 120 revolutions per 
minute, where : 



Vo ~ Vo,i 20 + *^Vc 



( 20 ) 



with 



Vo , 120 ~ 



0.98 - 0.55Cg 



0.90-0.33Cfl 



for a single screw 



for twin screws 



( 21 ) 



and for propellor rotating speed N revolutions per minute 



(5;/o = 



0.360 - 0.0029iV 



for a single screw 



( 22 ) 



0. 1 35 - 0.001 1 N for twin screws 

Thus we have all the required parameters for the Silverleaf and Dawson formula. 

This formula tends to be more accurate over the entire range of ship specific 
parameters. The parameter limitations are delineated in the beginning and cover a wide 
range of vessels. 
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C. COMPARISON OF THE METHODS 



A quick comparison of the two methods utilizing equations (4) and (12) as generated 
by a test code of the formulas and data from established merchant designs [Ref 4 : pp. 
137- 171] is presented in tabular form in Table 1. 



Table 1. COMPARISON OF POWER PREDICTION METHODS 



Ship Design 


Desiiin Power 
^hp) 


Admiralty 
Method (hp) 


S & D Method 
(hp) 


Large General Cargo 


24.000.0 


20.351.91 


23,780.55 


Container 'A' 


17,500.0 


14,127.90 


11.503.22 


Container 'B' 


32,000.0 


30,511.81 


29,480.57 


Roll On- Roll Off 


37,000.0 


28,573.58 


27,435.27 


LASH Barge Carrier 


32,000.0 


26,203.87 


23,787.66 


SEABEE Barge Carrier 


36,000.0 


26,772.95 


21,918.27 


Tanker 'A' 


15.000.0 


12,045.93 


11,077.69 


Tanker 'B' 


45,000.0 


53,068.73 


36,030.91 


LN'G Tanker 


43.000.0 


39,727.77 


37,910.40 


Bulk Carrier 


15.300.0 


10,940.75 


8,578.57 


Ore, Bulk Oil Carrier 


24.000.0 


21.700.58 


18,981.25 



Imnrediately we can draw some conclusions of the two methods used here to calcu- 
late propulsive power requirements; 

• In all cases except one the Admiralty Method has underestimated the power re- 
quirements, although after an allowance for a service and fouling margin, that 
underestimation is not excessive. The one case where this estimation did exceed the 
actual power requirements ( Tanker 'B' ) the result was not grossly excessive, and 
was in error of only 15 percent. 

• The Silverleaf and Dawson Method also will always underestimate the power re- 
quirements by a slight margin. This margin of error for underestimation is usually 
not excessive and can be accounted for by service, fouling, and other margins. 



A test program was used to generate the data for Table 1 since the TOOL BOX 
screen displays would consume excessive space in this paper. 

Also, the following should be noted : 

• These values are computed for essentially calm seaways only. No added resistance 
from waves in a rough sea is predicted by these methods. 
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• The figures shown are not the total plant requirement for propulsion, only that di- 
rectly needed at the propeller to drive the ship. Thus margins for gear train losses, 
service losses, fouling losses, etc must be added for a total plant size estimation 
figure. 

• And similarily, from Silverleaf & Dawson's formula there is also a correction for a 
ship design that includes a hydrodynamically well designed (WD) bulbous bow. 
For that case, the power requirement for propulsion is 95% of the power require- 
ment shown, since the improved bow shape lessens the hull form resistance. Those 
figures are also listed in Appendix C on the next line directly below the S & D 
calculation. 

Overall though, the two methods taken in combination do allow a fairly accurate 
prediction of a designs propulsive power requirements. The results generated by this 
module of TOOL BOX can thus be used to size and estimate the machinery box volume, 
estimate the manpower required to maintain and operate the machiner}^ plant, estimate 
the cost of the propulsion plant, and so forth. 

D. THE TOOL BOX POWER PREDICTION MODULE 

Implementation of this algorithm of code into the existing TOOL BOX framework 
constituted a fairly easy and straight forward process. An examination of the subroutine 
methodology of the previous modules was made and the same basic pattern was fol- 
lowed. The resulting code is shown in Appendix D and borrows many elements of the 
STATIC STABILITY module (whose flow diagram was utilized) in addition to the vis- 
ual display dialog boxes for HELP WINDOWS that present pertinent information to the 
user and then disappear. 

1. Operation of the Power Prediction Module. 

The procedure used to operate the Power Prediction IVIodule is essentially a 
logical progression through the available menu choices. We will progress through the 
screen displays in an example designed to show all the relevant points of the TOOL BOX 
Power Prediction Module. It is understood that when the user is instructed to enter a 
response to a program prompt, that they will follow-up that action by depressing the 
Return key on the VT220 LK201 keyboard. That action is needed to complete any data 
entry operation. 

2. An Instructional and Illustrative Example 

Shown in Figure 4 on page 19 is the TOOL BOX Opening Menu as the pro- 
gram is initiated by the user. This screen is common to the core code of the program, 
and all the program modules are initiated through this screen using the mouse device. 



18 



MAIN MINU 



STATIC STABILITY 



MANEUVEPING 



POWEP PREDICTION 



endurance 



RESERVE MODULE 



FILE UTILITIES 



GENERATE REPORT 



EXIT THE program 



VERSION 18.01 



A program for preliminary 
ANALYSIS OF 

SHIP characteristics 

PROFESSOR F A PAPOULIAS 
AND 

GERALD r . iCGOlilAN 
JAMES R PLOSAY 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 
1989/90 




Figure 4. TOOL BOX Opening Menu Screen Display 

The third mouse selection available to the user is now POWER PREDICTION 
and a click on this selection (when the mouse pointer is in that box and the text has 
become reverse highlighted) will initiate this module and call up its associated display 
screen. The screen resulting from this selection is shown as Figure 5 on page 20. 

As Figure 5 is displayed to the user, it should be noted that nothing is shown 
in the DATA DISPLAY AREA on the right hand side. That is because no data is 
passed in at initialization time. A mouse click on the first menu choice ENTER DATA 
FROM KB will initialize the data entry subroutines and the program defaults will dis- 
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POUIft mtOICTION WINOOU 






OPTIONS ^''£NU 



CNTCR DATA FROM KB 



input data from file 



STORE DISPLAYED DATA 



plot grapm to screen 



PLOT GRAPH TO DISK 



RETURN TO MAIN 



EXIT THE program 



RESPONSE hPEh 



INSTRUCTIONS 



DhTh DISP^hY hPEh 



] 



(SPUN TOOLBOX 



Figure 5. POWER PREDICTION MODULE Opening Menu Screen Display 

play in the data region as shown in Figure 6 on page 21. Otherwise the user can click 
the mouse pointer on the second box INPUT DATA FROM FILE and load a prede- 
fined data file, perhaps one used before or created by the instructor. This sequence is 
shown as Figure 7 on page 22, with an existing data file entered in the RESPONSE 
AREA to be loaded for data entry. 
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routft pkloictioh uinoou 



m 



OPTIONS h'ENIJ 



INPUT DftTft FROM FILE 



STORE DISPLAYED DPTft 



PLOT GRAPH TO SCREEN 



PLOT GRAPH TO DIS>' 



RETURN TO MAIN 



EXIT THE program 



RESPONSE hPEh 



INSTRUCTIONS 



ENTER ft LINE NUMBER 
OR [RETURN] TO EXIT 



DhTh DISPlhY hPEh 



1. 


SHIP NHME 




GENCARGO 


2. 


LENSTH (Lpp) 


. . . (FT) . . 


502.50 


3. 


DESIGN DRftFT 


. . . (FT) . . 


35.0 


4. 


BEAM 


. . . (FT) . . 


02.0 


5. 


DISPLACEMENT 


. . (LTONS) 


31995 0 


6. 


BLOCK COEFICIENT. . . 


. . . (Cb) . . 


0.66980 


7. 


SHIP SPEED 


. . . (KTS) . 


20.00 


B. 


PROPELLER SPEED 


. . . (RPM) . 


120.0 


9. 


NUMBER OF SHAFTS. . . 




10 



gp R£Q[) XHPJTS*^*^*’*’*^* 



•^♦♦♦horsepower . (HP) •♦•♦♦♦♦♦•♦ 

METHOD; 

HDMIRRLTY 

S.&. D.UJ/.NO BB 

5.6.D.U)/ BB 

BB= WELL DESIGNED BULB BOW. . . . 



20352 

23775 

22506 



ISRUN TOOLBOX 



Figure 6- DATA ENTRY Display Screen at Initialization of KB Entry 
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Figure 7. INPUT DATA FROM FILE Filename Entry 
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After the user clicks the mouse pointer on the ENTER DATA FROM KB se- 
lection, the program calls in the stored default values of ship parameters. This is done 
for several reasons: 

• It makes little sense to see all zeroes for data by a new user. Default values uill 
give the user an idea of what is required. 

• There are so many restricting data checks on allowable entries in the algorithm 
( eqns. 5- 10 ) that the resulting HELP displays would consume inordinate 
amounts of time (if all initial entries were zeroes). This is because as we progress 
from zero length, beam, draft, displacement, block coefTicient, ..., etc., to our final 
values the error checking codes would be called in with almost eveiy parameter 
change due to the large changes in magnitude from zero to any reahstic value. 

• It speeds up the processing of the initial calculations. 

Proceeding, note that, in the lower left hand INSTRUCTIONS box of 
Figure 6 on page 21, the user is prompted to enter a line number for whichever line it 
is desired to modify. This can be entered as an integer (I) or a real number (1.). The 
program code does not care specifically one way or another. Suppose we wish to modify 
the existing data from the set displayed, which comes from the LARGE GENERAL 
CARGO vessel of Table 1 on page 17, to another set from that same table for the 
TANKER B design. First we would need to change the data file name to a unique value 
for that design. We would enter a I for the line to change, then enter the name of our 
design. Notice that the user is prompted for the exact detail required for the corre- 
sponding Line number chosen, as in the INSTRUCTIONS box at the lower left corner 
of Figure 8 on page 24. So here we'll use BTANKER for our design name, as shown in 
Figure 8 above the INSTRUCTIONS box, where the user will also note that their data 
entry is mimicked in the RESPONSE AREA on the left middle screen box as shown. 

This is one feature that will help the user : to have all input echoed back to the 
screen as it is being entered, so that changes can be made on the spot if required. After 
that entr>’ is completed, this data will be written to the screen corresponding to the Line 
entry that was modified, in this case Line 1. 
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OPTIONS KENLI 



DhTh DISP-hY hPEh 



input d«t« from file 



STORE DISPLAYED DATA 



PLOT graph to screen 



PLOT GRAPH TO DISK 



RETURN TO MAIN 



EXIT THE PROGRAM 



RESPONSE HP Eh 



BTftNKER 



INSTPUCTIONE 



ENTER NEW SHIP NHME 
[RETURN] TO EXIT 



1 


SHIP NRME 




GENCARGO 


2. 


LENGTH (Lpp) 


. . . (FT) . . 


582.50 


3, 


DESIGN DRHFT 


. . . (FT) . . 


35.0 


4 . 


BEHM 


. . . (FT) . 


B2.0 


5. 


DISPLACEMENT 


. . (LTONS) 


31995.0 


6. 


BLOCK COEFICIENT, . . 


. . . (Cb) . . 


0 . 66930 


7. 


SHIP SPEED. 


. . . (K.T3) . 


20.80 


B. 


PROPELLER SPEED 


. . . (RPM) . 


120 0 


9. 


NUrvBER OF SHAFTS. . . 




1.0 



QP f^£QD 



♦^^RESULTS^^* 

4»»«.«.»#.|^QP5£PQ^£P (PIP) **«-*«•*»»#» 

METHOD: 

HDMIRftLTY 20352. 

S.&.D.UI/.NO BB 23775. 

3.&.D.W/ BB 225B6. 

BB= . UDELL . DESIGNED. BULB. BOUD 



l$RUN TOOLBOX 



Figure 8. DATA ENTRY Screen slioung Mirror of User Input Values 







Continuing then, we need to modify some other design parameters for our new 
design, so we'll move on to Line 2, design LENGTH (LPP). For this entry we'll enter 
a 2 for Line 2, then the value 1 143.00 that corresponds to our design length. 

Almost as soon as this operation is complete, a HELP WINDOW appears be- 
cause we have violated a range requirement of one of the two sets of calculations. 
Which specific equation we v'iolated is unknown at this point, but the HELP WINDOW 
tells us that the LENGTH/BEAM RATIO IS OUT OF RANGE for one of the calcu- 
lations, as shown in Figure 9 on page 26. This screen is a timed display screen of the 
type shown in Chapter 2 previously, so it will disappear after approximately 15 seconds. 
The user can keep this display on screen if desired by depressing the VT220 key FI which 
corresponds to HOLD SCREEN for the VAX system. Then the screen message can be 
copied or noted as to what problems might have arisen. 

Also notice that in Figure 9 on page 26 that the Results Block on the lower 
right corner of the display has been blanked to avoid confusion while this display is in 
view. After this display has been in view for 15 seconds, it will disappear from view and 
the algorithm will proceed. Upon completion of the calculations, the results will reap- 
pear and the screen will be updated to include our changes to Line 2. This is shown in 
Figure 10 on page 27. Notice in this case that the Silverleaf and Dawson result has been 
set to zero, since we exceeded the algorithms limits on Length-to-Beam ratio as deline- 
ated in Equation 8. The value we entered is perfectly reasonable; however, we made 
such a drastic change in one parameter that is linked to another parameter through these 
range limitations that we exceeded the range for — given by Equation 8. This display 
( Figure 10) now shows the value for the Admiralty Prediction, so we can easily guess 
that our data was invalid for the Silverleaf and Dawson Prediction. The Admiralty value 
is valid, we just have no second value to correlate it against as desired. Notice also in 
Figure 10 on page 27 that the Displacement value of Line 5 has been updated using 
Equation 1 1 , assuming a constant Q value. 

Proceeding further towards our goal of entering the TANKER B design data, 
we need to update the next item, Line 3, DESIGN DRAFT. Thus we'll enter a 3 for 
Line 3, then 74.00 for the TANKER B draft. Again, almost immediately, another HELP 
WINDOW appears, this time displaying a message that we have entered data that makes 
the BEA.M, DIUAFT RATIO OUT OF RANGE. Again, note that the Results Block has 
been cleared since our design data has changed, making the previous calculations invalid. 
We see this screen display in Figure 1 1 on page 28. This display shows that the same 
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mmm i>owi:r rredzctxon paxlure ••• 
TOOL BOX MAY FAIL TO ACCURATELY 
PREDICT POWER REQUIREMENTS FOR 
ADMIRALTY OR SAD METHOD DUE TO 
LENCTH/BCAM RATIO OUT OF RANGE 

PLEASE WAIT 

PROGRAM WILL RESUME 



PLOT GPAPH TO DISK 



RETURN TO MfilN 



EXIT THE PROGRAH 



P ESPOUSE hPEh 



INS TPUCT IONS 



ENTER LENGTH (LPP) IN FT 
[RETURN] TO EXIT 



9. NUrvBER OF SHRFT3 1.0 

Qf REQD INPUTS******** ♦****#**»**# 



***RE5ULTS*** 

******H0R3EP0UiJER. (HP) ********** 

METHOD: 

RDMIRRLTY 

3.&.D.W/.N0 00 

3.&.D.W/. . . .00 

00=.IjUELL.DE3IGNED 0UL0 0OUI . . . 



ISRUN TOOLBOX 



Figure 9. DATA ENTRY Help Screen Display (1) 

routine is used for all the error messages generated in the calculations section, making 
it an efTicient subroutine since it can display context sensitive help messages. 
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rautR iHirnicTioN uinddu 



OPTIONS f’EHU 



INPUT DATA FROM FILE 



STORE DISPLAYED DATA 



plot graph to screen 



PLOT GRAPH TO DISK 



RETURN TO MAIN 



EXIT THE PROGRAM 



RESPONSE AREA 



INSTRUCTIONS 



ENTER A LINE NUMBER 
OR [RETURN] TO EXIT 



DhTh DISP^hY HPnH 



1 


SHIP NAME 




BTANKER 


2 


LENGTH (Lpp) 


, , , (FT) . . 


1143. 00 


3. 


DESIGN DRAFT 


. . , <FT) , . 


35.0 


4. 


BEAM 


. . . (FT) . . 


82. 0 


5. 


DISPLACEMENT 


. . (LT0N3) 


6277B, 


6. 


BLOCK COEFICIENT, . . 


. . , (Cb) . . 


0.669B0 


7. 


SHIP SPEED 


. . . (KT5) , 


20. BO 


B. 


PROPELLER SPEED 


. . . (RPM) . 


120. 0 


9. 


NUMBER OF SHAFTS . . 
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»END OF REQD INPUTS******** 






***RESULT5*** 

METHOD: 

ADMIRALTY 31914, 

S,&,D.IjU/.NO BB O.OOOOOE^OO 

5.&.D.UJ/. . . .BB O.OOOOOE^OO 

BB= WELL DESIGNED BULB BOU) . . . 



ISRUN TOOLBOX 



Figure 10. DATA ENTRY Display after Update 
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HLLP WINnOW 






« « • POWER PREDICTION FAILURE • • « 

TOOL BOX MAY FAIL TO ACCURATELY 
PREDICT POWER REQUIREMENTS FOR 
ADMIRALTY OR SAD METHOD DUE TO 
REAM/DRAFT RATIO OUT OF RANCE 

PLEASE WAIT 

PROGRAM WILL RESUME 



♦ 06 



PLOT GRPPM TO DlSf^ 



RETURM TO r»WIN 



EXIT THE PROGRAM 



RESPONSE ftp Eh 






INSTRUCTIONS 



ENTER DRAFT (T) IN FT. 
[RETURN] TO EXIT 



• V 



U 'F RU F LILLR J ^ VLU ■■ ■ ■ ■' ■ ■f' lMH I V. 

9. NUf^ER OF SHAFTS 1.0 

*»***END OF REQD INPUTS******** ***' 



***RESULT5*** 

******H0RSEP0UIER. (HP) ******* 

METHOD: 

ADMIRALTY 

S.&.D.UU/.NO BB 

S.6.D.UU/. . . .BB 

BB: . WELL . DESIGNED BULB. BOW 



I SPUN TOOLBOX 



Figure 11. DATA ENTRY Help Screen Display (2) 
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After another 15 second delay, the HELP WINDOW in Figure 11 disappears, 
and we can continue. Proceeding along this route further, we must also change the data 
for Lines 4, 5, and 7 to make all the neccessary changes to get to the TANKER B design 
parameters. Notice here that we will not change Line 6 for the Block Coefficient C^, 
since the algorithm knows that we have uniquely specified Draft, Beam, Length, and 
Displacement, thus from Equation 1 1 again, the Block Coefficient is directly specified. 
At any time if we change either the Beam, Draft, or Length, then the code recalculates 
Displacement assuming Cg does not change. But if we force a new Displacement, then 
the code calculates a new based upon Equation 1 1 and the specified Length, Beam, 
and Draft. Finally, all the DATA is updated and our DATA DISPLAY AREA appears 
similar to Figure 12 on page 30, with all the required data in place and new calculations 
for Predicted Power appear in the Results Block. 

At this point we can enter a simple Return key to exit from KEYBO.A.RD 
ENTRY mode, and the screen appears as in Figure 12 on page 30. The mouse pointer 
is now enabled again, and we can proceed through the menu options. Here it would be 
prudent to save our new data set for later use or retrieval, so we can select Option 3 on 
the menu, STORE DISPLAYED DATA, to a file using the filename we gave on Line 
1. Clicking on that Option will bring up another HELP WINDOW, reininding us to use 
different and unique filenames for our data since saving files under the same name will 
simply overwrite the previous data, rendering it useless in the future. This screen ap- 
pears in Figure 13 on page 31 and uses another timed display dialog box developed in 
Chapter 2 Figure 2 on page 8. 

After the file save operation is complete, the program will inform the user of the 
actual filename saved by writing it to the screen display in the INSTRUCTIONS box 
as shown in Figure 14 on page 32. It is important to note here that the program uses 
a default .DAT filename extension and truncates filenames to an 8.3 format. That is, 
eight (8) letters plus the three (3) letters for the .DAT extension. Names longer than this 
will be truncated to this format. 

The user can now select the fourth option on the menu, PLOT GRAPH TO 
SCREEN to have a CA-DISSPLA plot of speed versus power required plotted to the 
VAX terminal screen. Selection of this option will show a screen display as seen in 
Figure 15 on page 33 after the graph is calculated and plotted by the system. This 
process takes a few seconds, so a message appears in the INSTRUCTIONS BOX for the 
user suggesting that they sit back and relax for a moment. The time delay is not all that 
lengthy, and the screen display will soon appear. 
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pouni PKcntciiiiN uiNuou 






OPTIONS KENU 



input dptp from file 



STORE DISPLAYED DATA 



PLOT graph TO screen 



PLOT graph to disk 



RETURN TO MAIN 



EXIT THE program 



RESPONSE hPEh 



INSTRUCTIONS 



ENTER R LINE NUMBER 
OR [RETURN] TO EXIT 



tP 

^■SPUN TOOLBOX 

WLsmmmmmmm 



DhTh DISP-hY hPEh 



1. 


SHIP NHME 




BTANK.ER 


2. 


LENGTH (Lpp) 


.. . (FT) . . 


1143.00 


3. 


DESIGN DRAFT 


.. . (FT) . . 


74.00 


4. 


BEAM 


.. . (FT) . . 


220.00 


5. 


DISPLACEMENT 


. . (LTONS) 


450910.0 


6. 


BLOCK COEFICIENT. . . 


. . . (Cb) . . 


0.01036 


7. 


SHIP SPEED 


. . . (KTS) . 


15.90 


0. 


PROPELLER SPEED 


. . . (RPM) . 


120.0 


9, 


NUMBER OF SHAFTS. . . 




1.0 



**»**END OF REQD INPUTS* 



***RESULTS*** 

•••••HORSEPOWER. (HP) •**•**• 

METHOD: 

RDMIRRLTV 

S.8.D.W/.N0 00 

S.S.O.W/ 00 

00:. WELL. DESIGNED. BULB. BOW. 



53069 

36031. 

34229. 



Figure 12. DATA ENTRY Screen Display >vith all Input Changes 



Here we encounter the first difficult point of the program, deleting the 
CA-DISSPLA plot and continuing our work. Although a number of methods have been 
tried, the easiest method of deleting the plot from the screen is to click the mouse pointer 
on the calling window just visible at the lower edge of the TOOL BOX windows, in order 
to bring this window to the foreground. Now a simple carriage return will delete the plot 
window, and the user can again click the mouse pointer on the TOOL BOX display 
window ( now in the background ) to bring it back to the foreground, and continue. 
This method is archaic, but it is the only consistent method of deleting the display this 
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Figure 13. DATA STORE Display using Timed Dialog Boxes 

author has found. It does not seem entirely intuitive, but it works. A better method 
must be found. 

Alternatively, the user can select the next screen option. Option 5, PLOT 
GRAPH TO DISK to plot the same screen display graph to a disk file. If we select this 
option, again, some time is involved in plotting the CA-DISSPLA graph, so the program 
informs the user of this in the INSTRUCTIONS area as before. 

This graph plotting operation takes approximately 3 minutes to complete with 
no progress updates for the user to follow, so some patience is required. When the plot 
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PDWfR PRCOICTION (tfINDOU 



OPTIONS MENU 



DhTh DISPlhY hPEh 



ENTER DATA FROM KB 



INPUT D«TA from file 



STORE OISPLRYED DRTR 



PLOT GRPPH TO SCREEN 



PLOT GRAPH TO DISK 



RETURN TO MAIN 



EXIT the PROGRAM 



RESPONSE APEh 



INSTRUCTIONS 



FILE 5RVED RS 
BTRNKER .DRT 



1 


SHIP NRME 




BTRNKER 

1143.00 


2 


LENGTH (Lpp) 


. . . (FT) . . 


3. 


DESIGN DRRFT 


. . . (FT) . . 


74.00 


4. 


BERM 


. . . (FT) . . 


22B 00 


5 


DISPLRCEMENT 


. . (LTONS) 


450910.0 


6. 


BLOCK COEFICIENT. . . 


. . . (Cb) . . 


0.B1B36 


7. 


SHIP SPEED 


. . . (K.TS) . 


15.90 


0 . 


PROPELLER SPEED. . . . 


. . . (RPM) . 


120.0 


9, 


NUNBER OF SHRFTS. . . 




1.0 



QP REQD iNpy »*^**P***»»** 



***RE5ULT5*** 

^ ^f^p^ *-*******-** »#*■**■*■**■**■** 

METHOD; 

RDMIRRLTV 

5.&.D.UJ/.N0 BB 

S.&.D.W/. . . .BB 

BB= . WELL . DESIGNED . BULB. BOW 



53069. 

36031. 

34229. 



ISRUN TOOLBOX 



Figure 14. DATA STORE Filename Screen Display 

is complete, the program informs the user via another HELP WINDOW dialog box 
shown in Figure 16 on page 34 and uses a timed dialog box routine similar to those 
developed previously. This window informs the user of the filename that the graph was 
saved under, STD00001.DAT by default. It also informs the user how to get a hardcopy 
printout of the plot using the high resolution of the VAX Laser printer. 

Upon completion of this display, the user can either EXIT in one of the modes 
discussed before or can continue to develop other ship designs using different parame- 
ters. One further note is needed though; CA-DISSPLA by default appends new plots 
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Figure 15. PLOT GRAPH Screen Display for Screen Plot 

onto old existing plot files under the same filename, so any new plot will not have unique 
VAX/ VMS filenames or version numbers. This can be something of a sore spot if you 
only desire the last plot in a large plot file; however, the file can be edited with the EDT 
EDITOR to remove old plots. But a bonus is that new plots will not delete old plots 
since the VAX, 'VMS system defaults to storing only the two most recent version of a file, 
so no data will be lost due to that method, it may just be at the end of a very large file. 
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IE 



HUP UIINOOU 



P ENT 



SAVING DEFINED GRAPH FILE ••• 

GRAPH IS BEING WRITTEN TO DISK 

USING THE FILENAME: 

*'STD00003.. DAT** 

YOU CAN OBTAIN HARDCOPY PLOTS 

BY THE DCL COMMAND: 

**PRINT/DUEUE»LASER STDOOOOX. DAT** 
AFTER TERMINATION OF TOOLBOX. 

PLEASE WAIT 

PROGRAM WILL RESUME 



EXIT THE PPOGPAM 



RESPONSE hPEh 



INSTRUCTIONS 



CftLCULRTING PLOT. . . 
PLEftSL RELAX & WAIT 



♦♦♦RESULTS*** 

METHOD: 

ADMIRALTY 

S.&.D.UU/.NO BB 

S.&.D.U)/ BB 

BB: . WELL. DESIGNED. BULB. BOUU 



53069. 

36031. 

34229. 



jWUN ON 9/5/90 USING SERIAL NUMBER 9132 AT NAUV POSTGRAO SCHOOL 
jpROPRIETARV SOFTWARE PROOUCT OF COMPUTER ASSOCIATES. INC. 

IM132 UIRTUAL STORAGE REFERENCES; 0 REAOS; O WRITES. 



Figure 16 . PLOT GRAPH Screen Display upon Completion of Disk Plot 

Hopefully, we have shown enough of a complete process here to allow a new 
user enough familiarity to accomplish preliminary design tasks using this Module. 
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3, Po«er Prediction Module Output 

There are four (4) primary' outputs of the Power Prediction Module; 

• The screen displays, for quick calculations or information. 

• The saved data file, used alone or in conjunction with the Report Module. 

• The Speed-Power Relationship plots. 

• The Report Module printout of the Power Prediction data, which is much more 
useful than the second item, but requires the use of another Module. 

a. SCREEN DISPLAYS 

These displays have been dealt with in some depth, and there is little else to 
say except to use the program, get familiar with it and its capabilities, and use it as a 
design tool. The speed at which new calculations can be performed by simply changing 
any variable make this type of program excellent for the What if..... type of calculations, 
without the boring repetitiveness of hand performed calculations. 

b. STORED DATA FILE 

Figure 17 is an example of a stored data file typically output from the pro- 
gram. In this case, it corresponds to the previous developmental example, the 
BTANKER.DAT file: 



BTANKER 

1143. 0 

74. 000 

228. 00 

0. 45091E-f06 
0. 81836 
15. 900 

120. 00 
1 . 0000 
53069. 
36031. 
34229. 



Figure 17. STORED DATA File (Typical); Output from Worked Example 



As it is easily seen that, unless the user is familiar with the order of the re- 
sult variables, very' little information is easily obtained from Figure 17, although all the 
information is there in the same order as the screen displays dealt with previously. 
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f. SPEED-POWER RELA TIONSHIP Plots 

Figure IS on page 37 is also from the BTAXKER development and is the 
disk plot that resulted from that data. The screen plot is identical to this plot. 

As can be seen from this figure, much useful information is displayed on this 
plot. In the lower righthand corner w'e have both the filename that the plot was devel- 
oped under and the time it was developed. Also, we have a full plot of speed in knots 
versus power in horsepower for the full range of ship speeds from zero (0.0) to the design 
operating speed (PW_VK). It is easily noted that all three plots trend upward in a cubic 
relationship to speed through the water. One interesting point however was the tend- 
ency of the Silverleaf and Dawson plots to start at very high power at low speeds, de- 
crease to a minimum, then trend upw'ards again similar to the Admiralty plots. This 
error can be ignored since our equations show a definite relationship of speed at or near 
the Boundary Speed Vg from Equation 13. At low speeds, the ship speed is far aw'ay 
from the required proximity to the Boundary Speed, so the prediction may not be en- 
tirely valid. That plus the speed-to-length ratio at low speed may preclude the ships 
power from even being predicted at these extremely low speeds. The solution then was 
to find the minimum Predicted Power level and its corresponding speed in the power 
matrix and set all Predicted Potver points at speeds slower than that speed equal to the 
minimum Predicted Power. This has the effect of saying to the user that some speed 
exists at which efficient operation will result in minimum power useage at speeds below 
that particular speed. At operations above that speed, the penalty is seen in the cubic 
growth of the power required curve. This is a side benefit of using this Predicted Power 
curve Relationship that would not have been obvious without the plots. 
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Figure 18. SPEED-POWER RELATIONSHIP Plot 
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d. REPORT MODULE Predicted Power Report 

It was also desired to use the existing Report Module to add an additional 
report using the data generated by the Power Prediction module. This subroutine source 
code is shown in Appendix D and is consistent with the display and data entry format 
of the other Report sections. Again, we will use our BTANKER data file and develop 
an illustrative example for the report generated by this section. 

After the user has finished generating the design data using the Power Pre- 
diction Module, the data file must be saved to user storage and the filename noted for 
future use. Then, exiting the Power Prediction module using the RETURN TO MAIN 
option, the user can now select the MAIN MENUs GENERATE REPORT option us- 
ing the mouse pointer. This selction will produce a new screen display for the Report 
Module as shown in Figure 19 on page 39. This report menu is common to all the 
sections of the Report Module, with the new addition of Option 3, CREATE POWER 
PREDICTION RPT. 

The user can select this option using the mouse device, and the screen dis- 
play will shift to appear as Figure 20 on page 40, with the users selection highlighted in 
reverse video and the program query filling the right hand side. Here we will enter 
BTANKER.DAT in response to the programs request for the data filename and exten- 
sion to be used in the report. Notice that what we type is also echoed to the screen 
display as we type it. 

Upon entering the data filename as requested, the right hand side of the 
screen display will clear, and again the user will be queried for another entry, this time 
for the filename to be given to the produced report, as shown in Figure 21 on page 41 
For this we can give any name desired, however a logical choice would be the same 
filename as the data with a .RPT or .PRN extension to indicate what it contains. We 
have chosen to use BTANKER.PRN, although any other selection would be adequate. 

• Again, the user is cautioned not to use the same name and extension as another 
pre-existing file (unless that data is to be discarded) since the VAX System will only 
save two copies of a file with the same name and extension. Thus it is possible for 
old data files to be deleted by the system if the user is not careful. 

Upon entering the output filename and extension as requested in 
Figure 21, the right hand side will again clear, and the report will be written to disk 
storage using the filename specified by the user. At this point, the user can enter EXIT 
and be returned to the MAIN MENU screen, or the Report we have just generated can 
be sent to an output device. This selection is option 5, PRINT A REPORT, and is 



38 




[ 




CREATE MANEUVERING RERORI 



CREATE STATIC STAB. REPORT 



CREATE ROWER RREDICTION RRT 



CREATE ENDURANCE RERORT 



RRINT A REPORT 



EXIT 



EXIT THE PROGRAM 



w — 

^■^RUN TOOLBOX 



Figure 19. REPORT MODULE Opening Screen Display 

shown in Figure 22 on page 42 after the selction has been made. Here the program will 
again query the user for the filename and extension of the report to be printed. Any 
existing filename and extension can be used, but it would be prudent to print out the 
report that we have just finished generating. So here we'll enter BTANKER.PRN as 
requested, which is shown in Figure 22, to complete the action. 

Finally, the program needs to know what output device to be used to print 
the report on, and it will queiy’ the user for this information as shown in Figure 23 on 
page 43. At this point we must know the VAX system names for output devices, either 
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RI TOUT MCNII 



a 



CPtATE NANEUVEPlHG PEPOPT 



CPEATE STATIC STAB PEPOPT 



CPEATE ENOUPANCE PEPOPT 



PPIHT A PEPOPT 



EXIT the PPOGPAM 



E 



*RUN TOOLBOX 



ENTER THE FILE NAME 
AND EXTENSION OF THE 
DATA FILE TO BE USED 
OR [RETURN] TO EXIT 



BTANKER.DAT 



Figure 20. REPORT MODULE Request for Data Entr)' (Input) 

LASER, LA210_i (i= 1,2), or LA75_i (i= 1,2), as valid printing devices. For this exam- 
ple, we’ll enter LASER as shown, and the report file will be sent to the named device 
as requested. At this point we can EXIT, or accomplish another print task, if we have 
generated multiple data files from the POWER PREDICTION Module. In this manner, 
multiple reports can be created and printed in one session, so long as unique files contain 
the input data required. 
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CREATE liANEUVERING REPORT 



CREATE STATIC STAB. REPORT 




CREATE ENDURANCE REPORT 



PRINT A REPORT 



EXIT 



ENTER THE FILE NAME 
AND EXTENSION TO BE 
ASSIGNED TO THE REPORT 
OR [RETURN] TO EXIT 



BTANKER.PRN 



EXIT THE program 



— 

*RUN TOOLBOX 



Figure 21. REPORT MODULE Request for Data Entr}’ (Output) 
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CREATE WNEUVEPING PEROPT 



CREATE STATIC STAB. REPORT 



CREATE POWER PREDICTION RPT 



create endurance report 




EXIT 



ENTER THE FILE NAME 
AND EXTENSION OF THE 
REPORT TO BE PRINTED 
OR [RETURN] TO EXIT 



BTANKER.PRN 



EXIT the program 



$RUN TOOLBOX 



Figure 22. 



REPORT MODULE Report Output File Request 



fill [IE 



Ilf PORT MI Nil 



CREATE MANEUVERING REPORT 



CREATE STATIC STAB REPORT 



CREATE ROWER PREDICTION RRT 



CREATE ENDURANCE REPORT 



EXIT THE program 



I 



SRUN TOOLBOX 



ENTER THE PRINTER TO 
USED TO PRINT THE FILE 



LASER 



Figure 23. REPORT MODULE Output Device Request 
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{]) PREDICTED POWER Output Report. Shown in Figure 24 on page 
45 and Figure 25 on page 46 is the Report file created by the Report Module, using the 
data w'e created from the BTANKER development. As can be seen, some other useful 
information is given in the Report that was not available in either the stored data file 
or via the screen displays. 

We have added a calculation for some of the Coefficients of Form 
that are so readily used in preliminary ship design at various points, and are always 
helpful to have handy. We have made them available here in the Report since they 
would not comfortably fit on the Screen Displays shown previously. 
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TOOL BOX 

POWER PREDICTION REPORT 
THIS REPORT WAS GENERATED USING THE PROGRAM 
TOOL BOX WHICH WAS DEVELOPED FOR THE NAVAL 
ENGINEERING DEPARTMENT OF THE 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 

PROFESSOR F. PAPOULIAS 
AND 

LT. GERALD MCGOWAN 
LT. JAMES PLOSAY 
1989/90 



PAGE 1 OF 2 / POWER PREDICTIOI^ 



Figure 24. REPORT MODULE Output Report, Page 1 of 2 



POWER PREDICTION REPORT 

THE INPUT SHIP PARAMETERS ARE AS FOLLOWS 

THE REPORT IS' LOCATED IN* FILE : BTANKER.PRN 

THE INPUT DATA FILE USED IS ; BTANKER.DAT 

************************************************** 

DESIGN PARAMETERS 

SHIP NAME IS BTANKER 



LENGTH BETWEEN PERPENDICULARS FT 


1143.00 


DESIGN DRAFT FT 


74.0000 


DESIGN BEAM FT 


228.000 


DESIGN DISPLACEMENT LTONS 


450910. 


DESIGN DISPLACEMENT LTONS 

COEFFICIENTS OF FORM******************** 

BLOCK COEFFICIENT: 0.818360 

SPEED-LENGTH RATIO: 0.470299 

BEAM-DRAFT RATIO; 3.08108 

LENGTH-BEAM RATIO: 5.01316 

***«***«*«*****«*****«***********«*«***« 


450910. 


DESIGN OPERATING SPEED (NOM) KNOTS 


15.9000 


PROPULSION SHAFT SPEED RPM 


120.000 


NUMBER OF PROPULSION SHAFTS 

************************************************** 

DESIGN RESULTS (HORSEPOWER) 


1.00000 


ADMIRALTY POWER 

SILVERLEAF & DAWSON POWER: 


53068.7 


DESIGN W/0 BULBOUS BOW 


36030.9 


DESIGN W/ BULBOUS BOW 


34229.3 


A************************************************* 

POWER ESTIMATIONS PRESENTED ARE BASED UPON FOLLOWING: 

(1) METHOD OF ADMIRALTY COEFFICTFMT- 

HARVALD, Sv.Aa . , "RESISTAPCF "1' ri -r... - i t 

SHIPS", JOHN WILEY & SOUS, yO»^' ' '0? 

(2) METHOD OF SILVERLEAF AND DM’-<:^il, AS MOD’ TIED BY 
ERICHSEN: ERICHSEN,S., REPORT Mo. 123. 

"OPTIMUM CAPACITY OF SHIPS A5JD ^OFT TEPnillALS" 
UNIV. OF MICHIGAN, ANN AP.P^^.M’.; ’?"i 

***************************♦-■-■#-•-■-■-- --^--**i»***-f 

PAGE 2 OF 2 / POWER PREDICTION 


S . 



Figure 25. REPORT MODULE Output report, Page 2 of 2 
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E. CONCLUSIONS 

Thus we conlude the development and examples for the Predicted Power calcu- 
lations section. Using the data displayed in Table 1 on page 17 we can conclude that 
our methods seem reasonably accurate enough for preliminary design work, and cer- 
tainly when both values are available. 

It would be possible to use this Module and the Report with only one of the three 
Power predictions available by the Module, either just the Admiralty Prediction, or the 
two Silverleaf and Dawson Predictions. But this gives us no real feel for the accuracy 
of the predicted value by cross checking it against another value from an independently 
developed method. But using these figures and some additional data from reference 
material [Ref 4 : p. 127, Ref 1 1 : p. 16], we can begin to get the preliminary data for 
Machinery plant size, manning requirements, costs, fuel storage requirements, and other 
vital pieces of the puzzle that we need to continue the iteration and design the vessel. 
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IV. ENDURANCE CALCULATIONS 



A. BACKGROUND 

Another useful calculation to perform during Preliminary Ship Design once the ex- 
pected Propulsive Power requirements are known, is the Voyage Endurance calculation. 
From this calculation, a number of useful figures arise: expected range of the vessel at 
design operating speed, "the amount of fuel required to be carried for a certain length 
voyage or speed, the added weight of that fuel to the vessels total displacement, and so 
on. These figures will give us the total amount of fuel that needs to be carried, and thus 
the tank volume required to store it, the additional dimensions and size required to 
contain that storage, the relative efficiency of using one type propulsion plant versus 
another, and so forth. Additionally, this type of programmed calculation is well suited 
to the fV/iat if.... type of iterations w'here design variables are modified to see their 
overall effect upon the design. Specifically, items such as: 

• IVhai if we speed up our transit from xx to yy knots? How will that effect the amount 
of fuel vvc have to carry? 

• What if we change from geared diesels to steam turbine propulsion equipment? What 
effect will that change in fuel economy have on the amount of fuel we are required to 
carry? 

As we can see here, this application is well suited to the applications of modern GUI 
programming that can reduce or eliminate the need for repetitious hand calculations, 
thus making it relatively simple to modify a parameter of the preliminarx' design at any 
point and check its effect upon the overall design. 

B. METHODOLOGY 

Basically, we can proceed along one of several directions: (1) we can set the design 
speed and voyage length, then see how much fuel would be required, or (2) we can set 
the amount of fuel carried and how fast we want to go, then see how far we can go, or 
(3) we can set the voyage distance and amount of fuel available, then see how fast we 
can go to get there. The first option, (1), seems to be the most likely candidate for cal- 
culation since it provides the most important factor in this stage of design: the total re- 
quirement of fuel in weight and volume that we are required to carry. The other options 
can be achieved through an iterative process by adjusting the voyage length or speed to 
get a desired weight of required fuel stowage or fixed volumetric limitations. 
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From several sources [Ref. 6 : pp. 173- 174, Ref 12 ; pp. 253- 282], there are al- 
gorithms for calculating endurance based upon all the relevant ship data that should be 
available up to this point in the preliminary design. The essential core parts of the al- 
gorithms are the same however, so we will only examine the most thorough and involved 
method researched. 

I. Method Of U.S. Navy, Design Data Sheets (DDS 9400-1) 

This method, [Ref 13], utilizes a tabular form provided by the Navy to its sur- 
face combatants to calculate all the essential parameters that we require here for en- 
durance calculations. This method and its form , which is shown in Figure 26 on page 
50, is dated from 01 November 1963, but is still valid in all its calculations. 

In order to complete the calculation, the following sets of data must be known: 

• Endurance Variables: 

(1) Endurance Required (miles) 

(2) Endurance Speed (knots) 

(3) Endurance Fuel Load (tons) 

• Design Parameters: 

(4) Full Load Displacement (Itons) 

(5) Rated Full Power (shp) 

(6) Design Endurance Power (shp; at 1,2, and 3) 

(7) Cruising Electric Load (kW) 

(8) Tailpipe Allowance Factor (%) 

• Engineering Parameters: 

(9) Propulsion Fuel Rate (Ibs/shp- hr) 

(10) Auxiliaiy Generator Fuel Rate (Ibs/kw- hr) 

(11) Other Fuel Consumption (Ibs/hr) 

Then, from these sets of numbers, we can proceed to calculate the estimated 
endurance in terms of the amount of fuel we are required to carry, in the following 
manner (using (a 2 ) as referring to the data parameter required from the above list): 

1. First we need the Average Endurance Power, estimated as (6) x 1.10. We will call 
this result (A). Then we can find the ratio of (A) to the rated Full Power, 
{A) j (5) and call this result (B). 

2. At (A), we also need to know the Propulsion Fuel Rate in (Ibs/shp- hr). We will 
call this item (C). 



49 



SECUTUTY CLASanCATJON _ 

(Iflcludia^ Gft>uD No.) 

DATE 

APPENDIX a 

SURFACE SHIP ENDURANCE CALCULATION FORM 

DESIGN ^ 

PREPARED 5Y 

CriZCXZD BY 



EXAMPLES 







S<earL 


1 Diesel Of Gas Tu.Tiine 


(1) 


Eaduri.nce Retired, .Miles 


3,000 


1.200 


(2) 


Eadurincc Speed, lOtoti 


15 


6 


(3) 


Full Load DispUccraeoc, Ton* 


3,000 


400 


(A) 


Rated Full Power, SHP 


50,000 


700 


(5) 


Design Endurance Power J (2)&C3),SHP 


3,000 


150 


(6) 


ATerage Eodurance Power, SHP: 


3.000 X 1.10 - 3,300 


150x I.IO- 165 




(5)x 1.10 






(7) 


Ratio, Avge. End. SHP/r*:cd F.P. SHP: 


.066 


0.24 




(6)/(4) 






(8) 


Cuising Elecric Load, )CT 


500 


30 


(9) 


CalraJaied P.^r^puiiioo Fuel RateP (6), 


— 


0.505 




Ibs/SHP-hr. ’ 






(20) 


Cilc. Prop. Fuel ConiuspdoD, lb»/nr: 


— 


C.505 X 165 » 83.4 




(9)x(6) , 






(11) 


Caic. Aui. Gen. Fuel Rare (S), 




0.690 




Ibs/XT-hr. 






(12) 


Caic. Aui. Gen. Fuel ConsurapdoD, 


— 


0.690 X 30 = 20.8 




Ibs/nr: (ll)x(S) 






(15) 


Caic. Fuel Consumpdon For Odner 


— 


15.0 (heating) 




Serrices, Ibs/hr. 






(U) 


Total Calc, all-punoose Fuel Consurtp- 


— 


83.4 X 20.8 - 15.0 = 119.2 




don. IbsAr: (10) - (12)- (13) 




119.2/165 « 0.722 


(15) 


Calc. All-purpose Fuel Rate, lbs/SH?-br: 


1.00 




(14)/(6) or Heat Balance 




1.04 


(16)* 


Fuel Race Correcdon Facrof Based on(7) 


1.04 


(17) 


Soecitied Fuel Rate, Ibs/SHP-hr; 


1.00 x 1.04 - 1.04 


C.722 X 1.04 = 0.750 




’(15)x (16) 




0.750 X 1.05 » 0.7S7 


(IS) 


Avge. Endurance Fuel Rare, lbs/SHP*hr; 


1.04 X 1.05 - 1.092 




(17)x 1.05 






(19) 


Endurance Fuel (Burnable). Tons: 
(l)x «5)x(18)/(2)x 2240 


3000 X 3300 X 1.C92 


12CC X 165 X 0.7S7 _ ^ 


1 5 x 2240 


6 X 2240 

0.95 

11.6/0.95 = 19-2 


(20) 


Tailpipe Allowance Factor 


0.9S 


CD 


Endurance .^uel Load, tons: (lP)/(20) 


322/0.98 - 329 



REFERENCES FOR SOURCE OF DATA 

Design EndurijiCe Power 

All-Purpose Fuel Rue 

Instilled Fuel Lo»d 



Figure 26. U.S.N. Design Data Sheet (typical), for Endurance Calculations 
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3. Using (C) we can estimate the Propulsion Fuel Consumption, (C) x (A), in (lbs, 'hr), 
and label this result (D). 

4. Also, we will need to know the Auxiliary (Electric) Generator Fuel Rate in 
(Ibs/kW-hr), and the 'Other' Fuel Rate in (Ibs'hr). We will label these as (E) and 
(F). The information required to make these calculations should become available 
after the preliminary Propulsion requirements have been estimated, and a plant 
type is selected from the available literature. 

5. Using (E), we can estimate the Generator Fuel Consumption as (£) x (7), w’^hich 
w'e will label as (G), and then (Z)) + (G) + (F) w^hich we will call the Total Fuel 
Consumption, (H). 

6. N'ow we need a Fuel Consumption Correction, based upon the magnitude of (B). 
This parameter will be labeled (I). 

7. Using (I), W'e can get the Specified Fuel Rate (/) x (A), and label this as (J). 

8. Finally, we can obtain the Average Endurance Fuel Rate, (f) x 1.05, labeled (K). 

9. This figure will give rise to the Theoretical Endurance Fuel required, or 



(1) X (A) X (K) 
(2) X 2240.0 



which we will label as (L). 

10. Before we can finish the calculation however, w'e must account for (8), the Tailpipe 
Allowance that indicates how' much fuel is carried that is below' the suction point 
of the fuel oil transfer system and is thus unaccessible for use. This varies with tank 
geometry, being 0.95 for broad shallow' tanks, and 0.98 for tall narrow' tanks. 

11. Using (8), W'e can calculate the Endurance Fuel Load as (L) / (8) and save this as 
the desired result. 

From the development above, we can easily have the user input the required 
values (1) through (11), and have a simple routine to calculate required fuel load in tons, 
as well as total trip time at the input values of speed and distance, volumetric storage 
requirements for the fuel using standard conversion factors, and total percent of full load 
devoted to fuel storage. These last parameters are useful for checking the validity of the 
design and whether or not it is reasonable in terms of standard design practices. 
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C. COMPARISON OF THE METHOD 



Using the same Design Data as in Chapter 3 from Reference 4 that was used to 
generate Table I on page 17, we can compare our estimated Endurance calculations 
against the actual Endurance Fuel Oil load values of some modern commercial vessels 
at their design endurance load and speed. 



Table 2. COMPARISON OF DDS9400-1 ENDURANCE ESTIMATION 



Ship Design 


' Actual Endurance 
(Tons Fuel Oil) 


DDS9400-1 Endurance 
(Tons Fuel Oil) 


Large General Cargo 


3596.00 


3200.8 


Container 'A' 


3380.00 


2435.1 


Container 'B' 


6943.00 


4032.4 


Roll On- Roll Off 


3465.00 


4617.7 


LASH Barge Carrier 


4928.00 


4088.7 


SE.ABEE Barge Carrier 


5997.00 


5169.9 


Tanker 'A' 


3624.00 


2662.9 


Tanker 'B' 


17857.00 


8427.1 


LNG Tanker 


4400.00 


6040.7 


Bulk Carrier 


1868.00 


2115.0 


Ore Bulk Oil Carrier 


4845.00 


4024.2 



The values used here in comparison are for a general voyage of 10,000.0 nautical 
miles at rated power, using fuel economy values from the literature [Ref. 11 : p.l6], and 
an estimated cruising electric load given by KW — Q.Q\5shp + 1.6iV+ 9jN + 80.0, where 
shp is cruising shaft horsepower and N is the number of crew and passengers onboard, 
with the resultant electic load in Kilowatts (k\V).l All other data is obtained from Ref- 
erence 4, pages 137- 171, and the calculations were performed using TOOL BOX screen 
displays.2 



1 Rule-of-Thumb equation for Preliminary' Electic Load calculations obtained from Naval Sea 
Systems Command, Washington D.C. 

2 Other standards used for unknown values include: Tailpipe Allowance of 0.95, Auxiliary 
generator fuel oil rate of 0.700 Ib/kW-hr, and 'other' fuel oil rate of 15.0 Ib/hr. 
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D. THE TOOL BOX ENDURANCE ESTIMATION MODULE 
1. Operation of the Endurance Estimation Module 

The major operations of the ENDURANCE Module function the same as all 
the other modules. A quick examination of Figure 27 on page 54, will show that the 
first three mouse selection functions: ENTER DATA FROM KB, READ DATA FROM 
FILE, and STORE DISPLAYED DATA are all consistent functions from the other 
modules. The fourth selection, CALCULATE ENDURANCE LIMIT, performs the ac- 
tual calculation of fuel load for the data as it is shown on the program screen. Notice 
from the screen display shown in Figure 27 that the last two required data points, hnes 
10. and 11. are highlighted by an asterisk. These two data points have the calculational 
routine tied in to them so that as this data is modified, the calculations are performed 
and the screen updated immediately. The reasoning for this is that these are the signif- 
icant variables that will be modified most often. Also, these values will most likely not 
be entered until all the other variables are set to non- zero values, saving the time in- 
volved in tracking through the HELP WINDOW error displays during the program in- 
itialization. Otherwise, if a higher line number is modified, the user will need to enter 
the new value, hit RETURN to exit keyboard data entry' mode, and click the mouse 
pointer on the fourth selection, CALCULATE ENDURANCE LIMIT. The function- 
ing methodology, data entry process, and HELP windows function similarily to the 
POWER ESTIMATION Module that we have just covered in some depth, so they will 
not be repeated here. 

The other new feature that has been added to this module is the fifth mouse 
selection shown in Figure 27, that is labeled CHANGE HELP LEVEL ON/OFF. This 
new feature has been added because of the limitations imposed on the programmer due 
to having only a small area available for information presentation to the user in the 
INSTRUCTIONS area at the lower left corner. This small area limits the available 
number of lines of text information that can be presented to the user to provide hints 
or examples of input data that is required. We can easily institute context sensitive 
HELP WINDOWS using our previously developed dialog box routines, but wanted the 
ability to provide the user the ablity to disable HELP information and prompting for 
each value if they desire once they become proficient. This has been accomplished using 
a simple subroutine named SET_HELP_LEVEL that is called from the ENDURANCE 
main menu, and resets a pointer that is an input value to the SET_HELP windowing 
routine used to create and display the help information. The HELP level setting is ON 
by default each time the ENDURANCE Module is called from the MAIN MENU, so 
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PPI llir TIfiN UINIIUU 



■ 1 


OPTIONS f'ENIJ 


□ 


IihTh IHSP-hY hPih 


□ 



ENTER DATA EROM KB 



INPUT DATA FROM fUE 



STORE DISPLAYED DATA 



CALC ENDURANCE LIMIT 



CHANGE HELP LEVEL ON/OFF 



RETURN TO MAIN 



EXIT THE PROGRAM 



P ESPOUSE hPEh 




SELECT ftN OPTION 
WITH THE MJUSE 



1. 


SHIP NAME 


CONT-A 


2. 


FULL LOAD DISPLACEMENT (TONS) 


22080. 


3. 


FULL RATED PLANT POWER (SHP) 


17500. 


4. 


AVG ENDURANCE POWER (SHP) . . . 


17500. 


5. 


CRUISING ELEC LOAD (KW) . 


468.00 


6. 


TAIL PIPE ALLOWANCE (ST) . . 


0.95000 


7. 


CRUISING MAIN FUEL RATE 


0.48000 


0. 


CRUISING ELEC FUEL RATE 


0.70000 


9. 


CRUISING OTHER FUEL RATE 


15.000 


10. 


•EST. ENDURANCE RANGE (MILES). 


10000. 


11. 


•EST. ENDURANCE SPEED (KNOTS). 


20.000 



• »»***»**»*«»»«H»»*»»«*»*«*»*«»*» •*»*•»»**»»* 

••RESULTS FROM CftLC** 

»•»**»»**•»*»*»«»**»**»**•»»*•»• *****•»**»»* 

EST. ENDURANCE FUEL LOAD (TONS) 2435.1 

ELAPSED TRIP TIME (HRS): 500.00 

FUEL STORAGE REQUIREMENTS (FT *3) 92532. 

FUEL WEIGHT ALLOWANCE CU-FL) ; . . 11.028 



RUN TOOLBOX 



Figure 27. ENDURANCE Module Main Data Display Screen 

an advanced user can simply reset the help level indicator upon entering this module if 
they so desire. Figure 28 on page 56 shows this process and the display window that 
informs the user of this selection. Alternatively, the user can set the help level indicator 
back to ON at any time during the progress of the program, without any loss of conti- 
nuity. For this case, the screen display shown would be exactly similar to Figure 28 
except that the instructions would be reversed. 

• NOTE: Setting the HELP_LEVEL indicator to OFF does not disable the presen- 
tation of HELP WINDOWS containing error messages or information messages. 
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I 








This feature solely pertains to the parameter definition displays presented during 
the process of entering data from the keyboard in this module. 

Thus the user has the ability to determine some of the functionality of their 
work environment by changing the way in w’hich information is presented. Once a de- 
sign begins to develop, there will be less need for the parameter definition windows, so 
this option will allow the user the flexibility of definitions only when they are desired. 
The user could also get a program parameter definition by entering that data line num- 
ber while HELP_LEVEL is ON, without changing the actual value entered. 

E. ENDURANCE ESTIMATION MODULE OUTPUT 

Again, the available output methods are essentially siinilar to those of the other 
modules, except that in this case there is no CA-DISSPLA graph output available. The 
screen displays, and stored data files are all closely related to those presented previously, 
so we won't devote space to reviewing them. 

1. REPORT MODULE Estimated Endurance report 

Examining the printed report available from the REPORT Module, we are pre- 
sented with a two page report detailing all the Endurance variables and the results 
obtainined from the ENDURANCE Module as shown in Figure 29 on page 57, and 
Figure 30 on page 58. The method of obtaining a final printout of the report is identical 
to that used by any of the other modules incorporated into TOOL BOX, and any 
printout device available may be used. 

F. CONCLUSIONS 

Thus we conclude the development and discussion of the Endurance Estimation 
Module. The usefulness of the information calculated here can be very helpful in de- 
termining the size of the storage tanks used for Fuel Oil, the predicted or estimated 
amount of Fuel Oil stowage required in tons, and the percent of full load that this 
amount constitutes. The parameters that we have calculated all become essential ele- 
ments of the iterative design process shown in Figure 1 on page 3, and thus the easier 
we can make the process of obtaining them, then the entire design process will become 
easier. 
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Hiir wi 




• ««» SYSTEM SETTINGS •••• 

CHANGING HELP LEVEL FROM 
HELP WINDOWS ON 
TO 

HELP WINDOWS OFF 

PLEASE WAIT 

PROGRAM WILL RESUME 



RETURN TO MR IN 



EXIT THE RROCRAM 




W LLW F ULL 

9. CRUISING OTHER FUEL RftTE 15.000 

10. *E3T.ENDURftNCE RftNGE (MILES). 10000. 

11. »EST. ENDURANCE SPEED (KNOTS). 20.000 

••RESULTS FROM CALC*» 

»»»***»**»»* 

EST. ENDURANCE FUEL LOAD (TONS) 2435.1 

ELAPSED TRIP TIME (HRS); 500.00 

FUEL STORAGE REQUIREMENTS (FT *3) 92532. 

FUEL UtIGHT ALLOWANCE («-FL) : . . 11.020 



Figure 28. User Display during Resetting of the Help Level 
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Figure 29. 



TOOL BOX 

ENDURANCE LIMIT REPORT 

********************** 

THIS REPORT WAS GENERATED USING THE PROGRAM 
TOOL BOX WHICH WAS DEVELOPED FOR THE NAVAL 
ENGINEERING DEPARTMENT OF THE 
NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 

PROFESSOR F. PAPOULIAS 
AND 

LT. GERALD MCGOWAN 
LT, JAMES PLOSAY 
1989/90 



PAGE 1 OF 2 / ENDURANCE LIMIT 



Endurance Estimation Report, page 1 of 2 
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ENDURANCE LIMIT REPORT 



THE INPUT SHIP PARAMETERS ARE AS FOLLOWS 



THE REPORT IS LOCATED IN FILE : dtanker.prn 

THE INPUT DATA FILE USED IS : dtanker.dat 

************************************************** 

DESIGN PARAMETERS 

SHIP NAME IS dtanker 

DESIGN FULL LOAD DISPLACEMENT: 22080.0 

(LTONS) 

DESIGN FULL POWER LEVEL: 17500.0 

(SHP) 

DESIGN ENDURANCE POWER LEVEL: 16000.0 

(SHP) 

DESIGN CRUISING ELECTRIC LOAD; 500.000 

(KW) 

DESIGN FUEL TANKS TAIL PIPE ALLOWANCE: 0.950000 

(%) 

DESIGN MAIN PROPULSION FUEL ECONOMY: 0.450000 

( LES/SHP-HR) 

DESIGN ELECTRIC PLANT FUEL USE RATE: 0.250000E-01 

( LBS/KW-HR) 

DESIGN OTHER FUEL USE RATES:..... O.lOOOOOE-01 

( LBS/HR) 

DESIGN RANGE: 3000.00 

(NMILES) 

DESIGN CRUISING SPEED: 21.9000 

( KNOTS ) 

CALCULATED ENDURANCE FUEL LOAD: 552.260 

(TONS-FUEL OIL) 

CALCULATED JOURNEY TIME ALLOWANCE: 1 36.990 

(HRS ^ DESIGN SPD & DISTANCE) 

CALCULATED FUEL STORAGE REQUIREMENTS: 209B6.0 

(FT"3 VOLUMETRIC STORAGE) 

CALCULATED FUEL RATIO OF FULL LOAD: 2.50120 

U OF DESIGN FULL LOAD) 



************ltltltltit**1tit**********li***it**itit****i,it*it** 

ENDURANCE LIMIT VALUES PRESENTED ARE BASED UPON 
FOLLOWING: 

(1) U.S. NAVY DESIGN DATA SHEET DDS9400-1 FORMAT 
CALCULATIONS USING ESTIMATION FACTORS FOF 
FOULING, MACHINERY INEFFICIENCIES, ETC 
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Figure 30 



Endurance Estimation Report, page 2 of 2 



V. CONCLUSIONS AND RECOMMENDATIONS 



A. WINDOWING GUI'S AND PRELIMINARY SHIP DESIGN 

We have seen in the original paper [ Ref. 3] that we could design and program a 
Windowing GUI to accomplish preliminary ship design calculations using the capabili- 
ties of the VAX/VMS User Interface Services (UIS) routines and the FORTRAN pro- 
gramming language. Further, we have shown here that we could use these capabilities 
to provide context sensitive HELP information to the user, and enhance the user inter- 
face by the addition of Dialog boxes that provide information at appropriate moments 
during the execution of the program, thus negating a requirement to have a Users 
Manual handy at all times, or to require specific and lengthy training in the use of the 
program. 

B. CONCLUSIONS 

These capabilities provide us with the ability to develop a powerful and easy to use 
in-house development environment for working with preliminar}’ design problems that 
is easily expanded with new Modules and calculational routines as the needs of the 
Students and Faculty evolve, and that places no limitations on the developer. Specif- 
ically, we can summarize the attributes of TOOL BOX \n the following manner: 

• Easily utilized by design engineers and students, without lengthy training or com- 
plicated users manuals, to provide basic ship design parameters and study the ef- 
fects of changing those parameters in an interactive design environment using 
modern graphical programming techniques. 

• Follows accepted progression along an Iterative Design Spiral commonly in use in 
the Marine Industry today, thus affording the developer a logical progression thru 
the Preliminary Design process. 

• Provides valuable information in a timely fashion to the developer in an easily ac- 
cessed manner, thus avoiding the constraints of current Computer Aided Design 
(CAD) programs that do not utilize modularity. 

• Easily enhanced by the addition or modification of any module. 

• Easily programmed in whatever development language is most familiar or appro- 
priate. Here we have used FORTRAN; the UTS routines however are non-specific, 
and can be utilized by any language supported on the VAX Network, including C 
and PASCAL. 
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C. RECOMMENDATIONS 



Throughout the design and programming of the additions to TOOL BOX that we 
have made, we have kept a fairly lengthy wish list of things that we would do better or 
different given the opportunity. Things that 1 feel need to be accomplished in the next 
version include the following: 

• Continue the progression along the Iterative Design Spiral, continuing to add 
functionality and utility to the core TOOL BOX code. 

• Figure out a method of having more AST functions enabled simulaneously than 
we now use. It has been reported by the original author that too many of these 
functions enabled simultaneously causes a system slow-down and sluggishness. 
There should be a means to get around this capability. 

• Reprogram the opening menu, and perhaps the Module menus, to accomodate a 
wider range of user choices. As it stands currently, the Opening Menu has a single 
RESERVE module left for use. This needs to be increased if we are going to 
accomodate the entire design spiral of Figure 1 on page 3 into TOOL BOX. The 
mouse boxes could be made smaller or grouped into two or more columns to 
accomodate this feature. 

• Another possibility would be to group the Module functions into sub-groups, per- 
haps having a single Main Menu box used to select the two functions we have ac- 
complished here, instead of each Module having its own selection box on the Main 
menu. Then there would need to be sub-menus to select which of the two functions 
was desired from those available to the sub-group. 

• Enhanced use of color would make the system more user friendly and visually at- 
tractive. The HELP windows could easily be made difierent in color from the rest 
of the program, thus distinguishing them from the main program. 

• Find a means of adding a FILE-NOT-FOUXD recovery capability to the READ 
DATA FRO.Vl FILE selections that appear on all of the Module menus, thus 
avoiding program crashes due to a simple spelling mistake in data entry procedures. 

• Find a means, perhaps thru some assembly language routines, to significantly in- 
crease the speed of the program. There are a few places where the system response 
has diminished as the code size grew. This program now accounts for roughly 8000 
lines of FORTRAN source code, a language that has never been known for its 
speed. 

Finally, our last recommendation is that the development of TOOL BOX continue 
along these lines and those already established. This program could become a unique 
and easily adaptable tool for the Total Ship Systems Engineering (TSSE) curriculum at 
NTS Monterey. This program would provide a friendly computer-based development 
environment for the TSSE students and greatly aid in preliminary' design problems, thus 
allowing the students to concentrate on new technologies and methods, rather than 
laboriously computing oft-required, and oft-changing variable parameters such as those 
provided here. 
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Complete copies of the TOOL BOX source code, including parts not covered in this 
document, may be obtained from Professor F.A. Papoulias, Code ME/Pa, Naval Post- 
graduate School, Monterey CA 94943-5000. 



61 



APPENDIX A. TYPICAL DIALOG BOX SOURCE CODE 



A. VISUAL DISPLAY INFORMATION BOX 

This is the source code for the pop-up information display boxes utilized in the 
program. 

SUBROUTINE SHOW_SAVE 

c * 

C THIS SUBROUTINE POPS UP A DIALOG BOX TO INSTRUCT THE USER NOT * 

C TO USE THE SAME DATA FILE NAME AS USED IN A PREVIOUS PROGRAM * 

C SECTION [SINCE IT WILL OVERWRITE THE OLD ONE CAUSING ERRORS * 

C OR LOST DATA] . IT IT CALLED BY THE SUBROUTINE THAT SAVES THE * 

C DATA FILES * 

C * 

Cycyryc*yr*)vypyryp*yryrypypyp***yryc**ypycycycyry-yp*“vyp**yryc’yr^v*yr*yryryr'>vyp**yryfyryr')vyp'}vyipyryfyryr*yr*yr**yr*yf 

c 

IMPLICIT INTEGER(A-Z) 

INCLUDE ' SYS$LIBRARY: UISENTRY' 

INCLUDE ' SYS$LIBRARY: UISUSRDEF ' 

CHARACTER*34 SAVEl/’*** SAVING DEFINED DATA FILE *** ’/ 
CHARACTER*34 SAVE2/' DO NOT USE SAME NAME & EXTENSION ’/ 
CHARACTER*34 SAVE3/'AS DATA FILE FOR OTHER SUB-SECTION'/ 

CHARACTER*34 SAVE4/ ' PLEASE WAIT '/ 

CHARACTER*34 SAVES/’ PROGRAM WILL RESUME'/ 

REAL Y_POSN 

C*****CREATE THE DISPLAY FOR THE DIALOG BOX 

VD_SAVE=UIS$CREATE_DISPLAY(0. 0,0. 0,10. 0,5. 0,20. 0,6. 0) 

C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 

WD_SAVE=UIS$CREATE_WINDOW(VD_SAVE, ’SYS$WORKSTATION' , 'HELP BOX' ) 
C*****COPY ATTRIBUTE BLOCK 'O' AS BLOCK '2l' AND CHANGE THE FONT SIZE 

CALL UIS$SET_F0NT(VD_SAVE,0,21, 'DTABER0R03WK00GG0001UZZZZ02A000' ) 
C*****siGNAL THE USER TO GET THEIR ATTENTION 
CALL UIS$SOUND_BELL( ’SYS$WORKSTATION' ,4) 

C**-'**WRITE THE TEXT INTO THE WINDOW 
Y_P0SN=4. 0 

CALL UIS$TEXT(VD_SAVE,21,SAVE1,0. 2,Y_POSN) 

Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE,21,SAVE2,0. 2,Y_POSN) 

Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE,21,SAVE3,0. 2,Y_P0SN) 

Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE,21,SAVE4,0. 2,Y_POSN) 

Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE, 21, SAVES, 0. 2,Y_POSN) 

C*****SPAWN A NEW PROCESS TO WAIT 15 SECONDS AND COMPLETE PROCESS 
CALL LIB$SPAWN( 'WAIT 00: 00: IS' ) 

CALL UIS$SOUND_BELL( 'SYS$VmKSTATION' ,4) 

CALL UIS$DELETE_DISPLAY(VD_SAVE) 

RETURN 

END 



62 



B. CONVERSATIONAL DIALOG BOX 

This segment of code produces the dialog boxes that require user input to continue. 

SUBROUTINE NOWHERE 



c * 

C THIS SUBROUTINE IS CALLED BY A MOUSE BUTTON SELECTION OF AN * 
C OPTION THAT IS NOT ENABLED YET. IT POPS UP A DIALOG 'HELP' * 
C BOX THAT INSTRUCTS THE USER AND GIVES THEM THE OPTION OF * 
C CONTINUING BACK TO THE PROGRAM AT THE POINT OF INTERRUPTION * 
C OR ENDING THE PROGRAM ENTIRELY. * 
C * 



c 

IMPLICIT INTEGER(A-Z) 

INCLUDE 'SYS$LIBRARY:UISENTRY' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

CHARACTER*34 HINTl/'*** THIS FEATURE NOT AVAILABLE ***'/ 
CHARACTER*34 HINT2/’ TYPE "CONTINUE" TO RESUME '/ 

CHARACTER*34 HINT3/’ OR "EXIT" TO END PROGRAM '/ 

REAL Y_POSN 

C*****cREATE THE DISPLAY WINDOW FOR THE DIALOG BOX 

VD_TEST=UIS$CREATE_DISPLAY(0. 0, 0. 0 , 10. 0 ,5. 0 , 20. 0,5. 0) 

C*****CREATE THE WINDOW TO DISPLAY THE HELP TEXT 

WD_TEST=UIS$CREATE_WIND0W(VD_TEST, ' SYS$WORKSTATION' , 'HELP BOX' ) 
C***,v*SET THE SCREEN FONTS WE DESIRE FROM ATTRIBUTE BLOCK 'O' TO BLOCK 
C*****»21' AND CALL THE FILENAME OF THE FONT 

CALL UIS$SET_FONT(VD_TEST,O,21,'DTABEROR03WK00GG00OlUZZZZ02A000' ) 
C*****siGNAL THE USER ON PROCESS PAUSE FOR DISPLAY OF HELP BOX 
CALL UIS$SOUND_BELL( 'SYS $ WORKSTATION' ,4) 

C*****PRINT THE TEXT AS WE WANT IT 
Y_P0SN=4. 0 

CALL UIS$TEXT(VD_TEST, 21, HINTl, 0. 4,Y_POSN) 

Y_POSN=Y_POSN-l 

CALL UIS$TEXT(VD_TEST,21,HINT2,0. 4,Y_POSN) 

Y_POSN=Y_POSN-l 

CALL UIS$TEXT(VD_TEST,21,HINT3,0.4,Y_POSN) 

PAUSE 

C*****siGNAL THE USER PROCESS RESTARTED 

CALL UIS$SOUND_BELL( 'SYS$WORKSTATION' ,4) 

c*****return to calling process at point of interrupt 

CALL UIS$DELETE_DISPLAY(VD_TEST) 

RETURN 

END 
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APPENDIX B. TOOL BOX POWER PREDICTION 

A. SUBROUTINE POWER SOURCE CODE 

This is the source code for the TOOL BOX Power Prediction module and all the 
subroutines called by the module. The INCLUDE File TOP_POWER.FOR follows this 
section. 



SUBROUTINE POWER 

* 



THIS SUBROUTINE IS THE MAIN PART OF PROPULSIVE POWER CALCULATIONS. * 



IT SETS UP SCREENS AND CONTROLS POWER PROGRAM ACTIONS. 




* 

•k 


IT IS CALLED 


FROM THE MAIN TOOL BOX MENU BY THE "POWER : 


PREDICTION" 


* 


SELECTION USING THE MOUSE. IT REQUIRES THE FOLLOWING INPUTS FOR 


* 


OPERATION AFTER INITIALIZATION: 




* 


PW FNAME 


SHIP NAME GIVEN TO SAVED DATA FILES 


CHAR 


* 


PW LPP 


SHIP LENGTH BE'^^EEN PERPENDICULARS (FT) 


REAL 


* 


PW NSHAFT 


NUMBER OF PROPULSION SHAFTS ( 1 OR 2 ) 


REAL 


* 


PW DISP 


DISPLACEMENT (TONS STD SEAWATER) 


REAL 


* 


PW VK 


SHIP SPEED (KTS) 


REAL 


* 


PW B 


SHIP BEAM (FT) 


REAL 


Vc 


PW T 


SHIP DRAFT (FT) 


REAL 


Vf 


PW CB 


SHIP BLOCK COEFFICIENT 


REAL 


* 


PW NRPM 


SHAFT SPEED (RPM) 


REAL 


* 

k 


AND CALCULATES THE FOLLOWING INTERMEDIATE VALUES: 




k 


PW NRPS 


SHAFT SPEED (RPS) 




k 


PW LPPM 


SHIP LENGTH BETWEEN PERPENDICULARS (METERS) 


k 


PW VKM 


SHIP SPEED (METERS PER SECOND) 




k 


PW DISPM 


SHIP DISPLACEMENT (METRIC TONS) 




k 


PW VB 


S&D BOUNDARY SPEED VALUE 




k 


PW SLR 


S&D SPEED-LENGTH RATIO 




k 


PW BTR 


S&D BEAM -DRAFT RATIO 




k 


PW.LBR 


S&D LENGTH-BEAM RATIO 




k 




AND OTHER CALCULATION-SPECIFIC VARIABLES 




k 

k 


AND OUTPUTS THE FOLLOWING VALUES: 




k 


PW ADME 


ADMIRALTY POWER PREDICTION 




k 


PW SDNB 


S&D PREDICTION FOR NO BULBOUS BOW 




k 


PW_SDBB 


S&D PREDICTION FOR SHIP INCLUDING A BULBOUS BOW 


k 




ALL INPUT DESIGN VARIABLES 




k 


Sc Sc St St 6c St St 


UPDATED VALUES OF INPUTS IF ANY PARAMETER 


CHANGES 


k 

k 


ALL PARAMETERS ARE PASSED IN COMMON BY THE INCLUDE 'TOP. 


_POWER. FOR’ 


k 


FILE THAT 


HAS ALL THE PARAMETER DEFINITIONS, TYPE STATEMENTS, 


k 


AND SET UP. 




k 

k 



CALLED BY MAIN MODULE ’T00L_B0X' 



* 
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C CALLS SUBROUTINES ' POWER_KB_DATA_IN' , ' PW READ_FILE ' PW_RECORD' , * 
C ' PW_CALC_ 1 ' , ' POWER.GRAPH ' , ' TO_MAIN^ , AND ' EXIT ' * 
C * 



INCLUDE 'GENERAL. FOR' 

INCLUDE 'TOP_POWER. FOR' 

EXTERNAL TO_MAIN , EXIT , DARK2 , LIGHT2 .NOWHERE , PW_CALC_2 
EXTERNAL POWER_KB_DATA_IN , PW_RECORD , PW_READ_FILE , PW_CALC_1 
EXTERNAL PW_SDI_LABEL , PW_SORT , PW_FAI L , POWER_GRAPH 
EXTERNAL G_DISK,G_SCREEN 

c*****initialize the parameters to avoid divide by zero errors, also 

C THIS INSURES THAT WHEN WE INITIALLY CALL THE KEYBOARD DATA ENTRY 

C ROUTINE THAT THE DATA PASSES ALL THE 'PW_CALC_2' FAILURE 

C CRITERION. OTHERWISE WE HAVE THE POTENTIAL TO TAKE "FOREVER" 

C WAITING THROUGH THE HELP SCREENS AS WE SLOWLY ENTER OUR DESIGN 

C DATA. ALSO, IT GIVES THE USER AN IDEA OF WHAT KIND OF SHIP 

C WILL GENERATE VALID DATA. 

IF ( PW_TIMES .EQ. 0 ) THEN 
PW_KB_UP( 1)=' GENCARGO' 

PW_INFILE=' GENCARGO ' 

PW_FNAME=' GENCARGO. DAT' 

PW_NSHAFT=1. 0 
PW_KB_UP(9)=' 1. o' 

PW_VK=20. 80 
PW_KB_UP(7)='20. 80' 

PW_LPP=582. 50 
PW_KB_UP(2) = '582. 50 ' 

PW_T=35. 0 
PW_KB_UP(3)='35. O' 

PW_B=82. 0 
PW_KB_UP(4)='82. o' 

PW_DISP=31995. 0 
PW_KB_UP(5)=' 31995. O' 

PW_CB=0. 6698 

PW_KB_UP(6)='0. 66980' ' 

PW_NRPM=120. 0 
PW_KB_UP(8)=' 120. O' 

ENDIF 

PW TIMES=PW TIMES+1 



C*****ERASE THE DATA AREA 



CALL UIS$ERASE(VD_ID,3. 6,10. 1,9. 9,18. 5) 

C*****REMOVE THE MOUSE POINTER TO SOME OTHER AREA 

STA = UIS$SET_P0INTER_P0SITI0N(VD_ID,WD_MAIN,9. 9,4. 5) 
C*****CREATE POWER WINDOW 

WD_POWER=UIS$CREATE_WINDOW(VD_ID, ' SYS$WORKSTATION' , 'POWER 
& PREDICTION WINDOW' ,-.5,9. 9,10. 1,19. ,40. ,30. ) 

DO 10 Y_COOR = 13. 1,18. 4,. 8 
DY = Y_COOR + . 2 

CALL UIS$SET_P0INTER_AST(VD_ID,WD_P0WER,DARK2, ,X0,Y_C00R,X1, 
& DY.LIGHT2) 

10 CONTINUE 
C*****SET MENU TITLES 

OPTION(ll) = 0PTI0N(21) 

0PTI0N(12) = OPTION(20) 

C*****WRITE SCREEN 

CALL UIS$TEXT(VD_ID,0,OPTION(11),. 3,15. 1) 
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CALL UIS$TEXT(VD_ID,0,OPTION(12),. 3,15. 9) 

CALL UIS$SET_BUTTON_AST(VD_ID,WD_POWER,EXIT, , ,X0,13. 0, 

& XI, 13. 4) 

CALL UIS$SET_BUTTON_AST(VD_ID,WD_POWER,TO_MAIN, , ,X0, 13. 8, 

& XI, 14. 2) 

CALL UIS$SET_BUTTON_AST(VD_ID,WD_POWER,G_SCREEN, , ,X0,14. 6, 

& XI, 15.0) 

CALL UIS$SET_BUTT0N_AST(VD_ID,WD_P0WER,G_DISK,,,X0,15.4, ■ 

& XI, 15. 8) 

CALL UIS$SET_BUTTON_AST(VD_ID,WD POWER, PW RECORD, , ,X0, 16. 2, 
& XI, 16. 6) 

CALL UIS$SET_BUTTON_AST(VD ID,WD_POWER,PW READ FILE,,,XO 
& ,17. , XI, 17. 4) 

CALL UI S $ SET_BUTTON_AST( VD_ID , WD_POWER , POWEILKB_DATA_IN , 

& ,,X0,17. 8, XI, 18. 2) 

19 WNDOW = WD.POWER 

20 RETURN 

END 

C 

C 

SUBROUTINE POWER_KB_DATA_IN 

Q*************^V**^V*********************************'A^************** 



c * 

C THE POWER_KB_DATA_IN ROUTINES LOAD IN KB DATA AND ALLOWS * 
C THE OPERATOR TO CHANGE INPUTS AND SEE REAL TIME EFFECTS * 
C OF THOSE CHANGES ON THE POWERING CHARACTERISTICS OF THE * 
C PRELIMINARY DESIGN. * 
C * 
C CALLED BY SUBROUTINE 'POWER' * 
C CALLS SUBROUTINES 'PW_SDI_LABEL' , 'KEY_READ' , 'PW.SORT' , * 
C 'PW_CALC_1' * 



c 

INCLUDE 'TOP. FOR' 

I NC LUDE ' TOP_POWE R. FOR ' 

REAL PW_LINE_NO 

STA=UIS$SET_P0INTER_P0SITI0N(VD_ID,WD_P0WER,9. 9,14. 5) 
KB_ID=UIS$CREATE_KB( ' SYS$W0RKSTATI0N' ) 

CALL UIS$ENABLE_KB(KB_ID,WD_POWER) 

C*****fIRST, we we write UP THE LINE LABLES AND FIRST INSTRUCTIONS 
CALL PW_SDI_LABEL(VD_ID) 

C*****sECOND, CALCULATE AND DISPLAY STARTING VALUES 
CALL PW_CALC_1 

C*****THIRD, ALLOW USER TO CHANGE INPUT THEN 

c*****recalculate results in real time 

1 CALL UIS$ERASE(VD_ID, -0.4, 10. 1,3.4,10.9) 

CALL UIS$TEXT(VD_ID,7,'ENTER A LINE NUMBER' ,0. 0,11. 0) 

CALL UIS$TEXT(VD_ID,7,'0R [RETURN] TO EXIT' ,0. 0 , 10. 6) 

3 CALL KEY_READ (PW_LINE,. TRUE. ,0. 0,12. 0, *300) 

READ (PW_LINE,FMT='(F2. 0)' ,ERR = 5) PW_LINE_NO 
PLINE_NO = INT(PW_LINE_NO) 

GO TO (10,20,30,40,50,60,70,80,90) PLINE_NO 
C 

C*****IF IT GETS TO HERE A MISTAKE HAS BEEN MADE 

c 

5 CALL UIS$ERASE(VD_ID,-0. 4,10. 1,3. 4,10. 9) 
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CALL UIS$TEXT(VD_ID, 7, 'IMPROPER LINE NUMBER’ , 0. 0 11. 0) 

CALL UIS$TEXT(VD_ID,7, 'ENTER A NEW NUMBER PLEASE \ 0. 0 , 10. 6) 
GOTO 3 

CALL PW_SORT( PLINE_NO , DEL_Y, PW_LINE , *300 ) 

PW_KB_UP(1) = PW_LINE 

PW_INFILE = PW_LINE 

READ (PW.LINE, FMT='(G)' ERR = 5) 

WRITE (PW_KB_UP(1), FMT='(A)') PW_INFILE 
CALL UIS$ERASE(VD_ID, 8. 2,17. 9,9. 9,18.3) 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(1),8. 2,18. 2) 

GOTO 1 

CALL PW_SORT( PLINE_NO , DEL_Y , PW_LINE , *300 ) 

PW_KB_UP(2) = PW_LINE 

READ (PW_LINE, FMT=' (G) ’ ,ERR = 5) PW_LPP 
PW.DISP = PW_CB*PW_LPP*PW_B*PW_T/35 
WRITE (PW_KB_UP(5) , FMT=' (G12. 5) ' ) PW_DISP 
CALL UIS$ERASE(VD_ID,8. 2,16. 3,9. 9,16. 7) 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(5),8. 2,16. 6) 

CALL PW_CALC_1 
GOTO 1 

CALL PW_SORT( PLINE_NO , DEL_Y , PW_LINE , *300 ) 

PW_KB_UP(3) = PW_LINE 

READ (PW_LINE, FMT=' (G) ’ ,ERR = 5) PW_T 
PW_DISP = PW_CB*PW_LPP*PW_B*PW_T/35 
WRITE (PW_KB_UP(5), FMT=' (G12. 5) ' ) PW_DISP 
CALL UIS$ERASE(VD_ID,8. 2 , 16. 3 , 9. 9 , 16. 7) 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(5),6. 2,16. 6) 

CALL PW_CALC_1 
GOTO 1 

CALL PW_SORT( PLINE_NO , DEL_Y , PW_LINE , *300 ) 

PW_KB_UP(4) = PW_LINE 

READ (PW_LINE, FMT=' (G) ' ,ERR = 5) PW_B 

PW_DISP = PW_CB*PW_LPP*PW_B*PW_T/35 

WRITE (PW_KB_UP(5) , FMT=' ( G12. 5) ' ) PW_DISP 

CALL UIS$ERASE(VD_ID,8. 2,16. 3,9. 9,16. 7) 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(5) ,8. 2,16. 6) 

CALL PW_CALC_1 
GOTO 1 

CALL PW_SORT( PLINE_NO , DEL_Y , PW.LINE , *300) 

PW_KB_UP(5) = PW_LINE 

READ (PW_LINE, FMT=' (G) ' ,ERR = 5) PW_DISP 
PW_CB=3 5*PW_D I SP/ ( PW_LPP*PW_B*PW_T ) 

WRITE (PW_KB_UP(6), FMT=' (G12. 5) ' ) PW_CB 
CALL UIS$ERASE(VD_ID,8. 2 , 15. 9 , 9. 9 , 16. 3) 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(6),8. 2,16. 2) 

CALL PW_CALC_1 
GOTO 1 

CALL PW_SORT( PLINE_NO , DEL_Y , PW_LINE , *300 ) 

PW_KB_UP(6) = PW_LINE 

READ (PW_LINE, FMT=' (G) ' ,ERR = 5) PW_CB 



nonnnoooo no 



PW_DISP = PW_CB*PW_LPP*PW_B*PW_T/35 
WRITE (PW_KB_UP(5), FMT=* (G12. 5) ' ) PW.DISP 
CALL UIS$ERASE(VD_ID,8. 2,16. 3,9. 9,16. 7) 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(5),8. 2,16. 6) 

CALL PW_CALC_1 
GOTO 1 
C 

70 CALL PW_SORT(PLINE_NO,DEL_Y,PW_LINE,*300) 
PW_KB_UP(7) = PW_LINE 

READ (PW.LINE, FMT='(G)' ERR = 5) PW_VK 
WRITE (PW_KB_UP(7), FMT=‘ (G12. 5) ' ) PW_VK 
CALL PW_CALC_1 
GOTO 1 
C 

80 CALL PW_SORT(PLINE_NO,DEL_Y,PW_LINE,*300) 

PW_KB_UP(8) = PW.LINE 

READ (PW_LINE, FHT='(G)' ERR = 5) PW_NRPM 
WRITE (PW_KB_UP(8), FMT=HG12. 5) ' ) PW_NRPM 
CALL PW_CALC_1 
GOTO 1 
C 

90 CALL PW_SORT(PLINE_NO,DEL_Y,PW_LINE,*300) 

PW_KB_UP(9) = PW_LI.NE 

READ (PW_LINE, FMT=' ( G) ' ,ERR = 5) PW_NSHAFT 
WRITE (PW_KB_UP(9), FMT=' (G12. 5) ' ) PW_NSHAFT 
CALL PW_CALC_1 
GOTO 1 
C 

300 CALL UIS$DISABLE_KB(KB_ID) 

CALL UIS$ERASE(VD_ID,-.4,10. 1,3.4,10.9) 

CALL UIS$TEXT(VD_ID, 7, 'SELECT A.N OPTION', 0. ,11. ) 
CALL UIS$TEXT(VD_ID,7, ' WITH THE MOUSE ',0. ,10.6) 
RETURN 
END 






SUBROUTINE PW_SORT(PLINE_NO,DEL_Y,PW_LINE , COUNT,*) 






SUBROUTINE PW_S0RT SORTS OUT THE LINE NUMBER TO THE 
VARIABLE INVOLVED AND CALLS THE APPROPRIATE SCREEN 
INSTRUCTION TO BE WRITTEN IN THE INSTRUCTION BOX. 



CALLS SUBROUTINE 'KEY READ' 






INCLUDE 'TOP. FOR’ 

INCLUDE 'TOP_POWER. FOR' 

REAL DEL_Y, DY, DOT_Y 
LOGICAL PW_FCHAR 
CHARACTER PW_KB_INST(9)*25 
C*****SET THE SCREEN INSTRUCTIONS STRINGS 
PW_KB_INST(1) = 'ENTER NEW SHIP NAME' 
PW_KB_INST(2) = 'ENTER LENGTH (LPP) IN FT' 
PW_KB_INST(3) = 'ENTER DRAFT (T) IN FT. ' 
PW_KB_INST(4) = 'ENTER BEAM (B) IN FT. ' 



68 



o o 



PW_KB_INST(5) = 'DISPLACEMENT IN LIONS' 

PW_KB_INST(6) = 'ENTER BLOCK COEFICIENT' 

PW_KB_INST(7) = 'ENTER SHIP SPEED IN KNOTS' 
PW_KB_INST(8) = 'ENTER PROPELLER RPM' 

PW_KB_INST(9) = 'ENTER NUMBER OF SHAFTS' 

C*****erasE the instruction area AND FIND ITS DATA POSITION 
CALL UIS$ERASE(VD_ID,-0. 4, 10. 1,3. 4,10. 9) 

DEL_Y = 18. 6 - . 4 * PLINE_NO 
C*****writE the appropriate INSTRUCTION STRING 

CALL U I S $TEXT( VD_ I D , 7 , PW_KB_ INST( PLI NE_NO ) ,0.0,11.0) 
CALL UIS$TEXT(VD_ID, 7, '[RETURN] TO EXIT' ,0. 0 , 10. 6) 
C*****DECIDE IF WE WANT THE FIRST LINE (NAME) OR A VALUE 
IF (PLINE.NO .EQ. 1) THEN 
PW.FCHAR = .FALSE. 

ELSE 

PW_FCHAR = .TRUE. 

ENDIF 

CALL KEY_READ(PW_LINE,PW_FCHAR,0. 0, 12. 0, *300) 

DY = DEL_Y - . 35 

CALL UIS$ERASE(VD_ID,8. 2,DY,9. 9,DEL_Y) 

CALL UIS$TEXT(VD_ID,7,PW_LINE,8. 2,DEL_Y) 

RETURN 

300 RETURN 1 
END 



SUBROUTINE PW_SDI_LABEL( VD_ID) 

c 

C PERFORM DISPLAY VARIABLES IS A ROUTINE WHICH READS UP THE 
C DATA ONTO THE DISPLAY AREA ALONG WITH THE SCREEN LABELS. 

C 

C CALLED BY SUBROUTINES ' POWER' ,' POWER_KB_DATA_ IN' 

C 

IMPLICIT INTEGER (A-Z) 

INCLUDE ' TOP.POWER. FOR' 

INCLUDE ' SYS $ LIBRARY; UI SENTRY' 

INCLUDE 'SYS$LIBRARY: UISUSRDEF' 

REAL DEL_Y, DOT_Y 
CHARACTER P_LINE_LBL(20)*35 
CALL UIS$ERASE(VD_ID,4. ,10. 1,9. 9,18. 3) 



r’jV’jWf* 

* 

* 

* 






P_LINE_LBL(1) = 


'1. 


SHIP NAME 




P_LINE LBL(2) = 


'2. 


LENGTH (Lpp) 


....(FT).. 


P_LINE_LBL(3) = 


'3. 


DESIGN DRAFT 


....(FT).. 


P LINE LBL(4) = 


'4. 


BEAM 


....(FT).. 


P_LINE_LBL(5) = 


'5. 


DISPLACEMENT 


. . . (LTONS) 


P_LINE_LBL(6) = 


■6. 


BLOCK COEFICIENT. . 


. . . . (Cb). . 


P_LINE_LBL(7) = 


'7. 


SHIP SPEED 


(KTS). 


P LINE LBL(8) = 


'8. 


PROPELLER SPEED. . . 


(RPM). 


P LINE LBL(9) = 


'9. 


NUMBER OF SHAFTS. . 




P_LINE_LBL(10) = 
P_LINE_LBL(11) = 
P_LINE LBL(12) = 


1 * 
t T 
1 r 


****end of REQD INPUTS'^******* 


P.LINE LBL(13) = 




*** RESULTS*** 




P_LINE_LBL(14) = 


' * 


*****h0RSEP0WER. (HP) 


Vf 'sV’sV'jWf Vr ^Wf * 
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10 



20 



P_LINE_LBL(15) = 
P_LINE_LBL(16) = 
P_LINE_LBL(17) = 
P_LINE_LBL(18) = 
P_LINE_LBL(19) = 
P_LINE_LBL(20) = 
DO 10 I = 1,20,1 
DEL Y = 18. 6 



METHOD: 

ADMIRALTY 

S.&.D.W/.NO BB 

S.&. D. W/ BB 

BB=. WELL. DESIGNED. BULB. BOW 

•>V‘5V^Wr**Vf**Vc*Vc’*'>Wc**^V***:>V*******‘>V'>V 

.4*1 



CALL UIS$TEXT(VD_ID,7,P_LINE_LBL(I),4. ,DEL_Y) 
CONTINUE 

DO 20 I = 1,12,1 

DEL_Y =18.6 - .4*1 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(I),8. 2,DEL_Y) 
CONTINUE 
RETURN 
END 



SUBROUTINE PW CALC 1 






•k 



SUBROUTINE PW_CALC_1 * 

PW_CALC_1 CALLS PW_CALC_2 WHICH DOES THE POWER (HP) * 

CALCULATIONS. ON RETURN, PW_CALC_1 DISPLAYS THE RESULTS * 

OF THE CALCULATIONS AND THEN RETURNS TO THE ROOT PROGRAM. * 

Vf 

CALLED BY SUBROUTINES ' POWER' ,’ P0WER_KB_DATA_IN' ,' P0WER_GRAPH' * 
CALLS SUBROUTINE 'PW_CALC_2' * 

VfVc^VVf*V?:tVVfVrVc^V')V'5VVr*Vc**Vc*:V*:fV**^V***Vv:iV***V»Vc*’>VVc*VcTV:iVVcVcVfVc’^V'5V'>V*^V')V^^V'5VVvVc’'5V'5VVfVcVf** 



INCLUDE 'TOP. FOR’ 

INCLUDE ' TOP_POWER. FOR ' 

REAL DY 

CALL UIS$ERASE(VD_ID,8. 1 , 9. 95 , 9. 9 , 13. ) 

CALL PW CALC 2 

WRITE( PW_KbIuP(16), FMT=’ (G12. 5) ’ ) PW_ADME 
WRITE( PW_KB_UP(17), FMT=’ (G12. 5) ’ ) PW_SDNB 
WRITE( PW_KB_UP(18), FMT=' (G12. 5) ' ) PW_SDBB 
PW_KB_UP (10) = ' ************ ' 

pw_kbIup(ii) 

PW_KB_UP(12) 

PW_KB_UP(13) 

PW_KB_UP(14) 

PW_KB_UP(15) 

PW_KB_UP(19) 

PW_KB_UP(20) 

DO 10 I = 10,20,1 
DY = 18. 6 -.4*1 

CALL UIS$TEXT(VD_ID,7,PW_KB_UP(I),8. 2,DY) 
CONTINUE 
RETURN 
END 






= ' VcVrsWfVcvV?***** * 



SUBROUTINE PW_CALC_2 

C**Vf**********************'sV**************************Vc************* 
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c * 

C SUBROUTINE PW_CALC_2 JUST DOES THE CALCULATION OF THE POWER * 
C REQUIREMENTS FROM THE INPUT DATA. IT IS BASED UPON TWO * 
C ALGORITHMS FOR POWER PREDICTION: THE ADMIRALTY COEFFICIENT, * 
C AND THE METHOD OF SILVERLEAF AND DAWSON AS MODIFIED BY STIAN * 
C ERICHSEN AT UNIV. OF MICHIGAN. IT IS ALSO USED IN THE REPORT * 
C SECTION FOR GENERATION OF THE WRITTEN REPORT OF THE DESIGN. * 
C * 
C CALLED BY SUBROUTINE ’PW_CALC_l' * 
C CALLS SUBROUTINE 'PW_FAIL' * 



c 

INCLUDE 'TOP. FOR' 

INCLUDE ' TOP_POWER. FOR ' 

C*****FIRST CONVERT ENGLISH UtHTS TO METRIC FOR ADMIRALTY CALCULATION 
C*****knOTS->M/S ; FT->M ; LTONS->MTONS ; 

PW_NRPS = PW_NRPM/60. 0 
PW_LPPM = PW_LPP/3. 280 
PW_VKM = PW_VK*0. 514444 
PW_DISPM = PW_DISP/1. 016 
I_PW_NSHAFT=INT( PW_NSHAFT) 

C*****beGIN error checking loops for ADM CALCULATION 
C THIS SECTION CHECKS THE ADM FORMULA DATA FOR VALIDITY BEFORE 
C THE CALCULATION BEGINS. IF A PARAMETER IS OUT OF RANGE, THE 

C ERROR CODE IS SET AS 'PW_CAUSE' AND PASSED TO THE SUBROUTINE 

C 'PW_FAIL' WHICH DISPLAYS A HELP WINDOW SHOWING THE SPECIFIC 

C REASON THE DATA WAS REJECTED. THE ADM VALUES ARE SET TO 

C ZERO AND THE PROGRAM CONTINUES TO THE S&D CALCULATION. 

IF ( PW_LPPM .GT. 1162.196 ) THEN 
PW_CAUSE=1 

CALL PW_FAIL(PW_CAUSE) 

PW_ADME=0. 0 
GOTO 100 
END IF 

IF ( PW_NRPS*(PW_LPPM**0. 5) . GT. 136.364 ) THEN 
PW_CAUSE=2 

CALL PW_FAIL(PW_CAUSE) 

PW_ADME=0. 0 
GOTO 100 
ENDIF 

C*****CALCULATE POWER BY METHOD OF ADMIRALTY COEFFICIENTS 

PW_ADM=5. 0*((PW_DISPM**(2. 0))**(1. 0/3. 0) )*(PW_VKM**3. 0) 

& *(33. 0-0. 017*PW_LPPM)/( 15000. 0-110. 0*PW_NRPS*PW_LPPM**0. 5) 

C*****CONVERT METRIC KILOWATTS BACK TO HORSEPOWER 
PW_ADME = PW_ADM/0. 7456999 

C*****NOW WE HAVE ADMIRALTY POWER IN HORSEPOWER (HP) 

C 

100 CONTINUE 
C 

C*****NOW CALCULATE S&D POWER AND REQUIRED CORRECTIONS 
C 

C*****CALCULATE THE BOUNDARY SPEED AND SOME COEFFICIENTS 
PW_VB = (1. 70-1. 40*PW_CB)*(PW_LPP)**0. 5 
PW_SLR = PW_VK/(PW_LPP**0. 5) 

PW_BTR = PW_B/PW_T 
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PW_LBR = PW_LPP/PW_B 

c*****begin error checking loops for s&d calculation 

THIS SECTION CHECKS THE S&D FORMULA DATA FOR VALIDITY BEFORE 
THE CALCULATION BEGINS. IF A PARAMETER IS OUT OF RANGE, THE 
ERROR CODE IS SET AS 'PW_CAUSE' AND PASSED TO THE SUBROUTINE 
'PW_FAIL' WHICH DISPLAYS A HELP WINDOW SHOWING THE SPECIFIC 
REASON THE DATA WAS REJECTED. THE S&D VALUES ARE SET TO 
ZERO AND THE PROGRAM CONTINUES. 

IF (I_PW_NSHAFT . EQ. 1 ) THEN 

IF (PW_SLR .GT. 1.20 .OR. PW_SLR . LT. 0.40 ) THEN 
PW_CAUSE=3 

CALL PW_FAIL(PW_CAUSE) 

PW_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
END IF 

IF (PW_CB .GT. 0.86 .OR. PW_CB . LT. 0.50 ) THEN 

pw rAn^F=4 

CALL PW_FA I L( PW_C AUSE ) 

PW_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
END IF 

IF (PW_BTR .GT. 4.50 .OR. PW_BTR . LT. 2.0 ) THEN 
PW_CAUSE=5 

CALL PW_FAIL(PW_CAUSE) 

PW_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
ENDIF 

IF (PW_LBR .GT. 9.50 .OR. PW_LBR . LT. 3.33 ) THEN 
PW_CAUSE=6 

CALL PW_FAIL(PV^_CAUSE) 

PW_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
ENDIF 

ELSEIF (I_PW_NSHAFT . EQ. 2 ) THEN 

IF (PW_SLR .GT. 1.20 .OR. PW_SLR . LT. 0.40 ) THEN 
PW_CAUSE=3 

CALL PW_FA I L( PW_C AUSE ) 

P\C_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
ENDIF 

IF (PW_CB .GT. 0.80 .OR. PW_CB . LT. 0.54 ) THEN 
PW_CAUSE=4 

CALL PW_FAIL(PW_CAUSE) 

PW_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
ENDIF 

IF (PW_BTR .GT. 4.50 .OR. PW_BTR . LT. 2.0 ) THEN 
PW_CAUSE=5 

CALL PW_FAIL(PW_CAUSE) 

PW_SDNB=0. 0 
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PW_SDBB=0. 0 
GOTO 200 
ENDIF 

IF (PW_LBR .GT. 11.50 .OR. PW_LBR . LT. 3.80 ) THEN 
PW_CAUSE=6 

CALL PW_FAIL(PW_CAUSE) 

PW_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
ENDIF 
ELSE 

PW_CAUSE=7 

CALL PW_FAIL(PW_CAUSE) 

PW_SDNB=0. 0 
PW_SDBB=0. 0 
GOTO 200 
ENDIF 

C*****C0RRECT FOR SHIP LENGTH 

IF ( PW_LPP .GT. 1000.0 ) PW_CX=0. 850 

IF ( PW_LPP . LE. 1000. 0 ) PW_CX=0. 850+0. 00185* 

& ((1000. 0-PW_LPP)/100. 0)**2. 5 

C*****cORRECT FOR BEAM/DRAFT RATIO 

IF ( PW_B/PW_T . LT. 2. 40 ) C_KBT=0. 982 
IF ( PW_B/PW_T . GE. 2. 40 ) C_KBT=0. 960+ 

& 0. 00054*10. 0**((2. 0*PW_B)/(3. 0*PW_T)) 

C*****C0RRECT FOR SHIP SPEED NOT EQUAL TO BOUNDARY SPEED 

IF (PW_VK/PW_VB .LE. 1.0 ) C_KV=2. 75-7. 25*(PW_VK/PW_VB)+ 

6c 5 . 5 0* ( ( PW_VK/ PW_VB ) **2 . 0 ) 

IF (PW_VK/PW_VB .GT. 1.0 ) C_KV=21. 2-43. 2*(PW_VK/PW_VB)+ 

6c 2 3 . 0* ( ( PW_VK/ PW_VB ) ** 2.0) 

C*****CORRECT FOR HYDRODYNAMIC EFFICIENCY AT LPP=400. 0 FT 
IF ( I_PW_NSHAFT . EQ. 1 ) C_H400=2. 60-0. 2917*(PW_VK/ 

6c (PW_DISP**(1. 0/6. 0))) 

IF ( I_PW_NSHAFT . EQ. 2 ) C_H400=2. 38-0. 2917*(PW_VK/ 

6c (PW_DISP**(1. 0/6. 0))) 

C*****CORRECT FOR LPP NOT EQUAL TO 400. 0 FT 
IF ( PW.LPP .EQ. 400.0 ) C_KLBP=1. 0 

IF ( PW_LPP .NE. 400.0 ) C_KLBP=0. 9196+(PW_LPP*2. 31E-04)- 

6c (7. 5E-08*(PW_LPP**2. 0)) 

C*****CORRECT FOR OPEN WATER EFFICIENCY AT PROP SPEED=120 RPM 
IF ( I_PW_NSHAFT . EQ. 1 ) C_ETA_0120=0. 98-0. 55*PW_CB 

IF ( I_PW_NSHAFT . EQ. 2 ) C_ETA_0120=0. 90-0. 33*PW_CB 

C*****CORRECT FOR ACTUAL PROP SPEED RPM 

IF ( I_PW_NSHAFT . EQ. 1 ) C_DETA_0=0. 360-0. 0029*PW_NRPM 

IF ( I_PW_NSHAFT . EQ. 2 ) C_DETA_0=0. 135-0. 001 1*PW_NRPM 

C*****ADD THE TWQ ABOVE FOR TOTAL PROP SPEED CORRECTION 
PW_ETA_0=C_ETA_0 1 20+C_DETA_0 
C*****CALCULATE S6cD FOR NO BULBOUS BOW 

T0P_ANS=PW_CX*(PW_DISP**(2. 0/3. 0) )*(PW_VB**3. 0)*C_KBT^'^C_KV 
B0T_ANS=427 . 1*PW_ETA_0*C_H400*C_KLBP 
PW_SDNB=TOP_ANS/BOT_ANS 
C*****ADJUST S6cD FOR SHIP WITH BULBOUS BOW 
PW_SDBB=PW_SDNB*0. 9500 

C*****NOW WE HAVE S6cD POWER IN HORSEPOWER (HP) FOR BOTH BOW SHAPES 
C AND THE ADMIRALTY POWER. 

C 
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200 CONTINUE 
C 

RETURN 

END 

C 

C 

SUBROUTINE PW_RECORD 

Q^V:?c’Vf^V*^V^V******Vc**^V**:?f:Sr**:Wc:?r^V^V:iV*^V*Vf**yf****:?fVoV****yc:?f*:?f***9V***^V****^V** 



C * 
C SUBROUTINE 'PW.RECORD* SAVES THE KBD INPUT DATA TO A USER * 
C SPECIFIED FILE. FILE TYPE IS '.DAT' AUTOMATICALLY * 
C * 
C CALLED BY SUBROUTINE 'POWER'. * 
C CALLS SUBROUTINE 'SHOW_SAVE' * 



C*‘5vyf*yfycyryfyr^v****ycyryfyryf*yrVrVp*yf*yfyryfyr^vyfyc^vyryr*yf*yryr*ycyfyryryovyfycycyr^vyr*vcyryryfyryrycyryryfyc 



c 



c 



INCLUDE 'TOP. FOR' 

INCLUDE 'TOP_POWER. FOR' 

STA = UIS$SET_P0INTER_P0SITI0N(VD_ID,WD_P0WER,0. ,10. ) 
CALL SHOW_SAVE 

0 PEN( 2 1 , F I LE=PW_I NF I LE , STATUS= ' UNKNOWN ' , b 1 ank= 'null') 

WRITE (21,200) PW_INFILE 

WRITE (21,210) PW_LPP 

WRITE (21,210) PW_T 

WRITE (21,210) PW_B 

WRITE (21,210) PW_DISP 

WRITE (21,210) PW_CB 

WRITE (21,210) PW_VK 

WRITE (21,210) PW_NRPM 

WRITE (21,210) PW_NSHAFT 

WRITE (21,210) PW_ADME 

WRITE (21,210) PW_SDNB 

WRITE (21,210) PW_SDBB 

ENDFILE (21) 

CL0SE(21) 



c*****and report out when done 

C 

PW_FNAME = PW_INFILE( : INDEX(PW_INFILE, ' '))//'. DAT' 
CALL UIS$ERASE(VD_ID,-. 4,10. 1,3. 4,10. 9) 

CALL UIS$TEXT(VD_ID,7, 'FILE SAVED AS',0. ,11. ) 

CALL UIS$TEXT(VD_ID,7,PW_FNAME,0. ,10.6) 

200 FORMAT (BN,2X,A) 

210 FORMAT (G12.5) 

RETURN 

END 

C 

c 



SUBROUTINE PW_READ FILE 



C 

C 

C 

C 

C 

C 

C 



.t. y > ^ .1. ^ y. ^ 

* 



READ FILE IS A SUBROUTINE WHICH READS THE DATA INPUT FILE * 
SPECIFIED BY KEYBOARD ENTRY * 



CALLED BY SUBROUTINE 'POWER' * 

CALLS SUBROUTINES ' PW_SDI_LABEL' , ' PW_CALC_1 ' * 
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(]***********w*iVVc*********i*r*'jV'!V**iWf***iWf**Vcyf**AiSf**»V*Vf**?Wf***Vc***5Wf*** 

C 

INCLUDE 'TOP. FOR’ 

INCLUDE ' TOP_POWER. FOR ' 

STA=UIS$SET_POINTER_POSITION(VD_ID,WD_POWER,0. ,10. 2) 
KB_ID=UIS$CREATE_KB( ' SYS$WORKSTATION' ) 

CALL U I S $ ENAB LE_KB ( KB_ I D , WD_P0WER ) 

CALL UIS$ERASE(VD_ID,-. 4,10. 1,2. 3,10. 9) 

CALL UIS$TEXT(VD_ID,7, 'ENTER THE FILE NAME’ ,0. ,11. ) 

CALL UIS$TEXT(VD_ID,7,'AND FILE EXTENSION' ,0. , 10. 6) 

CALL KEY_READ(PW_LINE, 'TRUE' ,0. ,12. ,*150) 

OPEN( 2 1 ,FILE=PW_LINE , STATUS= ' UNKNOWN ' ) 

READ (21,200) PW_INFILE 
READ (21,210) PW_LPP 
READ (21,210) PW_T 
READ (21,210) PW_B 
READ (21,210) PW_DISP 
READ (21,210) PW_CB 
READ (21,210) PW_VK 
READ (21,210) PW_NRPM 
READ (21,210) PW_NSHAFT 
REWIND 21 

C**** NOW READ IT IN AS CHARACTER DATA TO PRINT TO SCREEN 
READ (21,200) PW_KB_UP(1) 

DO 10 I = 2,9,1 

READ (21,220) PW_KB_UP(I) 

10 CONTINUE 

C 

200 FORMAT (2X,A) 

210 FORMAT (G12.5) 

220 FORMAT (2X,A) 

150 CLOSE(21) 

CALL PW_SDI_LABEL(VD_ID) 

CALL UIS$DISABLE_KB(KB_ID) 

CALL PW_CALC_1 

CALL UIS$ERASE(VD_ID,-. 4,10. 1,3. 4,10. 9) 

CALL UIS$TEXT(VD_ID, 7, 'SELECT AN OPTION' ,0. ,11.) 

CALL UIS$TEXT(VD_ID,7, ' WITH THE MOUSE ',0. ,10.6) 

RETURN 

END 



SUBROUTINE SAVE_GRAPH 



ycVfVrVc VcVr VcVc Vc V-T Vc Vc VcVc Vc VirVtVr Vc VrVc Vc Vc V? Vf VirVt VrVc Vc Vc VrVr Vc* 






THIS SUBROUTINE POPS UP A DIALOG BOX TO INSTRUCT THE USER THAT * 
THE POWER PREDICTION CURVE GRAPH PLOT HAS BEEN WRITTEN TO DISK * 



UNDER THE NAME ' POWER_PLOT. UIS ' AND INSTRUCTS THEM ON HOW TO * 
RETRIEVE THE FILE IN HARDCOPY FORM. IT IS CALLED BY THE SUB- * 
ROUTINE 'power' SUBMENU SELECTION 'GRAPH POWER’ AND IS * 

DISPLAYED AFTER THE GRAPH HAS BEEN CREATED, VIEWED, AND SENT * 
TO DISK. * 

* 

CALLED BY SUBROUTINE ’ POWER_GRAPH ' * 
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c 

IMPLICIT INTEGER(A-Z) 

INCLUDE ’ SYS$LIBRARY: UISENTRY' 

INCLUDE ' SYS$LIBRARY: UISUSRDEF' 

I NCLUDE ' TOP_POWER . FOR ' 

CHARACTER*34 PLOTA/'*** SAVING DEFINED GRAPH FILE ***'/ 
CHARACTER*34 PLOTB/' GRAPH IS BEING WRITTEN TO DISK '/ 
CHARACTER*34 PLOTC/' USING THE FILENAME: 7 

CHARACTER*34 PLOTD/ ' "STD00001.DAT" 7 

CHARACTER*34 PLOTE/' YOU CAN OBTAIN HARDCOPY PLOTS 7 
CHARACTER*34 PLOTF/' BY THE DCL COMMAND: 7 

CHARACTER*34 PLOTG/ ' "PRINT/QUEUE=LASER STDOOOOl. DAT" 7 
CHARACT£R*34 PLOTH/' AFTER TERMINATION OF TOOLBOX. 7 

CHARACTER*34 PLOTI/ ' PLEASE WAIT 7 

CHARACTER*34 PLOTJ/' PROGRAM WILL RESUME 7 

REAL Y_POSN 

c**,v**create the display for the dialog box 

VD_PLOT=UIS$CREATE_DISPLAY(-5. 0,-2.5,10. 0,5. 0,26. 0,14. 0) 
C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 

WD_PLOT=UIS$CREATE_WINDOW(VD_PLOT, 'SYS$WORKSTATION' , 'HELP WINDOW' ) 
C*****C0PY ATTRIBUTE BLOCK 'O' AS BLOCK '23' AND CHANGE THE FONT SIZE 

CALL UIS$SET_FONT(VD_PLOT,0,23, 'DTABER0R03WK00GG0001UZZZZ02A000' ) 
C*****COPY ATTRIBUTE BLOCK '23' AS BLOCK '24' AND CHANGE TO BOLD FONT 

CALL UIS$SET_FONT(VD_PLOT,23,24,'DTABER0R03WK00PG0001UZZZZ02A000' ) 

c*****signal the user to get their attention 

CALL UIS$SOUND_BELL( 'SYS $ WORKSTATION' ,4) 

C*****WRITE THE TEXT INTO THE WINDOW AND SPAWN PROCESS TO 'WAIT' 15 SEC 
Y_P0SN=4. 0 

CALL UIS$TEXT(VD_PL0T,24,PL0TA,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,23,PL0TB,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,23,PL0TC,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,24,PL0TD,-2.5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,23,PL0TE,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,23,PL0TF,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,24,PL0TG,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,23,PL0TH,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,23,PL0TI,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_PL0T,23,PL0TJ,-2. 5,Y_P0SN) 

CALL LIBSSPAWNC 'WAIT 00:00: 15') 

CALL UIS$S0UND_BELL('SYS$W0RKSTATI0N' ,4) 

CALL UIS$DELETE_DISPLAY(VD_PLOT) 

RETURN 

END 
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SUBROUTINE PW_FAIL( PW.CAUSE) 



c * 

C THIS SUBROUTINE POPS UP A DIALOG BOX TO INSTRUCT THE USER THAT * 
C THE POWER PREDICTION ALGORITHM WILL FAIL TO ACURATELY PREDICT * 

C PROPULSIVE POWER REQUIREMENTS DUE TO THE INPUTED DATA BEING * 

C OUTSIDE OF THE PARAMETERS ALLOWED BY THE SILVERLEAF & DAWSON * 

C OR ADMIRALTY FORMULAS. THE SPECIFIC CAUSE IS PASSED AS THE * 

C PARAMETER ’PW_CAUSE' TO THE SUBROUTINE FROM THE ERROR CHECKING * 

C LOOPS IN SUBROUTINE PW_CALC_2. THIS PARAMETER ENABLES THIS * 

C SUBROUTINE TO BE WRITTEN ONCE FOR ALL CASES AND STILL BE * 

C SPECIFIC FOR EACH INDIVIDUAL CASE. * 

C * 

C THE ALGORITHMS USED ARE DETAILED IN SUBROUTINE PW_CALC_2 * 

c * 

C CALLED BY SUBROUTINE ’ PW CALC 2' * 






IMPLICIT INTEGER(A-Z) 

INCLUDE 'SYS$LIBRARY: UISENTRY’ 

INCLUDE ’SYS$LIBRARY: UISUSRDEF' 

INCLUDE ’T0P_P0WER. FOR' 

CHARACTER*34 FAILA/'*** POWER PREDICTION FAILURE ***’/ 
CHARACTER*34 FAILB/' TOOL BOX MAY FAIL TO ACCURATELY '/ 
CHARACTER-'-34 FAILC/' PREDICT POWER REQUIREMENTS FOR '/ 
CHARACTER*34 FAILD/' ADMIRALTY OR S & D METHOD DUE TO '/ 

CHARACTER*34 FAILE/ ' PLEASE WAIT '/ 

CHARACTER*34 FAILF/' PROGRAM WILL RESUME'/ 

REAL Y_POSN 
CHARACTER*36 FAIL(7) 

c*****initialize the failure messages array 

FAIL(l) = ' SHIP LENGTH OUT OF RANGE ' 

FAIL(2) = ' SHIP LENGTH OR NRPM OUT OF RANGE ' 

FAIL(3) = ' SPEED/LENGTH RATIO OUT OF RANGE ' 

FAIL(4) = ' SHIP BLOCK COEFFICIENT OUT OF RANGE' 

FAIL(5) = ' BEAM/DRAFT RATIO OUT OF RANGE ' 

FAIL(6) = ' LENGTH/BEAM RATIO OUT OF RANGE ' 

FAIL(7) = ' NUMBER OF SHAFTS INCORRECT (1 OR 2)' 

C*****creATE THE DISPLAY FOR THE DIALOG BOX 

VD_FAIL=UIS$CREATE_DISPLAY(-5. 0,-2. 0,10. 0,4. 0,26. 0, 10. 0) 
C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 

WD_FAIL=UIS$CREATE_WINDOW(VD_FAIL, 'SYS$WORKSTATION' , 'HELP WINDOW' ) 
C*****C0PY ATTRIBUTE BLOCK 'O' AS BLOCK '25' AND CHANGE THE FO.NT SIZE 

CALL UIS$SET_FONT( VD_FAIL , 0,25,' DTABER0R03WK00GG000 1UZZZZ02A000 ' ) 
C*****C0PY ATTRIBUTE BLOCK '25' AS BLOCK '26' AND CHANGE TO BOLD FONT 

CALL UIS$SET_FONT( VD_FAIL, 25 , 26 , ' DTABER0R03WK00PG0001UZZZZ02A000 ’ ) 
C*****SIGNAL THE USER TO GET THEIR ATTENTION 

CALL UIS$SOUND_BELL( 'SYS$WORKSTATION' ,4) 

C*****writE the text into THE WINDOW AND SPAWN PROCESS TO 'WAIT' 15 SEC 
Y_POSN=2. 8 

CALL UIS$TEXT(VD_FAIL,26,FAILA,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_FAIL,25,FAILB,-2. 5,Y_POSN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_FAIL,25,FAILC,-2. 5,Y_POSN) 

Y POSN=Y POSN-.6 
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CALL UIS$TEXT(VD_FAIL,25,FAILD,-2. 5,Y_P0SN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_FAIL,26,FAIL(PW_CAUSE) ,-2. 5,Y_P0SN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_FAIL,25,FAILE,-2. 5,Y_P0SN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT( VD_FAIL , 25 , FAILF , -2. 5 , Y_POSN) 

CALL LIB$SPAWN( 'WAIT 00:00:15') 

CALL UIS$SOUND_BELL( 'SYS$W0RKSTATI0N' ,4) 

CALL UIS$DELETE_DISPLAY(VD_FAIL) 

RETURN 

END 



SUBROUTINE POWER_GRAPH(TYPE) 






* 



THIS SUBROUTINE CREATES A CA-DISSPLA METAFILE GRAPH OF THE * 
PROPULSIVE POWER CURVE OF THE SHIP DESIGN FROM ZERO (0) * 
KNOTS TO THE OPERATING SPEED 'PW_VK' GIVEN AT THE STAGE * 
OF DESIGN WHEN THE SUBROUTINE WAS CALLED. IT IS CALLED BY * 
A MOUSE CLICK ON THE MAIN POWER MENU 'GRAPH POWER' SO IT * 
IS IMPORTANT THAT THE DESIGN PARAMETERS BE FINALIZED AT * 
THAT POINT SINCE ALL THE REQUIRED PARAMETERS WILL BE PASSED * 
IN COMMON BLOCKS. IT ALSO CALLS THE SUBROUTINE ' SAVE_GRAPH'* 
THAT ANNOUNCES THE METAFILE NAME TO THE USER BEFORE IT * 
WRITES THE FILE TO DISK. * 

* 



CALLED BY SUBROUTINE 'POWER' 
CALLS SUBROUTINES 'SAVE GRAPH' 



* 



AND MANY CA-DISSPLA ROUTINES * 



INCLUDE 'TOP. FOR' 

INCLUDE 'TOP.POWER. FOR' 

CHARACTER*23 DATETIME , PWG_NAME 

REAL LM,LF,BF,TF,N,DELTAM,DELTAT,VM,VK(0: 15 ) , PHP_AD , VKP , 

& CB,C_VB,C_X,KBT,H400,KV,ETA0120,DETA0,C_ETA0,VMTK(0: 15), 

& GPW_SD(0: 15),GPW_SDB(0: 15) ,KLBP,C_VL,C_BT,C_LB,VMT(0: 15), 

& TOP_SD,TOP_SD1,BOT_SD,NRPM,GPW_ADM(0: 15) ,GPW_AD.MP,VMG, 

& MIN.POWERl ,MIN_P0WER2 ,NAX_PWR 

INTEGER I_NSHAFT , INDEXl , INDEX2 , I , J , K , L , PLOT_TYPE , TYPE 
INTEGER*4 STATUS 
PLOT_TYPE=TYPE 

C*****ZERO OUT ALL MATRICES AND VARIABLES 
DO 5 1=0,15,1 
VK(I)=0. 0 
VMTK(I)=0. 0 
GPW_SD(I)=0. 0 
GPW_SDB(I)=0. 0 
VMT(I)=0. 0 
GPW_ADM(I)=0. 0 
5 CONTINUE 

NIN_P0WER1=0. 0 
MIN_POWER2=0. 0 
INDEX1=0 
INDEX2=0 
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MAX_P0WER=0. 0 

c*****recalculate common values for accuracy 

CALL PW_CALC_1 

C*****teLL world to sit tight and wait 

CALL UIS$ERASE(VD_ID, -0.4,10. 1,3.4,10. 9) 

CALL UIS$TEXT(VD_ID, 7, 'CALCULATING PLOT. . . ' ,0. 0, 11. 0) 

CALL UIS$TEXT(VD_ID, 7, 'PLEASE RELAX & WAIT' , 0. 0 , 10. 6) 

c*****transfer values via common blocks 

PWG_NAME=PW_INFILE 
I_NSHAFT=INT( PW_NSHAFT) 

LF=PW_LPP 

BF=PW_B 

TF=PW_T 

DELTAT=PW_DISP 

VKP=PW_VK 

CB=PW_CB 

NRPM=PW_NRPM 

0*v?***set up ca-disspla environment******************************* 

C 'LN03I' SENDS OUTPUT TO A LASER PRINTER CAPABLE FILE; 'PAGE' * 
C SETS PAGE SIZE AS STANDARD, ' BLOWUP' INCREASES PICTURE SIZE; * 
C 'AREA2D' SETS 2D PLOT SIZE; 'TRIPLX' SELECTS THE FONT FILE; * 
C 'X/Y name's SET AXIS NAMES; 'HEADIN' SETS THE GRAPH TITLE; * 
C 'GRAF' SETS AXIS SIZES(ORIGIN,MAX X,Y ,UNITS , SCALE , ETC); * 

C 'THKFRM' DRAWS A THICK LINE AROUND THE SUBPLOT AREA; 'FRAME' * 
C FRAMES THE ENTIRE PAGE; 'GRID' SELECTS GRAPH GRIDDING: * 

C 'RASPLN' SELECTS TYPE OF CURVE SMOOTHING/ INTERPOLATION TO USE;* 
C 'HEIGHT' SETS CHARACTER HEIGHTS AS WE WANT THEM TO BE; * 

Q ^ y .j. .. y Vf Vf V;’ Vc* V' "ti y V fV "if “ti it it it it it it it it it it it it it it 

C*****deTERMINE if desire plot to disk (1) OR TO SCREEN (2) 

IF ( PLOT_TYPE .EQ. 1 ) THEN 
CALL PGPX 

ELSEIF ( PLOT_TYPE . EQ. 2 ) THEN 
CALL LN03I 

ELSE 

CONTINUE 

ENDIF 

CALL PAGE(8. 5,11. 0) 

CALL BLOWUPCl. 15) 

CALL AREA2D(6. 0,8. 0) 

CALL TRIPLX 
CALL HEIGHTCO. 200) 

CALL XNAMECSPEED (KTS)$',100) 

CALL YNAME( 'horsepower (HP)$',100) 

CALL HEIGHTCO. 375) 

CALL HEADINC 'SPEED-POWER RELATIONSHIP$' , -24, -1. , 1) 

CALL HEIGHTCO. 200) 

MAX_PWR=MAXC PW_ADME , PW_SDBB , PW_SDNB ) 

CALL GRAFCO. 0,2. 0,PW_VK,0. 0,5000. 0,CMAX_PWR+2000. 0)) 

CALL THKFRMCO. 04) 

CALL FRAME 
CALL GRIDC-2,-1) 

CALL RASPLNC2.0) 

C*****SET UP ADMIRALTY POWER ARRAY IN 15 INCREMENTS 
LM=LF/3. 28 
DELTAM=DELTAT/1. 016 
VM=PW_VK*0. 514444/15. 0 
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N=NRPM/60. 0 
DO 10 1=0,15,1 
VMT(I)=I*VM 
VMTK(I)=PW_VK*I/15. 0 

GPW_ADMP=5. 0*((DELTAM**2. 0)**( 1. 0/3. 0) )*( VMT(I)**3. 0)* 

& (33. 0-(0. 017*LM))/( 15000. 0-110. 0*N*LM**0. 5) 

GPW_ADM(I)=GPW_ADMP/0. 7456999 
10 CONTINUE 

C*****ALSO, ACCOUNT FOR CASE THAT ADMIRALTY METHOD HAS FAILED 
IF ( PW_ADME .LT. 1.0) THEN 
DO 12 J=l,15,l 

GPW_ADM( J)=0. 0 
12 CONTINUE 

ENDIF 

c*****graph admiralty power curve 

CALL MARKER(16) 

CALL DOT 

CALL CURVE (VMTK,GPW_ADM, 16,1) 

CALL HEIGHTCO. 070) 

CALL LINES( 'ADM. POWER’, ID, 1) 

C*****fILL in S&D array in 15 INCREMENTS TO PLOT CURVE 
DO 20 1=0,15,1 

VK(I)=PW_VK*I/15. 0 

VMTK( I )=PW_VK*I / 15.0 

C_VB=( 1. 70-1. 4*CB)*LF**0. 5 

IF ( LF .GT. 1000.0 ) C_X=0. 85 

IF ( LF .LE. 1000.0 ) C_X=0. 85+0. 00185* 

& ((1000. 0-LF)/100. 0)**2. 5 

IF ( BF/TF .LT. 2.40 ) KBT=0. 982 
IF ( BF/TF .GE. 2.40 ) KBT=0. 96+ 

& 0. 00054*10. 0**((2. 0*BF)/(3. 0*TF)) 

IF ( VK(I)/C_VB .LE. 1.0 ) KV=2. 75-7. 25*(VK( I)/C_VB)+ 

& 5. 50*((VK(I)/C_VB)**2. 0) 

IF ( VK(I)/C_VB .GT. 1.0 ) KV=21. 2-43. 2*( VK( I)/C_VB)+ 

& 23. 0*((VK(I)/C_VB)**2. 0) 

IF ( I_NSHAFT. EQ. 1) H400=2. 60-0. 2917*VK( I)/(DELTAT**( 1. 0/6. 0)) 
IF (I_NSHAFT. EQ. 2) H400=2. 38-0. 2917*VK( I )/(DELTAT**( 1. 0/6. 0) ) 
IF ( LF .EQ. 400.0) KLBP=1. 0 

IF ( LF .NE. 400.0) KLBP=0. 9196+2. 31E-04*LF-7. 50E-08*LF**2. 0 
IF ( I.NSHAFT .EQ. 1 ) ETA0120=0. 98-0. 55*CB 

IF ( I_NSHAFT .EQ. 2 ) ETA0120=0. 90-0. 33*CB 

IF ( I.NSHAFT .EQ. 1 ) DETA0=0. 360-0. 0029*NRPM 

IF ( I_NSHAFT .EQ. 2 ) DETA0=0. 135-0. 0011*NPRM 

C_ETAO=ETAO 1 2 O+DETAO 
C 

T0P_SD=C_X*(DELTAT^'*(2. 0/3. 0) )*(C_VB**3. 0)*KBT^--KV 
BOT_SD=427. 10*C_ETA0*H400*KLBP 
C 

GPW_SD(I) =TOP_SD/BOT_SD 
GPW_SDB(I) = 0. 9500*(T0P_SD/B0T_SD) 

C 

20 CONTINUE 

C*****SEARCH THE POWER ARRAYS TO FIND THE MINIMUMS AND THEIR INDICES. 
C THIS HAS TO BE DONE DUE TO THE INVALIDITY OF THE S&D FORMULAS 
C AT SPEEDS « BOUNDARY SPEED. THIS ERROR CAUSES THE UN- 

C CORRECTED GRAPHS TO START AT HIGH HORSEPOWER AT LOW SPEED 
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(AN UNREALISTIC ASSUMPTION), DROP TO A MINIMUM (THAT WE FIND), 
AND THEN GROW ASYMPTOTICALLY TO THE TRUE PREDICTION. THUS WE 
WILL SIMPLY FIND THE MINIMUM AND SET ALL LOW SPEED VALUES TO THE 
LEFT OF THIS POINT EQUAL TO THAT MINIMUM FOR EACH CURVE. 
C*****ALSO, ACCOUNT FOR CASE THAT S&D METHOD HAS FAILED 

IF ( PW_SDNB . LT. 1. 0 . AND. PW_SDBB . LT. 1. 0 ) THEN 
DO 22 J=l,15,l 

GPW_SD(J)=0. 0 
GPW_SDB(J)=0. 0 
22 CONTINUE 

ENDIF 

MIN_POWERl=GPW_SD( 0) 

MIN_P0WER2=GPW_SDB(0) 

DO 25 J=l,15,l 

IF ( GPW_SD(J) .LT. MIN_P0WER1 ) THEN 
MIN_P0WER1=GPW_SD(J) 

INDEX1=J 

ELSE 

CONTINUE 

ENDIF 

IF ( GPW_SDB(J) .LT. MIN_POWER2 ) THEN 
MIN_POWER2=GPW_SDB( J ) 

INDEX2=J 

ELSE 

CONTINUE 

ENDIF 

25 CONTINUE 

C*****sET ANY POWER TO LEFT OF S&D MINIMUM EQUAL TO MINIMUM 
DO 26 K=INDEX1,0,-1 
GPW_SD(K)=MIN_P0WER1 

26 CONTINUE 

DO 27 L=INDEX2,0,-1 
GPW_SDB(L)=MIN_POWER2 

27 CONTINUE 

C*****GRAPH BOTH S&D CURVES 
CALL MARKER(15) 

CALL DASH 

CALL CURVE (VMTK,GPW_SD, 16,1) 

CALL HEIGHT(0. 070) 

CALL LINES('S&D POWER W/0 BB',ID,2) 

CALL MARKER(17) 

CALL CHNDOT 

CALL CURVE (VMTK,GPW_SDB, 16,1) 

CALL HEIGHT(0. 070) 

CALL LINES('S&D POWER W/ BB',ID,3) 

CALL HEIGHT(0. 150) 

CALL LEGNAM( ’***POWER CURVE S***$ ', 18) 

CALL LEGEND(ID,3,1. 75,6. 8) 

C*****STAMP FILENAME ON GRAPH SO ORIGIN IS KNOWN 
CALL HEIGHT(0. 100) 

CALL MESSAG(PWG_NAME,23,4. 5,-0. 75) 

C*****GET SYSTEM DATE&TIME AND STAMP ON GRAPH 
STATUS=LI B $DATE_TIME ( DATETIME ) 

CALL HEIGHT(0. 100) 

CALL MESSAG(DATETIME,23,4. 5,-0. 90) 

C***v.'*doNE 
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CALL ENDPL(O) 

CALL DONEPL 

C**,v**deTERMINE if plot went to disk (1) OR TO SCREEN (2) 

IF ( PLOT_TYPE .EQ. 1 ) THEN 
CONTINUE 

ELSEIF ( PLOT.TYPE . EQ. 2 ) THEN 
CALL SAVE_GRAPH 

ELSE 

CONTINUE 
END IF 

C*****CLEAR THE INSTRUCTION BOX SO THE USERS KNOWS ITS DONE 
CALL UIS$ERASE(VD_ID,-0. 4,10. 1,3.4,10.9) 

50 RETURN 
END 
C 
C 
C 

SUBROUTINE G_DISK 



c 

c 

c 

c 

c 

c 

G 



* 

THIS SUBROUTINE CALLS THE CA-DISSPLA ROUTINE TO PLOT TO * 
DISK FILE. * 

* 



CALLED BY SUBROUTINE MAIN MENU 'POWER' * 

CALLS SUBROUTINE ' POWER_GRAPH' * 












INCLUDE 'TOP. FOR' 

TYPE = 1 

CALL POWER_GRAPH(TYPE) 

RETURN 

END 



C 

C 

C 



SUBROUTINE G_SCREEN 



C 

C 

C 

C 

C 

C 

C 

C 





















* 



THIS SUBROUTINE CALLS THE CA-DISSPLA ROUTINE TO PLOT TO 
SCREEN. * 

* 



CALLED BY SUBROUTINE MAIN MENU 'POWER' * 

CALLS SUBROUTINE 'POWER GRAPH' * 









INCLUDE 'TOP. FOR' 

TYPE = 2 

CALL POWER_GRAPH(TYPE) 

RETURN 

END 



C 

C 

C 

C 
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B. SUBROUTINE POWER VARIABLE DECLARATIONS FILE 

This file is the INCLUDE File TOP_POWER.FOR that is common to ever>' sub- 
routine in the Power Prediction module and Power Prediction Report subroutine. This 
commonality is made to ensure that all data is available to every subroutine without 
needing to pass every parameter as an argument. 



*****the declarations for the propulsive power section in one file 

c 

CHARACTER PW_LINE*12 , PW_INFILE*12 
CHARACTER PW_KB_UP( 20)*12 , PW_FNAME*12 

INTEGER PW_TIMES , PW.CAUSE , I_PW_NSHAFT, PLINE_NO 

REAL PW_LPP,PW_T,PW_B,PW_DISP,PW CB,PW_VK 
REAL PW.NRPM , PW_SDNB , PW_SDBB , PW_NSHAFT 
REAL PW_VB , PW_CX , C_KBT , C_KV , C_H400 , C_KLBP 
REAL C_ETA_0 120 , C_DETA_0 , PW_ETA_0 , PW.TEMP 
REAL PW_ADME , PW_ADM , PW_S LR , PW_BTR , PW_LBR 
REAL PW_NRPS , PW_LPPM , PW_VKM , PW_D I SPM 
REAL T0P_ANS,B0T_ANS 

COMMON /PW_NAMES/ PW_LPP,PW_T,PW_B,PW_DISP,PW_CB,PW_VK 
COMMON /PW_VARS/ PW_NRPM,PW_NSHAFT,PW_TIMES,PW_KB_UP 
COMMON /PW_C0EFF/ PW_VB,PW_CX,C_KBT,C_KV,C_H400,C_KLBP 
COMMON /C0EFFS2/ C_ETA_0120,C_DETA_0,PW_ETA_0,PW_TEMP 
COMMON / PW_RESULTS / PW_SDNB , PW_SDBB , PW_ADME , PW_ADM 
COMMON /PW.METRIC/ PW_NRPS ,PW_LPPM,PW_VKM,PW_DISPM 
COMMON /PW_FILES/ PW_FNAME,PW_INFILE 
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APPENDIX C POWER PREDICTION TEST ROUTINE 



A. SOURCE CODE 

This code was used to verify and correct the Silverleaf and Dawson prediction for- 
mulas. It also makes an excellent 'stand- alone' program for general power calculations 
outside of the TOOL BOX environment. 



PROGRAM TRIAL_POWER 



C * 
C THIS PROGRAM TESTS THE SILVERLEAF & DAWSON PROPULSIVE POWER * 
C PREDICTION FORMULA FOR SOME REALISTIC DESIGNS WITH KNOWN * 
C PROPULSION REQUIREMENTS. THE ADMIRALTY PREDICTION IS ALSO * 
C CALCULATED FOR COMPARISON AND THE DATA IS PRESENTED IN * 
C TABLE 1 OF CHAPTER 3. * 
C * 
C REQUIRED INPUTS ARE AS FOLLOWS: * 
C SHIP LENGTH(LPP) IN FEET * 
C SHIP BEAM(B) IN FEET * 
C SHIP DRAFT(T) IN FEET * 
C SHIP DISPLACEMENTC DELTA) IN TONS STD SEAWATER * 
C NUMBER OF PROPULSION SHAFTS ( 1 OR 2 ) * 
C PROPELLER BLADE SPEED(N) IN REV/ SEC 

C SHIP SPEED(V) IN KNOTS * 



C 

C 

C 

C 

C 

C 

C 

C 



THIS PROGRAM IS WRITTEN IN VAX FORTRAN AND WILL ALSO COMPILE * 
UNDER LAHEY FORTRAN LP77 V3. 0 AND MICROSOFT FORTRAN 4.01 FOR * 
MS-DOS BASED SYSTEMS WITH A CHANGE TO THE OUTPUT FILENAME TO * 
ADHERE TO MS-DOS CONVENTIONS. * 



Vf*-.V*-!V***VrV.-***Vf***')Wf**Vr'!VVr****;V*'sWf******' 



itirlrifirlrifitifltirisiHfiHdrirlrir 



•5V*Vf*>W.-* 



REAL LM , LF , BF , TF , N , DELTAM , DELTAT , VM , VK , PHP_AD , PKW , NRPM , 

& C B , C_VB , C_X , C_KBT , C_KV , C_H4 00 , ETAO 120, DETAO , C_ETA0 , 

& PHP_SD , C_KLBP , PHP ISD , C_SL , C_BT , C_LB , 

& TOP_SD , BOT_SD 

INTEGER NSHAFT.IANSWER 
C*****OPEN AN OUTPUT FILE 

OPEN ( UNIT=8 , FILE=' TRIAL_POWER. TXT' , STATUS= ' UNKNOWN ' ) 

1 CONTINUE 

C*****INITIALIZE ( RE-INITIALIZE ) 'CORRECTING* VARIABLES 
VK=0. 0 
C_KV=0. 0 
C 

C*****QUERY FOR REQUIRED DATA 
C 

PRINT *, 'ENTER SHIP LENGTH IN FEET : ’ 

READ *,LF 

PRINT *, 'ENTER SHIP BEAM IN FEET : ' 

READ *,BF 

PRINT *, 'ENTER SHIP DRAFT IN FEET : ' 
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READ *,TF 

PRINT *, 'ENTER DISPLACEMENT IN TONS : ' 

READ *,dELTAT 

PRINT *, 'ENTER NUMBER OF SHAFTS AS AN INTEGER : ' 

READ ’’SNSHAFT 

PRINT */eNTER propeller SPEED IN REV/SEC ; ' 

READ *,N 

PRINT *, 'ENTER SHIP SPEED IN KNOTS : ' 

READ *,VK 
LM=LF/3. 28 
DELTAM=DELTAT/1. 016 
VM=VK*0. 514444 
C 

C*****erR0R check LENGTH(M),N*SQRT[LENGTH(M)] for ADMIRALTY METHOD 
C 

IF ( LM .GT. 1162.196 ) THEN 

PRINT * ******* length out OF RANGE( ADM) *****************' 
GOTO 13 ’ 

ELSE 

CONTINUE 

ENDIF 

IF ( N*(LM**0. 5) .GT. 136.364 ) THEN 

PRINT *,'******PRODUCT N*L(l/2) OUT OF RANGE (ADM)********' 
GOTO 13 

ELSE 

CONTINUE 

ENDIF 



C 

C*****CALCULATE POWER BY ADMIRALTY METHOD 
C 

PRINT *,' STARTING ADMIRALTY CALCULATION' 



C 



PK\\’=5. 0*((DELTAM**2. )**(!. 0/3. 0))*(VM**3. 0)*(33. 0-(0. 017*LM) 
& )/( 15000. 0-110. 0*N*LM**0. 5) 

c*****translate to horsepower from kilowatts 



PHP_AD=PK'W/0. 7456999 
C*****OUTPUT TO FILE AND SCREEN 
WRITE (8,7) 

WRITE (8,28) 

WRITE (8,8) 

WRITE (8,9) PHP_AD 
WRITE (8,7) 

WRITE (8,91) 

WRITE (8,92) 

WRITE (8,10) LM,DELTAM 
WRITE (8,11) 

WRITE (8,12) VM,N 

PRINT *,' POWER IN HP BY ADMIRALTY METHOD = 

C 



13 WRITE (8,7) 
C 



, PHP_AD 



C*****ERROR CHECK BLOCK COEFFICIENT, SPEED LENGTH RATIO, BEAM/DRAFT RATIO 
C FOR SILVERLEAF AND DAWSON METHOD 
C 



NRPM=N*60. 0 

CB=35 . 0*DELTAT/ ( LF*BF*TF ) 



85 



C_SL=VK/LF**o. 5 

C_BT=BF/TF 

C_LB=LF/BF 



IF ( NSHAFT .EQ. 1 ) THEN 

IF ( C_SL .GT. 1.20 .OR. C_SL . LT. 0.40 ) THEN 

PRINT *,'*******vk/SQRT(L) ratio OUT OF RANGE ( 1) ******** 
GOTO 30 

ELSE 

CONTINUE 
END IF 

IF ( CB .GT. 0.86 .OR. CB . LT. 0.50 ) THEN 

PRINT */,v******bLOCK COEFFICIENT OUT OF RANGE ( 1 )******* 
GOTO 30 

ELSE 

CONTINUE 

ENDIF 

IF ( C_BT . GT. 4. 50 . OR. C_BT . LT. 2. 0 ) THEN 

PRINT *, '********BE AM/ DRAFT RATIO OUT OF RANGE ( 1 )******* 
GOTO 30 

ELSE 

CONTINUE 

ENDIF 

IF ( C_LB .GT. 9.50 .OR. C_LB . LT. 3.33 ) THEN 

PRINT *,'********LENGTH/BEAM RATIO OUT OF BOUNDS( 1 )***** 
GOTO 30 

ELSE 



CONTINUE 

ENDIF 

ELSEIF ( NSHAFT . EQ. 2 ) THEN 

IF ( C_SL .GT. 1.20 .OR. C_SL . LT. 0.40 ) THEN 

PRINT *, ’*******VK/SQRT(L) RATIO OUT OF RANGE(2)* 
GOTO 30 



ELSE 

CONTINUE 

ENDIF 

IF ( CB . GT. 0. 80 . OR. CB . LT. 0. 54 ) THEN 

PRINT * , ’ *******bLOCK coefficient OUT OF RANGE(2)**** 
GOTO 30 

ELSE 

CONTINUE 

ENDIF 

IF ( C_BT .GT. 4.50 .OR. C_BT . LT. 2.0 ) THEN 

PRINT * , ' ********BEAM/DRAFT RATIO OUT OF RANGE(2)**** 
GOTO 30 

ELSE 

CONTINUE 

ENDIF 

IF ( C_LB .GT. 11.50 .OR. C_LB . LT. 3.80 ) THEN 

PRINT *,'********LENGTH/BEAM RATIO OUT OF RANGE(2)***' 
GOTO 30 

ELSE 












CONTINUE 

ENDIF 

ELSE 

PRINT *, ’********nuMBER of SHAFTS OUT OF RANGE( 1 , 2)******' 



o o o o 



GOTO 30 
ENDIF 
C 

WRITE (8 141) 

PRINT */ BLOCK COEFFICIENT = ' ,CB 

WRITE (8,15) CB 

PRINT SPEED-LENGTH RATIO = ' ,C_SL 

WRITE (8,16) C_SL 

PRINT BEAM-DRAFT RATIO = ' ,C_BT 

WRITE (8,17) C_BT 

PRINT LENGTH-BEAM RATIO = ' ,C_LB 

WRITE (8,18) C_LB 
WRITE (8,141) 

C 

WRITE (8,7) 

C 

C*****NOW CALCULATE BY SILVERLEAF & DAWSON METHOD 
C 

PRINT STARTING S & D CALCULATION' 

WRITE (8,14) 

C 

C_VB=( 1.70-1. 4*CB)*LF**0. 5 
IF ( LF .GT, 1000.0 ) C_X=0. 85 

IF ( LF .LE. 1000.0 ) C_X=0. 85+0. 00185*(( 1000. 0-LF)/100. 0)**2. 5 
IF ( BF/TF . LT. 2. 40 ) C_KBT=0. 982 
IF ( BF/TF .GE. 2.40 ) C_KBT=0. 96+ 

& 0. 00054*10. 0**((2. 0*BF)/(3. 0*TF)) 

IF ( VK/C_VB .LE. 1.0 ) C_KV=2. 75-7. 25*( VK/C_VB) + 

& 5. 50*((VK/C_VB)**2. 0) 

IF ( VK/C_VB .GT. 1.0 ) C_KV=21. 2-43. 2*(VK/C_VB)+ 

& 23. 0*((VK/C_VB)**2. 0) 

IF ( NSHAFT .EQ. 1 ) C_H400=2. 60-0. 2917*VK/(DELTAT**( 1. 0/6. 0) ) 

IF ( NSHAFT .EQ. 2 ) C_H400=2. 38-0. 2917*VK/(DELTAT**( 1. 0/6. 0) ) 

IF ( LF . EQ. 400. 0 ) C_KLBP=1. 0 

IF ( LF .NE. 400.0 ) C_KLBP=0. 9196+2. 31E-04*LF-7. 50E-08*LF**2. 0 

IF ( NSHAFT .EQ. 1 ) ETA0120=0. 98-0. 55*CB 

IF ( NSHAFT .EQ. 2 ) ETA0120=0. 90-0. 33*CB 

IF ( NSHAFT .EQ. 1 ) DETA0=0. 360-0. 0029*NRPM 

IF ( NSHAFT .EQ. 2 ) DETA0=0. 135-0. 001 1*NRPM 

C_ETA0=ETA0120+DETA0 

*****CALCULATE S & D POWER NUMERATOR & DENOMINATOR FOR 
NORMAL AND 'CORRECTED' VERSIONS 

T0P_SD=C_X*(DELTAT**(2. 0/3. 0) )*( C_VB**3. 0)*C_KBT*C_KV 
B0T_SD=427. 10*C_ETA0*C_H400*C_KLBP 
C*****VALUE 

PHP_SD =TOP_SD/BOT_SD 
C*****OUTPUT VALUES 

PRINT *,' POWER IN HP BY S&D METHOD FOR NO B/B = ' ,PHP_SD 

WRITE (8,19) PHP_SD 

C*****CORRECT IF WELL-DESIGNED BULBOUS BOW IS PRESENT & OUTPUT 
PHPISD = 0. 95*PHP_SD 

PRINT *,' POWER IN HP BY S&D METHOD FOR B/B = ', PHPISD 

WRITE (8,20) PHPISD 
WRITE (8,7) 
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non 



WRITE (8,91) 

WRITE (8,23) 

WRITE (8,24) LF,BF,TF 
WRITE (8,25) 

WRITE (8,26) DELTAT,N,NRPM 
WRITE (8,27) 

WRITE (8,271) VK,C_VB,NSHAFT 
WRITE (8,28) 

WRITE (8,28) 

C 

C*****SET 'BEGIN AGAIN' CONDITIONS TO EXIT OR CONTINUE 
C 

29 PRINT *,' ' 

30 CONTINUE 

PRINT *,' DO YOU WANT TO PROCESS AGAIN(??): INSTRUCTIONS: ' 

PRINT ENTER ZERO (0) TO TERMINATE OR ONE (1) TO PROCEED: ' 

PRINT*,' RESPONSE: (???) : ' 

READ *,IANSWER 
IF ( lANSWER .EQ. 0 ) THEN 
CLOSE (UNIT=8) 

GOTO 31 

ELSEIF ( lANSWER . EQ. 1 ) THEN 
GOTO 1 

ELSE 

GOTO 29 
ENDIF 



*****FORMAT STATEMENTS 



7 

8 
9 

91 

92 
10 
11 
12 
141 

15 

16 

17 

18 
14 

19 

20 

23 

24 

25 

26 

27 
271 

28 
C 

31 



(2X, 

(2X, 

(lOX 

(2X, 

(lOX 

(2X, 



FORMAT(2X,' 
FORMAT (2X, 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT(2X,' 
FORMAT (2X, 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 
FORMAT 



(2X, 

(2X, 

(2X, 

(2X, 

(2X, 

(2X, 

(lOX 

(2X, 

(lOX 

(2X, 

(lOX 

(2X, 

(2X, 



') 

'ADMIRALTY CALCULATION') 

'power IN HORSEPOWER FOR ADMIRALTY METHOD ' ,F12. 2) 

'SHIP DESIGN CHARACTERISTICS') 

, ' LENGTH(M) DISPL(MTONS) ' ) 

2F19. 2) 

,'SPEED(M/S) PROP SPD(RPS)') 

2F19. 2) 

Vr V \ Vr ^ ^ 

'BLOCK COEFFICIENT ',F15.6) 

'SPEED LENGTH RATIO ' ,F15. 6) 

'BEAM/DRAFT RATIO ' ,F15. 6) 

'LENGTH/BEAM RATIO ' ,F15. 6) 

'SILVERLEAF & DAWSON CALCULATIONS') 

'POWER IN HP BY S&D METHOD FOR NO B/B ' ,F15. 2) 

'POWER IN HP BY S&D METHOD FOR WD B/B ' ,F15. 2) 

, ' LENGTH(FT) BEAM(FT) DRAFT(FT) ' ) 

3F19. 2) 

, 'DISPL(TONS) SHAFT RPS SHAFT RPM') 

3F19. 2) 

,' DESIGN SPD(KTS) BOUN SPD(KTS) # SHAFTS') 

2F19. 2,13X,I3) 

■) 



CONTINUE 

STOP 

END 
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B. TEST ROUTINE OUTPUT 



This output was used to verify the algorithms and is used in part in Table 1 of 
Chapter 3. 



■ LARGE GENERAL CARGO Design 






ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 20351.91 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

177.59 31491.14 

SPEED(M/S) PROP SPD(RPS) 

10.70 2.00 



BLOCK COEFFICIENT 0. 669842 

SPEED LENGTH RATIO 0. 861818 

BEAM/DRAFT RATIO 2.342857 

LENGTH/BEAM RATIO 7. 103659 



Vc Vc Vc Vc Vc Vc Vr Vc 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY SGD METHOD FOR NO B/B 23780. 55 

POWER IN HP BY S8D METHOD FOR WD B/B 22591.52 



SHIP DESIGN CHARACTERISTICS 

LENGTH(FT) BEAM(FT) DRAFT(FT) 

582.50 82.00 35.00 

DISPL(TONS) SHAFT RPS SHAFT RPM 

31995.00 2.00 120.00 

DESIGN SPD(KTS) BOON SPD(KTS) # SHAFTS 

20.80 18.40 1 






• CONTAINER A Design 









*****Vr ********** 



ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 14127.90 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

176.83 21732.28 

SPEED(M/S) PROP SPD(RPS) 

10.29 2.00 
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BLOCK COEFFICIENT 0. 542293 

SPEED LENGTH RATIO 0. 830455 

BEAM/DRAFT RATIO 2.476191 

LENGTH/BEAM RATIO 7.435897 






SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 11503. 22 

POWER IN HP BY S&D METHOD FOR WD B/B 10928. 06 

SHIP DESIGN CHARACTERISTICS 

LENGTH(FT) BEAM(FT) DRAFT(FT) 

580.00 78.00 31.50 

DISPL(TONS) SHAFT RPS SHAFT RPM 

22080. 00 2. 00 120. 00 

DESIGN SPD(KTS) BOUN SPD(KTS) # SHAFTS 

20. 00 22. 66 1 

Vf*’)V*”jV*Vc”jWc’>V*'}V*9V*’5WfVr’jWf9WcVr>V*’5V*'5WrsWf’>V’jV^V*’)V*Vf9V*’)V*’jV**’jV*yc”>V'jV 



■ CONTAINER B Design 



')Wc”jWr5Wc’Vc”5V'>V’}V*’jWrjV^V'}V’)V^V'>V“/V*^V'jWf^V?WoV'jV'5V?V'5V'5Wf’jV’jV'jW«’VfVr>V'jV'jV^V**^V'5V^Wf 

ADMIRALTY CALCULATION 

POWER IN HORSEPOWER , FOR ADMIRALTY METHOD 30511.80 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

206. 40 38090. 55 

SPEED(M/S) PROP SPD(RPS) 

11.73 2.00 



*** ** ir VfVr *>V x *}V'>Wf * *^VVf 



BLOCK COEFFICIENT 0. 619424 

SPEED LENGTH RATIO 0. 876275 

BEAM/DRAFT RATIO 2. 794118 

LENGTH/BEAM RATIO 7. 126316 



■jV'jV^Wr 'jV’jV'jV’jV * Vr 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 29480.57 

POWER IN HP BY S&D METHOD FOR WD B/B 28006.54 

SHIP DESIGN CHARACTERISTICS 

LENGTHC FT) BEAM( FT) DRAFT( FT) 

677.00 95.00 34.00 

DISPL(TONS) SHAFT RPS SHAFT RPM 

38700.00 2.00 120.00 

DESIGN SPD(KTS) BOUN SPD(KTS) # SHAFTS 

22.80 21.67 1 



• ROLL ON/ROLL OFF Design 
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^V^Wf5V*'>V^V*>Wf'5WcVf'jV*5Wc*5V*>V*******9V^V:^'5V^V^V**9V**'5V****9V******Vr** 

ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 28573.58 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

195.12 33233.27 

SPEED(M/S) PROP SPD(RPS) 

11.83 2.00 



* ic ic * * ic * Vf ****** * * * * 



BLOCK COEFFICIENT 0.565724 

SPEED LENGTH RATIO 0. 909155 

BEAM/DRAFT RATIO 3. 187500 

LENGTH/BEAM RATIO 6.274510 



***************** 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 27435.27 

POWER IN HP BY S&D METHOD FOR WD B/B 26063.51 



SHIP DESIGN CHARACTERISTICS 

LENGTH(FT) BEAM(FT). . . 

640. 00 102. 00 

DISPL(TONS) SHAFT RPS. . . 

33765.00 2.00 

DESIGN SPD(KTS) BOON SPD(KTS) 

23. 00 22. 97 






DRAFT(FT) 

32. 00 

SHAFT RPM 

120. 00 

// SHAFTS 

1 

' ******* 



Vc Vc Vv Vc Vc Vc W Vc Vc W Vv V? Vc Vc VcVc Vc Vc Vc VcVcVc Vc VcVc Vc Vc Vc VcV'T Vc^Wc 



• LASH BARGE CARRIER Design 



ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 26203. 87 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

220.73 32135.83 

SPEED(M/S) PROP SPD(RPS) 



11. 57 


2. 00 


•id’s **** * ***Vr *** * ** 




BLOCK COEFFICIENT 


0. 563709 


SPEED LENGTH RATIO 


0. 836206 


BEAM/DRAFT RATIO 


3. 571429 


LENGTH/ BEAM RATIO 


7. 240000 



Vc Vc Vc Vr Vc V c Vc V c Vc Vc 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 23787.66 

POWER IN HP BY S&D METHOD FOR WD B/B 22598. 28 



SHIP DESIGN CHARACTERISTICS 
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LENGTHC FT) BEAM( FT) DRAFT( FT) 

724.00 100.00 28.00 

DISPL(TONS) SHAFT RPS SHAFT RPM 

32650.00 2.00 120.00 

DESIGN SPD(KTS) BOUN SPD(KTS) # SHAFTS 

22.50 24,51 1 



« SEABEE BARGE CARRIER Design 



^*>V*jWcVfV«”>V*jV^VyoV^V^iV^V*jWf*)Wf*'}Wr*>VVf^:V^iVi'Wf'>Wc'jWnWfVr^:V‘>V’>V:V^*>V^:V:V:V:V:^r:V 

ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 26772. 95 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

220. 34 56387. 79 

SPEED(M/S) PROP SPD(RPS) 

10. 29 2. 00 



VoV^VVf *>V>Wc*>V* 



BLOCK COEFFICIENT 0.670064 

SPEED LENGTH RATIO 0. 743962 

BEAM/DRAFT RATIO 2. 708440 

LENGTH/BEAM RATIO 6. 824363 



Vf>VVfVcVc*}V">V*Vr/W->V*}V*jV:tV*‘V 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 21918. 27 

POWER IN HP BY S&D METHOD FOR WD B/B 20822. 36 



SHIP DESIGN CHARACTERISTICS 

LENGTHC FT) BEAM(FT) DRAFT(FT) 

722. 70 105. 90 39. 10 

DISPL(TONS) SHAFT RPS SHAFT RPM 

57290.00 2.00 120.00 

DESIGN SPD(KTS) BOUN SPD(KTS) # SHAFTS 

20.00 20.48 1 



Vc V? >V Vf Vr Vc ' 



:Wf'jWfVfVr>V’5V';V“WcVcVfVr:W«"sV“V’sV’>Wc*jWfVc’5V*)WcVcVr*>WcVr**)Wf*>VVfVv*>V‘5V 



*Vf*VcVc*jV>VVfVr*>V*>V‘)VVfV-'sV^VV-'>VVcVf'>VVfVfV(rVcV«’*>V/VVfV-*>V?V’5V:V>V*iV>r*>V*i>VVr*>V*'3V*>V:^’V«’^V^V 



« TANKER A Design 



V\ V* Vr Vc Vc Vr Vc Vc Vc Vc ' 



* Vc Vc Vc Vc Vc Vc Vr Vc Vc Vc Vc Vc Vc Vc Vc Vv Vc Vc 



ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 12045. 93 



SHIP DESIGN CHARACTERISTICS 

LENGTHC M) DISPLCMTONS) 

201.22 46536.41 

SPEEDCM/S) PROP SPDCRPS) 

8.23 2.00 



Vr «VVc Vr Vc /VVc V; 'fCici 



92 



BLOCK COEFFICIENT 0.795976 

SPEED LENGTH RATIO 0. 622799 

BEAM/DRAFT RATIO 2.571429 

LENGTH/BEAM RATIO 7. 333333 



i c Vr ‘V Vc * Vr Vc Vr Vr ‘V Vf ^ V Vr Vf * Vc * 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 11077.69 

POWER IN HP BY S&D METHOD FOR WD B/B 10523. 80 

SHIP DESIGN CHARACTERISTICS 

LENGTH(FT) BEAM(FT) DRAFT(FT) 

660.00 90.00 35.00 

DISPL(TONS) SHAFT RPS SHAFT RPM 

47281.00 2.00 120.00 

DESIGN SPD(KTS) BOON SPD(KTS) # SHAFTS 

16.00 15.05 1 



■ TANKER B Design 



'>V*^VV-Vc’**VrVrVnVVr*:V';VVc*5V**5V*Vr***3VVc'*5VVfVr*->V**>V*jV***->V**Vc*‘5V'jVVfycVf*)V*^V'5V 

ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 53068. 73 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

348.48 443809.00 

SPEED(M/S) PROP SPD(RPS) 

8.18 2.00 



*jV VrVr VrVc’jV'jWf Vf-sV ' jVVc'sV *“V 



BLOCK COEFFICIENT 0.818361 

SPEED LENGTH RATIO 0. 470299 

BEAM/DRAFT RATIO 3.081081 

LENGTH/ BEAM RATIO 5.013158 



^V*5Wc Vf Vc* ‘sWt'jV'jVVc'A’* “V 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 36030.91 

POWER IN HP BY S&D METHOD FOR WD B/B 34229. 36 

SHIP DESIGN CHARACTERISTICS 

LENGTH(FT) BEAM(FT) DRAFT(FT) 

1143.00 228.00 74.00 

DISPL(TONS) SHAFT RPS SHAFT RPM 

450910.00 2.00 120.00 

DESIGN SPD(KTS) BOUN SPD(KTS) SHAFTS 

15.90 18.74 1 

’5WnWfVr'>Wc’'5V'iW<r^V')WfVc’*’5Wc'5V:Wc’*'5WrtWr^V^V?V^Wc’';V*‘5V^Wc*:fV‘5V'5V^V:iV^Wc'5VV^’5V^V?V:iV'5V 

Vc'5Wr^Wr:Wr)Wc*'5Wr*****VfV«V*Vc*Vr:Wr‘jV’>V*^V*Vr>Wc’*5Wc*^V:V*Vc-***5V*Vf*“Wr* 



« LNG TANKER Design 
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*'>WoWc*'>WoWcVf****')Wo’oWc‘>V'>V')WcVc’***VfVc*>V****iV*VrjWc*yf>V'?V***VciWf 

ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 39727. 77 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

273.48 93110.23 

SPEED(M/S) PROP SPD(RPS) 

10. 49 2. 00 



* Vc *** Vf * Vf Vf ***** 



BLOCK COEFFICIENT 0.717015 

SPEED LENGTH RATIO 0.681136 

BEAM/DRAFT RATIO 3. 972222 

LENGTH/BEAM RATIO 6.272727 



*********Vr**Vr**** 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 37910.40 

POWER IN HP BY S&D METHOD FOR WD B/B 36014. 88 



SHIP DESIGN CHARACTERISTICS 

LENGTH(FT) BEAM(FT) DRAFT(FT) 

897.00 143.00 36.00 

DISPL(TONS) SHAFT RPS SHAFT RPM 

94600.00 2.00 120.00 

DESIGN SPD(KTS) BOUN SPD(KTS) # SHAFTS 

20. 40 20. 85 1 

y- yj. y- y- 

****Vc****Vc’***********************VfVc*************** 



■ BULK CARRIER Design 



ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 10940. 75 



SHIP DESIGN CHARACTERISTICS 

LENGTH(M) DISPL(MTONS) 

178.05 31594.49 

SPEED(M/S) PROP SPD(RPS) 



8. 69 



2. 00 



* * * * * * * * * * * * * * * * * 

BLOCK COEFFICIENT 0.645051 

SPEED LENGTH RATIO 0.699327 

BEAM/DRAFT RATIO 2.912500 

LENGTH/BEAM RATIO 6. 266095 

* * * * * *- * * * * * * * * * Vc * 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 8578.57 

POW’ER IN HP BY S&D METHOD FOR WD B/B 8149.64 



SHIP DESIGN CHARACTERISTICS 
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LENGTHC FT) BEAM( FT) DRAFT( FT) 

584.00 93.20 32.00 

DISPL(TONS) SHAFT RPS SHAFT RPM 

32100.00 2.00 120.00 

DESIGN SPD(KTS) BOUN SPD(KTS) # SHAFTS 

16.90 19.26 1 



■ ORE/BULK/OIL CARRIER Design 
ADMIRALTY CALCULATION 

POWER IN HORSEPOWER FOR ADMIRALTY METHOD 21700.58 



SHIP DESIGN CHARACTERISTICS 

LENGTHC M) DISPL(MTONS) 

260.67 97647.63 

SPEEDCM/S) PROP SPD(RPS) 

8.49 2.00 



* ^ V Vf * * * * - jV * * * * * * * * * 



BLOCK COEFFICIENT 0.838120 

SPEED LENGTH RATIO 0.564288 

BEAM/DRAFT RATIO 2. 310044 

LENGTH/BEAM RATIO 8.081285 



Vf * * * * it it it it it it it it i t it it 



SILVERLEAF & DAWSON CALCULATIONS 

POWER IN HP BY S&D METHOD FOR NO B/B 18981. 25 

POWER IN HP BY S&D METHOD FOR WD B/B 18032. 19 

SHIP DESIGN CHARACTERISTICS 

LENGTHC FT) BEAMCFT) DRAFTCFT) 

855.00 105.80 45.80 

DISPLCTONS) SHAFT RPS SHAFT RPM 

99210.00 2.00 120.00 

DESIGN SPDCKTS) BOUN SPDCKTS) # SHAFTS 

16.50 15.40 1 



it it it it it it it it it it i t it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it it 
itititititititititititititititititititititititititititititititititititititititititititititititititit 
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APPENDIX D. POWER PREDICTION REPORT MODULE SOURCE 

CODE 



This is the source code for the subroutines that were added to the Report Module 
to produce the Power Prediction Reports. 



SUBROUTINE POWER_RPT 



c * 

C POWER REPORT SECTION READS IN A DATA FILE CREATED BY THE * 
C POWER PREDICTION MODULE OF THE MAIN PROGRAM AND WRITES A * 
C DETAILED HARD COPY REPORT OF THE RESULTS. THE FIRST PART OF * 
C THE REPORT ECHOS THE INPUT DATA AND THE CALCULATES THE * 
C RESULTS FOR THE REPORT. * 
C * 
C CALLED BY SUBROUTINE 'REPORT' * 
C CALLS SUBROUTINE ’ PW CALC 2' * 









INCLUDE 'TOP. FOR' 

INCLUDE ' TOP_POWER. FOR ' 

CHARACTER PW_INPUT_DATA*12 , PW_PRFILE*12 
C*****MOVE THE MOUSE POINTER AWAY FROM THE SELECTION 

STA = UIS$SET_POINTER_POSITION(VD_ID,WD_RPT,9. 9,23. 5) 
C*****SET UP ’THE FIRST SCREEN INSTRUCTIONS 



CALL UIS$ERASE(VD_ID, 4.1,20.1,8.4,24.9) 

CALL UIS$TEXT(VD_ID,7, 'ENTER THE FILE NAME ' ,4. 5 , 24. ) 
CALL UIS$TEXT(VD_ID,7, 'AND EXTENSION OF THE' ,4. 5,23. 6) 
CALL UIS$TEXT(VD_ID,7, 'DATA FILE TO BE USED' ,4. 5 , 23. 2) 
CALL UIS$TEXT(VD_ID,7, 'OR [RETURN] TO EXIT' ,4. 5 ,22. 8) 

c***,wcread the data input file name 



CALL KEY_READ (PW_INPUT_DATA, . FALSE. ,5. ,21. 2, *30) 
C*****NOW WE NEED THE NAME OF THE OUTPUT REPORT FILE 
CALL UIS$ERASE(VD_ID, 4.1,20.1,8.4,24.9) 

CALL UIS$TEXT(VD_ID,7, ' ENTER THE FILE NAME ' ,4. 5 , 24. ) 
CALL UIS$TEXT(VD_ID,7, ' AND EXTENSION TO BE ',4.5,23.6) 
CALL UIS$TEXT(VD_ID, 7, 'ASSIGNED TO THE REPORT' , 4. 5 , 23. 2) 
CALL UIS$TEXT(VD_ID,7, ' OR [RETURN] TO EXIT' ,4. 5 , 22. 8) 
CALL KEY_READ ( PW_PRFILE , . FALSE. ,5. ,21. 2, *30) 

C*****READ THE DATA IN 

1 OPEN (21, FILE = PW_INPUT_DATA, STATUS = 'OLD') 

OPEN (22, FILE = PW_PRFILE, STATUS = ' NEW' ,BLANK=' NULL' ) 

READ (21,100) PW_INFILE 

READ (21,110) PW_LPP 

READ (21,110) PW_T 

READ (21,110) PW_B 

READ (21,110) PW_DISP 

READ (21,110) PW_CB 

READ (21,110) PW_VK 

READ (21,110) PW_NRPM 
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o o 



READ (21,110) PW.NSHAFT 
READ (21,110) PW_ADME 
READ (21,110) PW.SDNB 
READ (21,110) PW_SDBB 
REWIND(21) 

CLOSE (21) 

C*****RE -CALCULATE FOR ACCURACY 
CALL PW_CALC_2 

PW_SLR = PW_VK/(PW_LPP**0. 5) 

PW_BTR = PW_B/PW_T 
PW_LBR = PW_LPP/PW_B 

c***,v*generate the cover page of the report 



DO 5 I = 1,18,1 

WRITE (22,130) ' ' 

5 CONTINUE 

WRITE (22,120) ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' 

& REPORT' 

WRITE (22,120) ' 

WRITE (22,120) ' 

&PROGRAM' 

WRITE (22,120) ' ' 

WRITE (22,120) ' 

6eNAVAL ' 

WRITE (22,120) ' 

WRITE (22,120) ' 

&THE ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' 

WRITE (22,120) ' 

WRITE (22,120) ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' ' 

WRITE (22,120) ' ' 

DO 6 I = 1,13,1 

WRITE (22,120) ' ' 

6 CONTINUE 

WRITE (22,120) ' ' 

WRITE (22,120) ' PAGE 1 OF 
WRITE (22,120) ' ' 

C*****IN THIS SECTION, THE ECHO ' 
DATA IS PRESENTED. 



TOOL BOX ' 

POWER PREDICTION 

THIS REPORT WAS GENERATED USING THE 
TOOL BOX WHICH WAS DEVELOPED FOR THE 

I 

ENGINEERING DEPARTMENT OF 

NAVAL POSTGRADUATE SCHOOL’ 
MONTEREY, CALIFORNIA' 

PROFESSOR F, PAPOULIAS ' 
AND ' 

LT. GERALD MCGOWAN ' 

LT. JAMES PLOSAY' 

1989/90 ' 



/ POWER prediction' 

SHIP PARAMETERS AND CALCULATED 



DO 10 I = 1,5,1 

WRITE (22,130) ' ' 
10 CONTINUE 
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POWER PREDICTION REPORT ' 



WRITE(22,120) 

WRITE(22,130) 

WRITE(22,120) 

WRITE(22,130) 

WRITE(22,150) 

WRITE(22,150) 

& PW_INPUT_DATA 
WRITE(22,160) 
WRITE(22,130) 
WRITE(22,150) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,160) 
WRITE(22,140) 
WRITE(22,140) 
WRITE(22,140) 
WRITE(22,140) 
WRITE(22, 160) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,130) 
WRITE(22,130) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,130) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,140) 
WRITE(22,130) 
WRITE(22,130) 
WRITE(22,130) 
WRITE(22,130) 
&NG: ’ 

WRITE(22,130) 

WRITE(22,130) 

&F' 

WRITE(22,130) 

&83' 

WRITE(22,130) 

& BY S. ' 
WRITE(22,130) 
WRITE(22,130) 



t 



THE INPUT SHIP PARAMETERS ARE AS FOLLOWS' 

I 

THE REPORT IS LOCATED IN FILE : ' ,PW_PRFILE 

THE INPUT DATA FILE USED IS : ' , 



DESIGN PARAMETERS' 

SHIP NAME IS ',PW_INFILE 

I 

LENGTH BETWEEN PERPENDICULARS FT ' ,PW_LPP 

DpiGN DRAFT FT ' ,PW_T 

DpiGN BEAM FT ' ,PW_B 

DpiGN DISPLACEMENT LTONS. . ' ,PW_DISP 

COEFFICIENTS OF FORM*******^’^****’^******* * 

BLOCK COEFFICIENT: ' ,PW_CB 
SPEED-LENGTH RATIO: ' ,PW_SLR 
BEAM-DRAFT RATIO: ' ,PW_BTR 

LENGTH -BEAM RATIO: ' ,PW_LBR 

*^V^WcVfVr^V^Wc'*'jVVf**^V^V^VVr^V^Wf^V*Vr^V^Wf^V^Wc^V^V^Wc'jV^V^Wc’^V'jV * 



DESIGN OPERATING SPEED (NOM) KNOTS. .' ,PW_VK 

f 

PROPULSION SHAFT SPEED RPM. . ' ,PW_NRPM 

t 

NUMBER OF PROPULSION SHAFTS ' ,PW_NSHAFT 

I 

Vr^VVc'5Vy*^VVfVc?V“VVr<V'iVVrVc*^V*5V*VVfVf*yrVr‘)VVv'jV^V*)V^VV?‘)VVcyrVc‘5V*'5VVfVcVc?VVcVr'iV‘5VVr'5V'3V * 



DESIGN RESULTS (HORSEPOWER)' 

I 

ADMIRALTY POWER ' ,PW_ADME 

I 

SILVERLEAF & DAWSON POWER: ' 

I 

DESIGN W/0 BULBOUS BOW ' ,PW_SDNB 

I 

DESIGN W/ BULBOUS BOW ' ,PW_SDBB 

t 



■sV**Vr*>V**5’f')V>V*5Wf'iViWc'*iV**tiWr')V****5V****'>Wc'*5Wf ***»>■)>******* ' 

POWER ESTIMATIONS PRESENTED ARE BASED UPON FOLLOWI 



(1) METHOD OF ADMIRALTY COEFFICIENTS' 

HARVALD, Sv. Aa. /'RESISTANCE AND PROPULSION 0 

SHIPS", JOHN WILEY & SONS, NEW YORK,N. Y. ; 19 

(2) METHOD OF SILVERLEAF AND DAWSON, AS MODIFIED 

ERICHSEN: ERICHSEN,S. , REPORT No. 123 , ' 
"OPTIMUM CAPACITY OF SHIPS AND PORT TERMINAL 
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&s"' 

WRITE(22,130) ' UNIV. OF MICHIGAN, ANN ARBOR, MI. ; 1971' 

WRITE( 22 , 130) ' ************************************************** ' 
WRITE (22,120) ' ' 

WRITE (22,120) ' PAGE 2 OF 2 / POWER PREDICTION' 

CLOSE (22) 

100 FORMAT (A) 

110 FORMAT (G12.5) 

120 FORMAT (10X,A) 

130 FORMAT (10X,A) 

135 FORMAT (lOX, A,' ' ,A) 

140 FORMAT (lOX, A,' ',G14.6) 

150 FORMAT (10X,2A) 

160 FORMAT (10X,A) 

c*****exit from the report and enable mouse 

30 CALL UIS$ERASE(VD_ID, 4.1,20.1,8.4,24.9) 

CALL UIS$TEXT(VD_ID,7,' SELECT AN OPTION WITH' ,4. 5 , 24. ) 

CALL UIS$TEXT(VD_ID,7, ' THE MOUSE' ,4. 5, 23. 6) 

RETURN 

END 



C 
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APPENDIX E. TOOL BOX ENDURANCE ESTIMATION 



A. SUBROUTINE ENDURANCE SOURCE CODE 

This program is the TOOL BOX Endurance Estimation module that is called from 
the Main Menu, and all of its associated subroutines. 



SUBROUTINE ENDURANCE 



c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



* 

THIS SUBROUTINE IS THE MAIN PART OF THE ENDURANCE CALCULATIONS. * 

IT SETS UP SCREENS AND CONTROLS ENDURANCE PROGRAM ACTIONS. * 

* 

IT CALCULATES ENDURANCE USING THE USN DDS9400-1 FORMAT 
IT IS CALLED FROM THE MAIN TOOL BOX MENU BY THE "ENDURANCE" * 

SELECTION USING THE MOUSE. IT REQUIRES THE FOLLOWING INPUTS FOR * 
OPERATION AFTER INITIALIZATION: * 

* 

ND_DISP FULL LOAD DISPLACEMENT (TONS) REAL * 

ND_PWR CRUISING POWER (HP) REAL * 

ND.NDPWR AVERAGE ENDURANCE POWER REQD FOR CRUISING * 

REAL * 

ND_ELEC CRUISING ELECTRIC LOAD (KW) REAL * 

ND_MFR MAIN PROPULSION FUEL RATE (LBS/HR) REAL * 

ND_EFR CRUISING ELECTRIC PLANT FUEL RATE (LBS/KW-HR) 

REAL * 

ND_OTH OTHER FUEL CONSUMPTION RATES REAL * 

ND_TPA FUEL STORAGE TAIL PIPE ALLOW (%) REAL * 

* 



USING THE FOLLOWING ITERATIVE VARIABLES: 

ND_RANGE ENDURANCE RANGE (NMILES) REAL 

ND_SPD ENDURANCE SPEED (KNOTS) REAL 

ND_FUEL VOYAGE FUEL CAPACITY (TONS) REAL 

AND CALCULATES THE FOLLOWING SOLUTION VALUES: 

ND_TIME TIME (HRS) FOR JOURNEY OF ND_RANGE AT ND_SPD 

REAL 

ND_VOL FUEL STORAGE CAPACITY REQD REAL 



ND_PCT FUEL PERCENT OF FULL LOAD DISPL REAL 

ALL PARAMETERS ARE PASSED IN COMMON BY THE INCLUDE ' TOP_ENDU. FOR' 
FILE THAT HAS ALL THE PARAMETER DEFINITIONS, TYPE STATEMENTS, 
AND SET UP. 

CALLED BY MAIN MODULE 'TOOL_BOX' 

CALLS SUBROUTINES 'ND_KB_DATA_IN' , ' ND_READ_FILE' ,ND_REC0RD' , 
'ND_CALC_1' , 'TO_MAIN' ,AND 'EXIT' 



* 

Vc 



Vc 

Vc 



Vc 

* 

Vc 

Vc 

* 

if 

if 

if 



if 



if 

if 

if 

if 

if 



if if if if if if if if if if if if if if if • 



if if if if if if if if if if if if if if ifif if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if if 



INCLUDE 'general. FOR' 

INCLUDE 'TOP_ENDU.FOR' 

EXTERNAL TO_MAIN , EX IT, DARK2 ,LIGHT2 , NOWHERE ,ND_CALC_2 
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EXTERNAL ND_KB_DATA_IN,ND_RECORD,ND_READ_FILE ,ND_CALC_1 
EXTERNAL ND_SDI_LABEL , ND_SORT , MUCH_TOO_MUCH , TOO_MUCH , ND_LI ST 
EXTERNAL KEY_READ,SET_HELP_LEVEL 
C*****SET HELP LEVEL FOR PARAMETER WINDOWS TO 'ON' UPON ENTRY 
ND_HELP=1 

C*****initIALIZE the parameters to avoid divide by ZERO ERRORS. ALSO 
THIS INSURES THAT WHEN WE INITIALLY CALL THE KEYBOARD DATA ENTRY 
ROUTINE THAT THE DATA PASSES W/0 DIVIDE BY ZERO ERRORS 
IF ( ND_ENTERS . EQ. 0 ) THEN 
C*****INPUT VARIABLES 
ND_SOLVE=l 
ND_INFILE=' GENERIC' 

ND_KB_UP(1)=' GENERIC' 

ND_DISP=1. 0 
ND_KB_UP(2) = ' 1. O' 

ND_PWR=0. 0001 
ND_KB_UP(3)='0. 0001' 

ND_NDPWR=0. 0001 
ND_KB_UP(4)='0. 0001' 

ND_ELEC=0. 0001 
ND_KB_UP(5) = ' 0. 0001' 

ND_TPA=0. 0001 
ND_KB_UP(6)='0. 0001' 

ND_MFR=0. 0001 
ND_KB_UP(7)='0. 0001' 

ND_EFR=0. 0001 
ND_KB_UP(8)='0. 0001' 

ND_OTH=0. 0001 
ND_KB_UP(9)='0. 0001' 

ND_RANGE=0. 0001 
ND_KB_UP(10) = '0. 0001' 

ND_SPD=0. 0001 
ND_KB_UP(11) = '0. 0001' 

ND_FUEL=0. 0001 
ND_KB_UP(12) = '0. 0001' 

ND_WGT=0. 0001 
ND_KB_UP(18)='0. 0001' 

ND_VOL=0. 0001 
ND_KB_UP(19)='0. 0001' 

ND_PCT=0. 0001 
ND_KB_UP(20) = '0. 0001' 

ENDIF 

ND_ENTERS=ND_ENTERS+1 
C*****ERASE THE DATA AREA 

CALL UIS$ERASE(VD_ID,3. 6,10. 1,9. 9,18. 5) 

G*****REMOVE THE MOUSE POINTER TO SOME OTHER AREA 

STA = UIS$SET_P0INTER_P0SITI0N(VD_ID,WD_MAIN,9. 9,4. 5) 
C*****CREATE POWER WINDOW 

WD_NDUR=UIS$CREATE_WINDOW(VD_ID, ' SYS$WORKSTATION' , 'ENDURANCE 
& PREDICTION WINDOW' ,-. 5,9. 9,10. 1,19. ,40. ,30. ) 

DO 10 Y_COOR = 13. 1,18. 4,. 8 
DY = Y_COOR + . 2 

CALL UIS$SET_P0INTER_AST(VD_ID,WD_NDUR,DARK2, ,XO,Y_COOR,X1, 

& DY,LIGHT2) 

10 CONTINUE 

c*****set menu titles 
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OPTION(ll) = OPTION(23) 

0PTI0N(12) = 0PTI0N(22) 

C*****wRITE SCREEN 

CALL UIS$TEXT(VD_ID,0,OPTION(11),. 3,15. 1) 

CALL UIS$TEXT(VD_ID,0,OPTION(12),. 3,15. 9) 

CALL UIS$SET_BUTTON_AST(VD_ID,WD_NDUR,EXIT, , ,X0,13. 0, 

& XI, 13. 4) 

CALL UIS$SET_BUTTON_AST(VD_ID,WD_NDUR,TO_MAIN, , ,X0, 13. 8, 

& XI, 14. 2) 

CALL UIS$SET_BUTT0N_AST(VD_ID,WD_NDUR,SET_HELP_LEVEL, , ,X0, 14. 6, 
& XI, 15.0) 

CALL UIS$SET_BUTT0N_AST(VD_ID,WD_NDUR,ND_CALC_1, , ,X0, 15. 4, 

& XI, 15. 8) 

CALL UIS$SET_BUTT0N_AST(VD_ID,WD_NDUR,ND_REC0RD, , ,X0, 16. 2, 

& XI, 16. 6) 

CALL UIS$SET_BUTT0N_AST(VD_ID,WD_NDUR,ND_READ_FILE, , ,X0 
& ,17. , XI, 17. 4) 

CALL UIS$SET_BUTTON_AST( VD_ID ,WD_NDUR, ND_KB_DATA_IN, 

& , ,X0,17. 8, XI, 18. 2) 

WNDOW = WD_NDUR 
20 RETURN 

END 



«v* 



SUBROUTINE ND_KB_DATA_IN 

' / Wf Vr V? Vr V? “ Wc Vc Vf Vf Vr Vf Vf Vf Vf Vp V* Vf Vf V* Vf Vf Vf Vf Vf Vr Vf Vf “Wf Vf Vf * W> y? 



* 

THE ND_KB_DATA_IN ROUTINES LOAD IN KB DATA AND ALLOWS * 
THE OPERATOR TO CHANGE INPUTS AND SEE REAL TIME EFFECTS * 
OF THOSE CHANGES ON THE ENDURANCE CHARACTERISTICS OF THE * 
PRELIMINARY DESIGN. * 






CALLED BY SUBROUTINE 'ENDURANCE' 

CALLS SUBROUTINES 'ND_SDI_LABEL' , 'KEY_READ' , 'ND_S0RT' , 
'ND_CALC_1', 'ND_LIST' 






* 

* 



* 



INCLUDE 'TOP. FOR' 

INCLUDE 'TOP_ENDU. FOR' 

REAL NT)_LINE_NO 

STA=UIS$SET_P0INTER_P0SITI0N(VD_ID,WD_NDUR,9. 9,14. 5) 
KB_ID=UIS$CREATE_KB( ' SYS$WORKSTATION' ) 

CALL UIS$ENABLE_KB(KB_ID,WD_NDUR) 

*,v***FiRST, WE WE WRITE UP THE LINE LABLES AND FIRST INSTRUCTIONS 
CALL ND_SDI_LABEL(VD_ID) 

*****sECOND, CALCULATE AND DISPLAY STARTING VALUES 
CALL ND_CALC_1 

*****THIRD, ALLOW USER TO CHANGE INPUT THEN 
*****recALCULATE RESULTS IN REAL TIME 
1 CALL UIS$ERASE(VD_ID,-0. 4,10. 1,3. 4,10. 9) 

CALL UIS$TEXT(VD_ID,7,'ENTER A LINE NUMBER' , 0 . 0 , 11. 0) 

CALL UIS$TEXT(VD_ID,7, 'OR [RETURN] TO EXIT' ,0. 0, 10. 6) 

3 CALL KEY_READ(ND_LINE,. TRUE. , 0 . 0 , 12 . 0 , * 300 ) 

READ (ND_LINE,FMT='(F2.0)' ,ERR = 5) ND_LINE_NO 
NLINE_NO = INT(ND_LINE_NO) 
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GO TO (10,20,30,40,50,60,70,80,90,100,110) NLINE_NO 
*****iF IT GETS TO HERE A MISTAKE HAS BEEN MADE ********************** 



5 CALL UIS$ERASE(VD_ID,-0. 4,10. 1,3. 4,10. 9) 

CALL UIS$TEXT(VD_ID, 7, 'IMPROPER LINE NUMBER' , 0. 0 . 11. 0) 

CALL UIS$TEXT(VD_ID,7, 'ENTER A NEW NUMBER PLEASE* ,0. 0, 10. 6) 
GOTO 3 
C 

10 CALL ND_S0RT(NLINE_N0,DEL_Y,ND_LINE,*300) 

ND_KB_UP(1) = ND_LINE 

ND_INFILE = ND_LINE 

READ (ND_LINE, FMT=' (G) ' .ERR = 5) 

WRITE (ND_KB_UP(1), FMT=*(A)') ND_INFILE 
CALL UIS$ERASE(VD_ID,8. 2 , 17. 9 , 9. 9 , 18. 3) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(1),8. 2,18. 2) 

GOTO 1 
C 

20 CALL ND_SORT(NLINE_NO,DEL_Y,ND_LINE,*300) 

ND_KB_UP(2) = ND_LINE 

READ (ND_LINE, FMT=' ( G) ' ,ERR = 5) ND_DISP 
WRITE (ND_KB_UP(2), FMT=' (G12. 5) ' ) ND_DISP 
CALL UIS$ERASE(VD_ID,8. 2, 17. 5,9. 9,17. 9) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(2),8. 2,17. 8) 

GOTO 1 
C 

30 ND_NDX=1 

IF ( ND_HELP .EQ. 1 ) CALL ND_LIST(ND_NDX) 

CALL ND_SORT( NLINE_NO , DEL_Y , ND_LINE , *300 ) 

ND_KB_UP(3) = ND_LINE 

READ (ND_LINE, FMT=' (G) ' ,ERR = 5) ND_PWR 
WRITE (ND_KB_UP(3), FMT=' (G12. 5) ' ) ND_PWR 
CALL UIS$ERASE(VD_ID,8. 2,17. 1,9. 9,17. 5) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(3),8. 2,17. 4) 

GOTO 1 
C 

40 ND_NDX=2 

IF ( ND_HELP .EQ. 1 ) CALL ND_LIST(ND_NDX) 

CALL ND_SORT( NLINE_NO , DEL_Y , ND_LINE , *300 ) 

ND_KB_UP(4) = ND_LINE 

READ (ND_LINE, FMT=' ( G) ' ,ERR = 5) ND_NDPWR 
WRITE (ND_KB_UP(4), FMT=' ( G12. 5 ) ' ) ND_NDPWR 
CALL UIS$ERASE(VD_ID,8. 2,16. 7,9. 9,17. 1) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(4),8. 2,17. 0) 

GOTO 1 
C 

50 ND_NDX=3 

IF ( ND_HELP .EQ. 1 ) CALL ND_LIST(ND_NDX) 

CALL ND_SORT( NLINE_NO , DEL_Y , ND_LINE ,*300 ) 

ND_KB_UP(5) = ND_LINE 

READ (ND_LINE, FMT=' ( G) ' ,ERR = 5) ND_ELEC 
WRITE (ND_KB_UP(5), FMT=' ( G12. 5 ) ' ) ND_ELEC 
CALL UIS$ERASE(VD_ID,8. 2,16. 3,9. 9,16. 7) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(5),8. 2,16. 6) 

GOTO 1 
C 
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60 ND_NDX=4 

IF ( ND_HELP .EQ. 1 ) CALL ND_LIST(ND_NDX) 
CALL ND_S0RT( NLINE_N0 , DEL_Y , ND_LINE , *300 ) 
ND_KB_UP(6) = ND_LINE 

READ (ND_LINE, FMT='(G)' ERR = 5) ND_TPA 
WRITE (ND_KB_UP(6), FMT=' (G12. 5) ' ) ND_TPA 
CALL UIS$ERASE(VD_ID,8. 2,15. 9,9. 9,16. 3) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(6),8. 2,16. 2) 
GOTO 1 
C 

70 ND_NDX=5 

IF ( ND_HELP .EQ. 1 ) CALL ND_LIST(ND_NDX) 
CALL ND_S0RT(NLINE_N0,DEL_Y,ND_LINE,*300) 
ND_KB_UP(7) = ND_LINE 

READ (ND_LINE, FMT=’(G)' ERR = 5) ND_MFR 
WRITE (ND_KB_UP(7), FMT=^G12. 5) ' ) ND_MFR 
CALL UIS$ERASE(VD_ID,8. 2,15. 5,9. 9,15. 9) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(7),8. 2,15. 8) 
GOTO 1 
C 

80 ND_NDX=6 

IF ( ND.HELP .EQ. 1 ) CALL ND_LIST(ND_NDX) 
CALL ND_SORT( NLINE_NO , DEL_Y , ND_LINE , *300 ) 
ND_KB_UP(8) = ND_LINE 

READ (ND_LINE, FMT='(G)' ERR = 5) ND_EFR 
WRITE (ND_KB_UP(8) , FMT=^G12. 5) ' ) ND_EFR 
CALL UIS$ERASE(VD_ID,8. 2,15. 1,9. 9,15. 5) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(8),8. 2,15. 4) 
GOTO 1 
C 

90 ND_NDX=7 

IF ( ND_HELP .EQ. 1 ) CALL ND_LIST(ND_NDX) 
CALL KD_SORT( NLINE_NO , DEL_Y , ND_LINE , *300 ) 
ND_KB_UP(9) = ND_LINE 

READ (ND_LINE, FMT=’ (G) ' .ERR = 5) ND_OTH 
WRITE (ND_KB_UP(9) , FMT=' (G12. 5) ' ) ND_OTH 
CALL UIS$ERASE(VD_ID,8. 2,14. 7,9. 9,15. 1) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(9) ,8. 2,15. 0) 
GOTO 1 
C 

100 CALL ND_SORT(NLINE_NO,DEL_Y,ND_LINE,*3O0) 
ND_KB_UP(10) = ND_LINE 

READ (ND_LINE, FHT=' (G) ' ,ERR = 5) ND_RANGE 
WRITE (ND_KB_UP(10), FMT=’ (G12. 5) ' ) ND_RANGE 
CALL UIS$ERASE(VD_ID,8. 2 , 14. 3 , 9. 9 , 14. 7) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(10),8. 2,14. 6) 
CALL ND_CALC_1 
GOTO 1 
C 

110 CALL ND_SORT(NLINE_NO,DEL_Y,ND_LINE,*3O0) 
ND_KB_UP(11) = ND_LINE 

READ (ND_LINE, FMT=' ( G) ’ ,ERR = 5) ND_SPD 
WRITE (ND_KB_UP(11) , FMT=’ (G12. 5) ' ) ND_SPD 
CALL UIS$ERASE(VD_ID,8. 2,13. 9,9. 9,14. 3) 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(11),8. 2,14. 2) 
CALL ND_CALC_1 



104 



GOTO 1 
C 
C 

300 CALL UIS$DISABLE_KB(KB_ID) 

CALL UIS$ERASE(VD_ID,-. 4,10. 1,3. 4,10. 9) 

CALL UIS$TEXT(VD_ID, 7, 'SELECT AN OPTION' ,0. ,11. ) 
CALL UIS$TEXT(VD_ID,7, ' WITH THE MOUSE ',0. ,10.6) 
RETURN 
END 
C 
C 

SUBROUTINE ND_SORT( NLINE_NO , DEL_Y , ND_LINE , COUNT ,* ) 



c * 

C SUBROUTINE ND_SORT SORTS OUT THE LINE NUMBER TO THE * 
C VARIABLE INVOLVED AND CALLS THE APPROPRIATE SCREEN * 
C INSTRUCTION TO BE WRITTEN IN THE INSTRUCTION BOX. * 
C * 
C CALLS SUBROUTINE 'KEY_READ' * 



c 

INCLUDE 'TOP. FOR' 

INCLUDE 'TOP_ENDU. FOR' 

REAL DEL_Y, DY, DOT_Y 
LOGICAL ND_FCHAR 
CHARACTER ND_KB_INST( 11)*25 
C*****SET THE SCREEN INSTRUCTIONS STRINGS 

ND_KB_INST(1) = 'ENTER NEW SHIP NAME ' 
ND_KB_INST(2) = 'ENTER FULL LOAD DISPL ' 
ND_KB_INST(3) = 'ENTER FULL RATED POWER ' 
ND_KB_INST(4) = 'ENTER ENDURANCE POWER ' 
ND_KB_INST(5) = 'ENTER NOMINAL ELEC LOAD ' 
ND_KB_INST(6) = 'ENTER TPA AS O.XY% ' 

ND_KB_INST(7) = 'ENTER FUEL RATE (LBS/HR)' 
ND_KB_INST(8) = 'ENTER FUEL USE(LBS/KwHR) ' 
ND_KB_INST(9) = 'ENTER FUEL USE (LBS/HR) ' 
ND_KB_INST(10) = 'ENTER ENDURANCE RANGE ' 

ND_KB_INST( 11) = 'ENTER ENDURANCE SPEED ' 

c*****erase the instruction area and find its data position 

CALL UIS$ERASE(VD_ID,-0.4,10. 1,3.4,10. 9) 

DEL_Y = 18. 6 - . 4 * NLINE_NO 
C****, '.-WRITE THE APPROPRIATE INSTRUCTION STRING 

CALL UIS$TEXT(VD_ID,7,ND_KB_INST(NLINE_N0) ,0. 0,11. 0) 
CALL UIS$TEXT(VD_ID,7,' [RETURN] TO EXIT' , 0. 0 , 10. 6) 
C*****DECIDE IF WE WANT THE FIRST LINE (NAME) OR A VALUE 
IF (NLINE_NO .EQ. 1) THEN 
ND_FCHAR = .FALSE. 

ELSE 

ND_FCHAR = .TRUE. 

ENDIF 

CALL KEY_READ(ND_LINE,ND_FCHAR,0. 0,12. 0, *300) 

DY = DEL_Y - . 35 

CALL UIS$ERASE(VD_ID,8. 2 , DY, 9. 9 ,DEL_Y) 

CALL UIS$TEXT(VD_ID,7,ND_LINE,8. 2,DEL_Y) 

RETURN 

300 RETURN 1 
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END 

C 

C 

SUBROUTINE ND_SDI_LABEL( VD_ID) 



C * 
C PERFORM DISPLAY VARIABLES IS A ROUTINE WHICH READS UP THE * 
C DATA ONTO THE DISPLAY AREA ALONG WITH THE SCREEN LABELS. * 
C * 
C CALLED BY SUBROUTINES ' ENDURANCE' ND_KB_DATA_IN' * 






c 



10 



20 



IMPLICIT INTEGER (A-Z) 

INCLUDE 'TOP_ENDU. FOR' 

INCLUDE ' SYS$LIBRARY: UISENTRY ' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

REAL DEL_Y, DOT_Y 

CHARACTER N_LINE_LBL( 20)*35 

CALL UIS$ERASE(VD_ID,4. ,10. 1,9. 9,18. 3) 



N_LINE_LBL(1) 
N_LINE_LBL(2) 
N_LINE_LBL(3) 
N_LINE_LBL(4) 
N_LINE_LBL(5) 
N_LINE_LBL(6) 
N_LINE_LBL(7) 
N_LINE_LBL(8) 
N_LINE_LBL(9) 
N_LINE_LBL(10) 
N_LINE_LBL(11) 
N_LINE_LBL(12) 
N_LINE_LBL(13) 
N_LINE_LBL(14) 
N_LINE_LBL(15) 
N_LINE_LBL( 16) 
N_LINE_LBL(17) 
N_LINE_LBL(18) 
N_LINE_LBL(19) 
N_LINE_LBL(20) 
DO 10 I = 1,20, 
DEL Y = 18. 6 



— I 



_ t 

_ I 



' 1 . 

' 2 . 

' 3. 
'4. 
'5. 
6 . 

7. 

8 . 

9. 

10 . 
11 . 



SHIP NAME 

FULL LOAD DISPLACEMENT(TONS) 
FULL RATED PLANT POWER (SHP) 
AVG ENDURANCE POWER (SHP).. 

CRUISING ELEC LOAD (KW) 

TAIL PIPE ALLOWANCE (%). 

CRUISING MAIN FUEL RATE 

CRUISING ELEC FUEL RATE 

CRUISING OTHER FUEL RATE. . . 
*E ST. ENDURANCE RANGE (MILES) 
*EST. ENDURANCE SPEED (KNOTS) 



_ I t 



_ I 



**RESULTS FROM CALC**' 

I 



J. JIm 



_ ! 

_ ! 



EST. ENDURANCE FUEL LOAD (TONS)' 

ELAPSED TRIP TIME (HRS); ' 

FUEL STORAGE RE13l)1REMENTS(FT 
' FUEL WEIGHT ALLOWANCE (%-FL): . . ' 

' *7V*>V****5V'***?V7Wf*****Vr**^V5V****Vc*i> ' 

- .4*1 



CALL UIS$TEXT(VD_ID,7,N_LINE_LBL(I),4. ,DEL_Y) 
CONTINUE 
DO 20 I = 1,9,1 



DEL_Y = 18. 6 - . 4*1 

CALL UIS$TEXT( VD_ID , 7 ,ND_KB_UP( I ) , 8. 2 ,DEL_Y) 
CONTINUE 
RETURN 
END 



C 

C 



SUBROUTINE ND_CALC_1 

C * 
C ND_CALC_1 CALLS ND_CALC_2 WHICH DOES THE ACTUAL ENDUAPvNCE * 
C CALCULATIONS. ON RETURN, ND_CALC_1 DISPLAYS THE RESULTS * 
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C OF THE CALCULATIONS AND THEN RETURNS TO THE ROOT PROGRAM. * 
C * 
C CALLED BY SUBROUTINES ' ENDURANCE ND_KB_DATA_IN' * 
C CALLS SUBROUTINE 'ND_CALC_2' * 



c 



10 



FOR' 



ND. 

ND 



RANGE 

SPD 



INCLUDE 'TOP. FOR' 

INCLUDE 'TOP_ENDU. 

REAL DY 

CALL UIS$ERASE(VD_ID,8. 1 , 9. 95 ,9. 9 , 14. 7) 

CALL ND_CALC_2 

WRITE( ND_KB_UP(10), FMT=' (G12. 5) ' ) 

WRITE( ND_KB_UP(11) , FMT=' (G12. 5) ' ) 

ND KB UP(12) — '************' 

NdIkbIuP(13) 

ND_KB_UP( 14) 

ND_KB_UP(15) 

ND_KB_UP(20) 

WRITE( ND_KB_UP(16), 

WRITE ( ND_KB_UP(17), 

WRITE ( ND_KB_UP(18), 

WRITE( ND_KB_UP(19), 

DO 10 I = 10,20,1 
DY = 18. 6 -. 4*1 

CALL UIS$TEXT(VD_ID,7,ND_KB_UP(I),8. 2,DY) 
CONTINUE 
RETURN 
END 



*?V?V********Vc ' 

FMT='(G12.5)') 
Fbfr='(G12.5)') 
FMT='(G12.5)') 
FMT='(G12.5)') 



ND_FUEL 
ND.TIME 
ND_V0L 
ND PCT 



C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 



C 



SUBROUTINE ND_CALC_2 



SUBROUTINE ND_CALC_2 JUST DOES THE CALCULATION OF THE ENDURANCE* 
REQUIREMENTS FROM THE INPUT DATA. IT IS BASED UPON THE NAVY 
DDS9400-1 ALGORITHM FOR ENDURANCE ESTIMATION. IT IS ALSO USED * 
IN THE REPORT SECTION FOR GENERATION OF THE WRITTEN REPORT OF * 
THE DESIGN. * 

CALLED BY SUBROUTINE 'ND_CALC_1' * 

CALLS SUBROUTINE 'ND FAIL' 



INCLUDE 'TOP. FOR' 
INCLUDE 'TOP_ENDU.FOR' 



C*****CALCULATE AVERAGE ENDURANCE POWER 



ND_C_AVG = ND_NDPWR * 1. 10 
C*****CALCULATE POWER RATIO FOR USE LATER 
ND_RATI0 = ND_C_AVG / ND_PWR 
C*****CALCULATE PROPULSION F/0 CONSUMPTION RATE 
ND_PFRATE = ND_MFR * ND_C_AVG 
C*****CALCULATE AUXILIARY GENERATOR F/0 RATE 
ND_GFRATE = ND_EFR * ND_ELEC 
C*****CALCULATE TOTAL F/0 USEAGE RATE 
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ND_TFR = ND_PFRATE + ND_GFRATE + ND_OTH 
C***-v*CALCULATE ALL PURPOSE F/0 RATE 
ND_APFR = ND_TFR / ND_C_AVG 
C*****SET FUEL RATE CORRECTION FACTOR 

IF ( ND_C_AVG .LE. 0. 3333*ND_PWR ) THEN 
FRCORR = 1. 04 

ELSEIF ( ND_C_AVG . GT. 0. 3333*ND_PWR . AND. 

& ND_C_AVG .LE. 0. 6667*ND PWR ) THEN 

FRCORR = 1. 03 

ELSEIF ( ND_C_AVG . GT. 0. 6667*ND_PWR .AND. 

& ND_C_AVG .LE. ND_PWR ) THEN 

FRCORR = 1. 02 

ELSE 

FRCORR = 1. 03 
ENDIF 

c*****calculate specified f/0 rate 

ND_SFR = FRCORR * ND.APFR 
C*****CALCULATE AVERAGE ENDURANCE F/0 RATE 
ND_AVGFR = ND_SFR * 1. 05 
C*****DETERMINE BURNABLE FUEL LOAD 

ND_WGT = ( ND_RANGE * ND_C_AVG * ND_AVGFR ) 

& / ( ND_SPD * 2240. 0 ) 

c*****adjust for unuseable fuel below suction 

ND.FUEL = ND_WGT / ND_TPA 

c*****calculate final output variables 

ND_VOL = ND_FUEL * 38. 00 

ND_PCT = ND_FUEL * 100. 0 / ND_DISP 

ND_TIME = ND_RANGE / ND_SPD 

C*****whAT if we are over DRASTICALLY OR JUST A LITTLE? 
IF ( ND_FUEL . GT. 0. 150*ND_DISP ) CALL TOO_MUCH 
IF ( ND_FUEL .GT. ND_DISP ) CALL MUCH_TOO_MUCH 
RETURN 
END 
C 
C 



SUBROUTINE ND.RECORD 

C * 

C SUBROUTINE 'ND_REC0RD' SAVES THE KBD INPUT DATA TO A USER * 

C SPECIFIED FILE. FILE TYPE IS '.DAT' AUTOMATICALLY * 

C * 

C CALLED BY SUBROUTINE 'ENDURANCE' * 

C CALLS SUBROUTINES 'SHOW_SAVE' * 

C 



INCLUDE 'TOP. FOR' 

INCLUDE 'TOP_ENDU. FOR' 

STA = UIS$SET_POINTER_POSITION(VD_ID,WD_NDUR,0. ,10. ) 
CALL SHOV_SAVE 

OPEN(23,FILE=ND_INFILE,STATUS='UNKNOWN' ,blank=' null ' ) 

WRITE (23,200) ND_INFILE 

WRITE (23,210) ND_DISP 

WRITE (23,210) ND_PWR 

WRITE (23,210) ND_NDPWR 

WRITE (23,210) ND_ELEC 

WRITE (23,210) ND_TPA 
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WRITE (23,210) ND_MFR 
WRITE (23,210) ND_EFR 
WRITE (23,210) ND_0TH 
WRITE (23,210) ND_RANGE 
WRITE (23,210) ND_SPD 
WRITE (23,210) ND_FUEL 
WRITE (23,210) ND_TIME 
WRITE (23,210) ND_V0L 
WRITE (23,210) ND_PCT 
ENDFILE (23) 

REWIND 23 
CL0SE(23) 

C 

C*****anD REPORT OUT WHEN DONE 
C 

ND_FNAME = ND_INFILE( : INDEX(ND_INFILE, ' ’))//'. DAT' 
CALL UIS$ERASE(VD_ID,-. 4,10. 1,3. 4,10. 9) 

CALL UIS$TEXT(VD_ID,7, 'FILE SAVED AS',0. ,11.) 

CALL UIS$TEXT(VD_ID,7,ND_FNAME,0, ,10.6) 

200 FORMAT (BN,2X,A) 

210 FORMAT (G12.5) 

RETURN 

END 



SUBROUTINE ND READ FILE 






READ FILE IS A SUBROUTINE WHICH READS THE DATA INPUT FILE * 
SPECIFIED BY KEYBOARD ENTRY 

Vf 

CALLED BY SUBROUTINE 'ENDURANCE' * 

CALLS SUBROUTINES ' ND_SDI_LABEL' , ' ND_CALC_1 ' * 



INCLUDE 'TOP. FOR' 

INCLUDE 'TOP_ENDU. FOR' 

STA=UIS$SET_POINTER_POSITION(VD_ID,WD_NDUR,0. ,10. 2) 
KB_ID=UIS$CREATE_KB( 'SYS$WORKSTATION' ) 

CALL UI S $ ENABLE.KB ( KB_ID , WD_NDUR ) 

CALL UIS$ERASE(VD_ID,-. 4,10. 1,2. 3,10. 9) 

CALL UIS$TEXT(VD_ID,7, 'ENTER THE FILE NAME' ,0. ,11.) 
CALL UIS$TEXT(VD_ID,7, 'AND FILE EXTENSION' , 0. ,10.6) 
CALL KEY_READ(ND_LINE, 'TRUE' ,0. ,12. ,*150) 

OPEN( 23 ,FILE=ND_LINE , STATUS=' UNKNOWN ' ) 

READ (23,200) ND_INFILE 
READ (23,210) ND_DISP 
READ (23,210) ND_PWR 
READ (23,210) ND_NDPWR 
READ (23,210) ND.ELEC 
READ (23,210) ND_TPA 
READ (23,210) ND_MFR 
READ (23,210) ND_EFR 
READ (23,210) ND_OTH 
READ (23,210) ND.RANGE 
READ (23,210) ND_SPD 
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READ (23,210) ND_FUEL 
READ (23,210) ND_TIME 
READ (23,210) ND_V0L 
READ (23,210) ND_PCT 
REWIND 23 

C**** NOW READ IT IN AS CHARACTER DATA TO PRINT TO SCREEN 
READ (23,200) ND_KB_UP(1) 

DO 10 I = 2,12,1 

READ (23,220) ND_KB_UP(I) 

10 CONTINUE 

C 

200 FORMAT (2X,A) 

210 FORMAT (G12.5) 

220 FORMAT (2X,A) 

150 CLOSE(23) 

CALL ND_SDI_LABEL(VD_ID) 

CALL UIS$DISABLE_KB(KB_ID) 

CALL ND_CALC_1 

CALL UIS$ERASE(VD_ID,-. 4,10. 1,3. 4,10. 9) 

CALL UIS$TEXT(VD_ID, 7, 'SELECT AN OPTION' ,0. , 11. ) 
CALL UIS$TEXT(VD_ID,7, ' WITH THE MOUSE ',0. ,10.6) 
RETURN 
END 



SUBROUTINE ND_LIST(ND_NDX) 

f'sV’sWrfrVc”;’ 



* 



THIS SUBROUTINE POPS UP A DIALOG BOX TO INSTRUCT THE USER OF * 
THE AVAILABLE OPTIONS FOR EACH '...NOT QUITE OBVIOUS... ' PARA- * 
METER THAT IS REQUESTED BY THE 'ENDURANCE' SUBROUTINE. THE 
SPECIFIC PARAMETER IS PASSED AS 'ND_CAUSE' AND INDICATES WHICH * 
LINES OF INSTRUCTION ARE TO BE DISPLAYED TO THE USER. * 



CALLED BY SUBROUTINE 'ND_CALC_2' * 

c 

IMPLICIT INTEGER(A-Z) 

INCLUDE ' SYS$LIBRARY; UISENTRY' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

INCLUDE 'TOP_ENDU.FOR' 

CHARACTER*34 TOP_A/ '***** PARAMETER DEFINITION *****'/ 
CHARACTER*34 TOP_B/ ' PLEASE ENTER FOLLOWING INFORMATION'/ 

CHARACTER*34 GO_A / ' PLEASE WAIT ' / 

CHARACTER*34 GO_B /' PROGRAM WILL RESUME'/ 

CHARACTER*34 G0_C /'OR PRESS [HOLD SCREEN] TO SAVE '/ 

REAL Y_POSN 

CHARACTER*36 TELL( 7 ) , LI ST( 7 ) , MLI ST( 7 ) , XLIST( 7 ) 

C*****INITIALIZE THE HEADINGS MESSAGES ARRAY 



TELL( 1) 
TELL(2) 
TELL(3) 
TELL(4) 
TELL(5) 
TELL(6) 



THE FULL P0V7ER RATING AT 
THE CRUISING POWER (TYPICAL) 

THE CRUISING ELECTRICAL LOAD 
STORAGE TANK TAIL PIPE ALLOWANCE 
MAIN PROPULSION MACH FUEL USEAGE 
GENERATOR FUEL CONSUMPTION RATE 



no 



TELL(7) = 'OTHER (MI SC) FUEL CONSUMPTION RATE* 

c*****initialize the instructions messages array 



_ f 



LIST(l) 
LIST(2) 
LIST(3) 
LIST(4) 
LIST(5) 
LIST(6) 
LIST(7) 

C*****INITIALIZE 
MLIST(l) = 
MLIST(2) = 
MLIST(3) = 
MLIST(4) = 
MLIST(5) = 
MLIST(6) = 
MLIST(7) = 

c*****initialize 

XLIST(l) = 
XLIST(2) = 
XLIST(3) = 
XLIST(4) = 
XLIST(5) = 
XLIST(6) = 
XLIST(7) = 

c*****create the 



DESIGN MAXIMUM FULL POWER (SHP) 

'FOR STEADY CRUISING, UNDER NORMAL 
'TO OPERATE NORMAL SHIPS FUNCTIONS 
'FOR %-FUEL ABOVE THE SUCTION POINT' 

'FOR NORMAL CONDITIONS AT CRUISING 
'AT CRUISING ELEC LOAD, WITH NORMAL 
'FOR AUX. EQUIP, BOATS, WINCHES, ETC' 

THE SECONDARY INSTRUCTIONS MESSAGES ARRAY 
'NEEDED FOR PROPULSION REQUIREMENTS' 
'WEATHER CONDITIONS & SEA STATE. 

'AND REQUIRED AUX/HAB/DECK GEAR. 

'OF XFER SYS THAT IS USEABLE. 

'SPEED AND POWER SELECTED ABOVE. 
'AUXILIARY MACHINERY IN USE. 

'NORMALLY REQUIRED BY THE CREW. 

EXAMPLES ARRAY 



'EX: ADMIRALTY POWER + 25% MARGIN 
'EX: 0.75 * FULL POWER (TYPICAL) 

'EX: 500KW (NOMINAL) CRUISING 
'EX: 0.95 FOR BROAD, SHALLOW TANKS 
'EX: 0.40 FOR GAS TURBINE @ 20K SHP 
'EX: 0.02 FOR DIESEL GENERATOR 
'EX: 0.05 FOR OTHER SHIP SERVICES 
DISPLAY FOR THE DIALOG BOX 
VD_SCR=UIS$CREATE_DISPLAY(-5. 0,-2. 6,10. 0,4. 0,26. 0,10. 0) 
C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 



WD_SCR=UIS$CREATE_WINDOW(VD_SCR, 'SYS$WORKSTATION' , 'HELP WINDOW' ) 
C*****C0PY ATTRIBUTE BLOCK 'O' AS BLOCK 'll' AND CHANGE THE FONT SIZE 
CALL UIS$SET_FONT( VD_SCR ,0,27,' DTABER0R03WK00GG0001UZZZZ02A000 ' ) 
C*****C0PY ATTRIBUTE BLOCK 'll' AS BLOCK '28' AND CHANGE TO BOLD FONT 

CALL UIS$SET_FONT( VD_SCR ,11 ,1& ,' DTABER0R03WK00PG0001UZZZZ02A000 ' ) 
C*****siGNAL THE USER TO GET THEIR ATTENTION 



CALL UIS$SOUND_BELL( 'SYS$WORKSTATION' ,4) 

C*****WRITE THE TEXT INTO THE WINDOW AND SPAWN PROCESS TO 'WAIT' 10 SEC 
Y_POSN=3. 0 

CALL UIS$TEXT(VD_SCR,28,TOP_A,-2. 5,Y_POSN) 

Y POSN=Y POSN-.6 



CALL UIS$TEXT( VD_SCR , 28 ,TOP_B ,-2.5 , Y_POSN) 
Y_POSN=Y_POSN-. 6 



CALL UIS$TEXT(VD_SCR,27,TELL(ND_NDX),-2. 5,Y_POSN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SCR,27,LIST(ND_NDX),-2. 5,Y_P0SN) 
Y POSN=Y POSN-.6 



CALL UIS$TEXT(VD_SCR,27 ,MLIST(ND_NDX) , -2. 5 , Y_POSN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SCR,27,XLIST(ND_NDX),-2.5,Y_POSN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SCR,28,GO_A,-2. 5,Y_POSN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SCR,28,GO_B,-2. 5,Y_POSN) 
Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SCR,28,GO_C,-2. 5,Y_P0SN) 

CALL LIB$SPAWN( 'WAIT 00: 00: lO' ) 

CALL UIS$SOUND_BELL( ' SYS$W0RKSTATI0N' ,4) 
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CALL UIS$DELETE_DISPLAY(VD_SCR) 

RETURN 

END 



SUBROUTINE MUCH TOO_MUCH 






* 

THIS SUBROUTINE POPS UP A DIALOG BOX TO INSTRUCT THE USER * 

THAT THE CALCULATED WEIGHT OF FUEL IS LARGER THAN THE VESSELS * 
OWN DISPLACEMENT AND THAT CLEARLY SOME CHANGES HAVE TO BE * 

MADE TO THE INPUT DATA. * 

* 

CALLED BY SUBROUTINE 'ND_CALC_2' * 

Q***^V*^V*^V**Vc'^V*yf^V*****Vf*y«V***^V**^Wf**^V^Wc****yf^V***^V**^V*^VVr**‘A’**yfVf‘A'9V**’:?p**** 

c 

IMPLICIT INTEGER(A-Z) 

INCLUDE ' SYS$LIBRARY: UISENTRY' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

INCLUDE ’TOP_ENDU. FOR' 

CHARACTER*34 MTM_A/'** PARAMETER ERROR HAS OCCURRED **'/ 
CHARACTER*34 MTM_B/' PLEASE NOTE THAT THE SHIPS FUEL '/ 
CHARACTER*34 MTM_C/' REQUIREMENTS EXCEED ITS STORAGE '/ 
CHARACTER*34 MTM_D/ ' ABILITY. ADJUST ITEMS 2 THRU 12 '/ 

CHARACTER*34 MTM_E/ ' PLEASE WAIT '/ 

CHARACTER*34 MTM_F/ ' PROGRAM WILL RESUME'/ 

REAL Y_POSN 

C*****CREATE THE DISPLAY FOR THE DIALOG BOX 

VD_MTM=UIS$CREATE_DISPLAY(-5. 0,-2. 0,10. 0,4. 0,26. 0,10. 0) 
C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 

WD_MTM=UIS$CREATE_WINDOW(VD_MTM, 'SYS$W0RKSTATI0N' , 'HELP WINDOW' ) 
C*-'-***c0PY ATTRIBUTE BLOCK 'O' AS BLOCK '29' AND CHANGE THE FONT SIZE 
CALL UIS$SET_FONT(VD_MTM,0,29, 'DTABER0R03WK00GG0001UZZZZ02A000' ) 
C*****C0PY ATTRIBUTE BLOCK '29' AS BLOCK '30' AND CHANGE TO BOLD FONT 
CALL UIS$SET_FONT(VD_MTM,29,30, 'DTABER0R03WK00PG0001UZZZZ02A000' ) 
C*****SIGNAL THE USER TO GET THEIR ATTENTION 

CALL UIS$SOUND_BELL( 'SYS$WORKSTATION' ,4) 

C*****WRiTE THE TEXT INTO THE WINDOW AND SPAWN PROCESS TO 'WAIT' 10 SEC 
Y_P0SN=2. 8 

CALL UIS$TEXT(VD_MTM,30,MTM_A,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_MTM,29,MTM_B,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_MTM,29,MTM_C,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_MTM,29,MTM_D,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_MTM,30,MTM_E,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TE.XT(VD_MTM,30,MTM_F,-2. 5,Y_P0SN) 

CALL LIB$SPAWN('WAIT 00:00: lO') 

CALL UIS$SOUND_BELL( ' SYS$W0RKSTATI0N' ,4) 

CALL UIS$DELETE_DISPLAY(VD_MTM) 

RETURN 

END 
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c 

c 



SUBROUTINE TOO.MUCH 









c * 

C THIS SUBROUTINE POPS UP A DIALOG BOX TO INSTRUCT THE USER * 

C THAT THE CALCULATED WEIGHT OF FUEL EXCEEDS THE NORMAL ALLOW- 
C ANCE FOR SHIPS OF APPROXIMATELY 15% OF TOTAL DISPLACEMENT. * 

C THUS, CLEARLY SOME CHANGES MUST BE MADE TO THE INPUT DATA. * 

C * 

C CALLED BY SUBROUTINE 'ND _CALC_2' * 






c 

IMPLICIT INTEGER(A-Z) 

INCLUDE ' SYS$LIBRARY: UISENTRY' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

INCLUDE 'TOP_ENDU. FOR' 

CHARACTER*34 TM_A/'** PARAMETER ERROR HAS OCCURRED **'/ 
CHARACTER*34 TM_B/' PLEASE NOTE THAT THE SHIPS FUEL '/ 
CHARACTER*34 TM_C/' REQUIREMENTS EXCEED THE NORMAL '/ 
CHARACTER*34 TM_D/ ' RANGE OF 15% TOTAL DISPLACEMENT. '/ 

CHARACTER*34 TM_E/' PLEASE WAIT '/ 

CHARACTER*34 TM_F/ ' PROGRAM WILL RESUME'/ 

REAL Y_POSN 

c*****create the display for the dialog box 

VD_TM=UIS$CREATE_DISPLAY(-5. 0 , -2. 0 , 10. 0 ,4. 0,26. 0,10. 0) 
C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 

WD_TM=UIS$CREATE_WINDOW(VD_TM, ' SYS$WORKSTATION' , 'HELP WINDOW' ) 
C*****C0PY ATTRIBUTE BLOCK 'O' AS BLOCK '3l' AND CHANGE THE FONT SIZE 
CALL UIS$SET_FONT(VD_TM,0,31, 'DTABER0R03WK00GG0001UZZZZ02A000' ) 
C*****coPY ATTRIBUTE BLOCK '3l' AS BLOCK '32' AND CHANGE TO BOLD FONT 
CALL UIS$SET_F0NT(VD_TM,31,32, 'DTABER0R03WK00PG0001UZZZZ02A000' ) 
C*****siGNAL THE USER TO GET THEIR ATTENTION 
CALL UIS$SOUND_BELL( ’SYS$WORKSTATION' ,4) 

C*****WRITE THE TEXT INTO THE WINDOW AND SPAWN PROCESS TO 'WAIT' 10 SEC 
Y_P0SN=2. 8 

CALL UIS$TEXT(VD_TM,32,TM_A,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_TM,31,TM_B,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_TM,31,TM_C,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_TM-,31,TM_D,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_TM,32,TM_E,-2. 5,Y_P0SN) 

Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_TM,32,TM_F,-2. 5,Y_P0SN) 

CALL LIB$SPAWN( 'WAIT 00: 00: 10' ) 

CALL UIS$SOUND_BELL( ' SYS$WORKSTATION' ,4) 

CALL UIS$DELETE_DISPLAY(VD_TM) 

RETURN 

END 

C 

C 

C 
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SUBROUTINE SET_HELP 

Q*VcVoV*Vir’>V^V'5V***’)V**VrVeVc'>V’)V>V’)V'>V')V'A'^V’>V’>V')V’>V’jV'>VVc*iV***>VVfVr'»V*^V**'>V>V9V'>V'>V*:Wr*Vc'jV*‘jV*‘jVVc^^ 

* 



SET_HELP ALLOWS THE USER TO DEFINE WHETHER OR NOT PARAMETER * 
DEFINITION WINDOWS ARE DISPLAYED DURING KB DATA ENTRY BY * 
SELECTING THIS OPTION ON THE MAIN 'ENDURANCE' MENU, THIS * 
FEATURE SIMPLY TOGGLES THE HELP LEVEL TO ON OR OFF AND CALLS * 
THIS WINDOW DISPLAY ACCORDINGLY. * 

* 



CALLED BY SUBROUTINE ' SET_HELP_LEVEL' * 



c 



IMPLICIT INTEGER(A-Z) 

INCLUDE ' SYS$LIBRARY: UISENTRY' 
INCLUDE 'SYS$LIBRARY:UISUSRDEF' 
INCLUDE 'TOP_ENDU. FOR' 
CHARACTER*34 HL(2) 



CHARACTER*34 SH_A/'**** HELP SYSTEM SETTINGS ****'/ 
CHARACTER*34 SH_B/' CHANGING HELP LEVEL FROM '/ 
CHARACTER*34 SH_C/' TO '/ 

CHARACTER*34 SH_D/' PLEASE WAIT '/ 

CHARACTER*34 SH_E/' PROGRAM WILL RESUME'/ 

REAL Y_POSN 

HL(1) = ' HELP WINDOWS ON ' 

HL(2) = ' HELP WINDOWS OFF ' 

C*****CREATE THE DISPLAY FOR THE DIALOG BOX 

VD_SH=UIS$CREATE_DISPLAY(-5, 0,-2. 0,10. 0,4. 0,26, 0,10. 0) 
C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 

WD_SH=UIS$CREATE_WINDOW(VD_SH, ' SYSSWORKSTATION' , 'HELP WINDOW' ) 
C*****COPY ATTRIBUTE BLOCK 'O' AS BLOCK *33' AND CHANGE THE FONT SIZE 
CALL UIS$SET_FONT( VD_SH ,0,33,' DTABER0R03WK00GG0001UZZZZ02A000 ' ) 
C*****COPY ATTRIBUTE BLOCK '33' AS BLOCK '34' AND CHANGE TO BOLD FONT 
CALL UIS$SET_F0NT(VD_SH,33,34, 'DTABER0R03WK00PG0001UZZZZ02A000' ) 
C*****SIGNAL THE USER TO GET THEIR ATTENTION 
CALL UIS$SOUND_BELL( 'SYS$WORKSTATION' ,4) 

C*****WRITE THE TEXT INTO THE WINDOW AND SPAWN PROCESS TO 'WAIT' 10 SEC 
Y_POSN=2. 5 

CALL UIS§TEXT(VD_SH,34,SH_A,-2. 5,Y_POSN) 

Y POSN=Y POSN-.6 



CALL UIS$TEXT(VD_SH,33,SH_B,-2. 5,Y_POSN) 

Y_POSN=Y_POSN-. 6 

IF ( ND_HELP .EQ. 0 ) THEN 

CALL UIS$TEXT(VD_SH,34,HL(l),-2. 5,Y_POSN) 
ELSEIF ( ND.HELP . EQ. 1 ) THEN 

CALL UIS$TEXT(VD_SH,34,HL(2),-2. 5,Y POSN) 
END IF 



Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SH,33,SH_C,-2. 5,Y_POSN) 

Y_POSN=Y_POSN-. 6 

IF ( ND_HELP .EQ. 0 ) THEN 

CALL UIS$TEXT(VD_SH,34,HL(2) ,-2. 5,Y_POSN) 
ELSEIF ( ND_HELP . EQ. 1 ) THEN 

C ALL UI S $TEXT( VD_SH , 34 , HL( 1 ) , - 2 . 5 , Y_POSN ) 
ENDIF 



Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SH,33,SH_D,-2. 5,Y_POSN) 
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Y_POSN=Y_POSN-. 6 

CALL UIS$TEXT(VD_SH,33,SH_E,-2. 5,Y_P0SN) 

CALL LIB$SPAWN( 'WAIT 00:00:10') 

CALL UIS$S0UND_BELL( 'SYS$WORKSTATION' ,4) 

CALL UIS$DELETE_DISPLAY(VD_SH) 

RETURN 

END 

C 

C 

C 

SUBROUTINE SET_HELP_LEVEL 

Q*********************************Vc^V*>V***************^V‘?V^V***Vc’5V*****^V*V«’* 



C * 
C THIS SUBROUTINE ALLOWS THE USER TO HAVE THE SYSTEM PARAMETER * 
C DEFINITION WINDOWS DISPLAYED OR NOT DISPLAYED, DEPENDING UPON * 
C THEIR KNOWLEDGE AND PREFERENCES AND FAMILIARITY WITH TOOL BOX. * 
C * 
C CALLED BY SUBROUTINE 'ENDURANCE' MAIN MENU * 
C CALL SUBROUTINE 'SET HELP' * 



c 

IMPLICIT INTEGER(A-Z) 

INCLUDE 'SYS$LIBRARY:UISENTRY' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

INCLUDE 'TOP_ENDU. FOR' 

IF ( ND_HELP .EQ. 1 ) THEN 
ND_HELP=0 
CALL SET_HELP 

ELSEIF ( ND_HELP . EQ. 0 ) THEN 
ND_HELP=1 
CALL SET.HELP 
ELSE 

ND_HELP=1 
CALL SET_HELP 
ENDIF 
RETURN 
END 
C 
C 

c 

B. SUBROUTINE ENDURANCE VARIABLES DECLARATIONS FILE 

This file is INCLUDE'd in each of the subroutines of the ENDURANCE program 
to porovide commonality of variable declarations. 

C*****THE DECLARATIONS FOR THE ENDURANCE CALCULTIONS SECTION IN ONE FILE 
C 

CHARACTER ND_LINE*12 ,ND_INFILE*12 
CHARACTER ND_KB_UP( 20)*12 ,ND_FNAME*12 

INTEGER ND_ENTERS ,NLINE_NO,ND_NDX,ND_HELP 

REAL ND.RANGE ,ND_SPD ,ND_DISP ,ND_PWR ,ND_ELEC ,ND_MFR 

REAL ND_EFR , ND_FUEL , ND_TPA , ND_WGT , ND_VOL , ND_PCT 

REAL ND_NDPWR , ND_OTH , ND_C_AVG , ND_RATI 0 , ND.PFRATE , ND_GFRATE 
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REAL ND_TFR , ND_APFR , FRCORR , ND_SFR , ND_AVGFR , ND_TIME 
COMMON /NDUR_VARS 1 / ND_TIMES , ND_KB_UP , ND_HELP 

COMMON /NDUR_VARS2/ ND_RANGE , ND_SPD , ND_DISP , ND_PWR , ND_ELEC , ND_MFR 
COMMON /NDUR_NAMES/ND_EFR , ND_FUEL , ND_TPA , ND_WGT , ND_VOL , ND_PCT 
COMMON /NDUR_RESULTS/ND_NDPWR , ND_OTH , ND_C_AVG , ND_RATIO , ND.PFRATE 
COMMON /NDUR_FILES/ ND_FNAME ,ND_INFILE ,ND_GFRATE ,ND_TIME 
COMMON /NDUR_COEFF/ND_TFR , ND_APFR , FRCORR , ND_SFR , ND_AVGFR 

c*****keep power results handy to pass in incase we need them 

COMMON /PW_RESULTS/ PW_SDNB,PW_SDBB,PW_ADME,PW_ADM 
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APPENDIX F. TOOL BOX ENDURANCE ESTIMATION REPORT 



A. ENDURANCE ESTIMATION REPORT SOURCE CODE 

This program contains the Source Code for the Endurance Estimation report. This 
code has been added to the overall TOOL BOX program REPORT. FOR, and can be 
called from the main Report menu. 



SUBROUTINE ENDU_RPT 

* 

SUBROUTINE ENDU_RPT READS IN A STORED DATA FILE FROM THE * 
ENDURANCE SECTION AND PRINTS A DETAILED REPORT OF THE DESIGN * 
PARAMETERS AND THEIR UNITS. * 









CALLED BY SUBROUTINE 'REPORT' * 

INCLUDE 'TOP. FOR' 

INCLUDE 'TOP_ENDU. FOR' 

CHARACTER ND_INPUT_DATA*12 , ND_PRFILE*12 
MOVE THE MOUSE POINTER AWAY FROM THE SELECTION 
STA = UIS$SET_P0INTER_P0SITI0N(VD_ID,WD_RPT,9. 9,23. 5) 

SET UP THE FIRST SCREEN INSTRUCTIONS 
CALL UIS$ERASE(VD_ID, 4.1,20.1,8.4,24.9) 

CALL UIS$TEXT(VD_ID,7, 'ENTER THE FILE NAME ' .4. 5 , 24. ) 

CALL UIS$TEXT(VD_ID,7, 'AND EXTENSION OF THE^ 4. 5 , 23. 6) 

CALL UIS$TEXT(VD_ID,7, 'DATA FILE TO BE USED' ,4. 5 , 23. 2) 

CALL UIS$TEXT(VD_ID,7, 'OR [RETURN] TO EXIT' ,4. 5 ,22. 8) 

READ THE DATA INPUT FILE NAME 

CALL KEY_READ (ND_INPUT_DATA ,. FALSE. ,5. ,21. 2, *30) 

NOW WE NEED THE NAME OF THE OUTPUT REPORT FILE 
CALL UIS$ERASE(VD_ID, 4.1,20.1,8.4,24.9) 

CALL UIS$TEXT(VD_ID,7, ' ENTER THE FILE NAME' ,4. 5,24. ) 

CALL UIS$TEXT(VD_ID,7, ' AND EXTENSION TO BE ,4.5,23.6) 

CALL UIS$TEXT(VD_ID, 7, 'ASSIGNED TO THE REP0RT^4. 5 ,23. 2) 

CALL UIS$TEXT(VD_ID,7,' OR [RETURN] TO EXIT' , 4. 5 , 22. 8) 

CALL KEY_READ (ND_PRFILE , . FALSE. ,5. ,21. 2, *30) 

READ THE DATA IN 

OPEN (23, FILE = ND_INPUT_DATA, STATUS = 'OLD') 

OPEN (24, FILE = ND.PRFILE, STATUS = 'NEW' ,BLANK=' NULL' ) 

READ (23,100) ND_INFILE 
READ (23,110) ND_DISP 
READ (23,110) ND_PWR 
READ (23,110) ND.NDPWR 
READ (23,110) ND_ELEC 
READ (23,110) ND_TPA 
READ (23,110) NDJIFR 
READ (23,110) ND_EFR 
READ (23,110) ND_OTH 
READ (23,110) ND_RANGE 
READ (23,110) ND_SPD 



1 
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READ (23,110) ND_FUEL 
READ (23,110) ND_TIME 
READ (23,110) ND_VOL 
READ (23,110) ND_PCT 
REWIND(23) 

CLOSE (23) 

C*****RE -CALCULATE FOR ACCURACY 
C*****gENERATE the cover page of the REPORT 



DO 5 I = 1,18,1 

WRITE (24,130) ' ’ 

5 CONTINUE 

WRITE (24,120) ' 

WRITE (24,120) ’ ' 

WRITE (24,120) ' 

WRITE (24,120) ’ 

WRITE (24,120) ' 
fiePROGRAM’ 

WRITE (24,120) ' ' 

WRITE (24,120) ' 

&NAVAL ' 

WRITE (24,120) ' 

WRITE (24,120) ’ 

&THE ' 

WRITE (24,120) ’ ' 

WRITE (24,120) ’ 

WRITE (24,120) ’ ' 

WRITE (24,120) ’ 

WRITE (24,120) ' ' 

WRITE (24,120) ' ' 

WRITE (24,120) ' 

WRITE (24,120) ' ' 

WRITE (24,120) ’ 

WRITE (24,120) ' ' 

Vv’RITE (24,120) ' 

WRITE (24,120) ’ 

WRITE (24,120) ' 

WRITE (24,120) ’ ' 

WRITE (24,120) ' ' 

WRITE (24,120) ' ' 

DO 6 I = 1,13,1 

WRITE (24,120) ' ' 

6 CONTINUE 

WRITE (24,120) ' ’ 

WRITE (24,120) ' PAGE 1 OF 
WRITE (24,120) ' ’ 

*****IN THIS SECTION, THE ECHO i 
DATA IS PRESENTED. 



TOOL BOX ' 

ENDURANCE LIMIT REPORT' 

THIS REPORT WAS GENERATED USING THE 

TOOL BOX WHICH WAS DEVELOPED FOR THE 

ENGINEERING DEPARTMENT OF 

NAVAL POSTGRADUATE SCHOOL' 
MONTEREY, CALIFORNIA' 

PROFESSOR F. PAPOULIAS ' 
AND ' 

LT. GERALD MCGOWAN ' 

LT. JAMES PLOSAY' 

1989/90 ' 



/ ENDURANCE LIMIT' 

SHIP PARAMETERS AND CALCULATED 
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DO 10 I = 1,5,1 
WRITE (24,130) 
CONTINUE 
WRITE(24,120) ' 
WRITE(24,130) ' ' 
WRITE(24,120) ' 
WRITE(24,130) ' ' 
WRITE(24,150) 'THE 



ENDURANCE LIMIT REPORT ' 

THE INPUT SHIP PARAMETERS ARE AS FOLLOWS' 
REPORT IS LOCATED IN FILE ; ' ,ND_PRFILE 



US 



100 

110 

120 

130 

135 

140 

150 



WRITE(24,150) 
&ND_INPUT_DATA 
WRITE(24,160) 
WRITE(24,130) 
WRITE(24,135) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
KRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,140) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE(24,130) 
WRITE (24,120) 
WRITE (24,120) 
CLOSE (24) 
FORMAT (A) 
FORMAT (G12.5) 
FORMAT (10X,A) 
FORMAT (10X,A) 
FORMAT (10X,A,' 
FORMAT (10X,A,' 
FORMAT (10X,2A) 



THE INPUT DATA FILE USED IS : ’ , 

DESIGN PARAMETERS' 

smp NAME IS ' ,ND_INFILE 

DESIGN FULL LOAD DISPLACEMENT; ’ ,ND_DISP 

(LTONS)' 

DESIGN FULL POWER LEVEL: ’ ,ND_PWR 

(SHP)' 

DESIGN ENDURANCE POWER LEVEL: ' ,ND_NDPW’R 

(SHP)' 

DESIGN CRUISING ELECTRIC LOAD; ' ,ND_ELEC 

(KW)’ 

DESIGN FUEL TANKS TAIL PIPE ALLOWANCE:,. ’ ,ND_TPA 

(%) ' 

DESIGN MAIN PROPULSION FUEL ECONOMY; ' ,ND_MFR 

(LBS/SHP-HR)' 

DESIGN ELECTRIC PLANT FUEL USE RATE; ' ,ND_EFR 

(LBS/KW-HR)' 

DESIGN OTHER FUEL USE RATES: ' ,ND_OTH 

(LBS/HR)' 

DESIGN RANGE: ' ,ND_RANGE 

(NHILES)' 

DESIGN CRUISING SPEED: ' ,ND_SPD 

(KNOTS)' 

CALCULATED ENDURANCE FUEL LOAD: ' ,ND_FUEL 

(TONS-FUEL OIL)' 

CALCULATED JOURNEY TIME ALLOWANCE: ' ,ND_TIME 

(HRS @ DESIGN SPD & DISTANCE)' 

CALCULATED FUEL STORAGE REQUIREMENTS: . . . ' ,ND_VOL 
(FT 3 VOLUMETRIC STORAGE)' 

CALCULATED FUEL RATIO OF FULL LOAD: ' ,ND_PCT 

i7o OF DESIGN FULL LOAD)' 

I 

I 

I 

I 

ENDURANCE LIMIT VALUES PRESENTED ARE BASED UPON' 

FOLLOWING: ' 

(1) U. S. NAVY DESIGN DATA SHEET DDS9400-1 FORMAT' 
CALCULATIONS USING ESTIMATION FACTORS FOR ' 
FOULING, MACHINERY INEFFICIENCIES, ETC' 

Vc VcVc Vc VcVc VcVc Vc tfWcVrVc Vc Vc Vc VcWVc Vc Vc Vc ^ 

f t 

' PAGE 2 OF 2 / ENDURANCE LIMIT' 



’,A) 

' ,G14. 6) 
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160 FORMAT (10X,A) 

c*****exit from the report and enable mouse 

30 CALL UIS$ERASE(VD_ID, 4.1,20.1,8.4,24.9) 
CALL UIS$TEXT(VD_ID,7, ' SELECT AN OPTION 
CALL UIS$TEXT(VD_ID,7, ’ THE MOUSE', ^ 

RETURN 
END 



WITH’ ,4. 5,24. ) 
.5,23. 6) 
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APPENDIX G. OTHER SOURCE CODE 



A. TOOL BOX MAIN PROGRAM 

This source code is presented for completeness, although much of it belongs to the 
original author, [Ref 3], it is still included here since some of the UIS function calls had 
to be changed to support the new Modules added in this development. 



C 

C 

C* 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c* 



PROGRAM TOOLBOX 
VERSION B18. 01 
VERSION DATE; 09/03/90 

Vc Vc Vc Vc Vr ViT VcVc Vr Vc Vc Vc Vc 

THE GOAL OF VERSION 18 IS TO ADD IN THE SCREEN DISPLAY PLOT 

TO THE PROPULSIVE POWER SECTION, AND TO CLEAN UP THE REMAINING 
CODE SECTIONS FOR FINALIZATION 



REVISION HISTORY 
VERSION B13. 6 
VERSION B14. 0 



VERSION B15. 0 
VERSION B16. 0 
VERSION B17.0 

VERSION B18. 0 

Vc Vr Vr V* Vr Vr ^ 






AS RECEIVED FROM G. MCGOWAN 
ADDED SEVERAL HELP WINDOWS TO EXISTING CODE 
FINALIZED IN B14. 03 
ADDED POWER PREDICTION MODULE 
FINALIZED IN B15. 88 
ADDED POWER PREDICTION REPORT 
FINALIZED IN B16. 09 
ADDED ENDURANCE MODULE 
ADDED USER SELECTION OF HELP LEVEL 
ADDED ENDURANCE ESTIMATION REPORT 
FINALIZED IN B17.42 
ADDED POWER GRAPH TO SCREEN OPTION 
FINALIZED IN B18. 12 

GENERAL HEADER INFORMATION 
INCLUDE 'GENERAL. FOR/LIST' 

EXTERNAL NOWHERE , REPORT , POWER , ENDURANCE 
EXTERNAL DARKl ,DARK2 , LIGHTl , STATIC , EXIT, MANEUVER, UTIL 
REAL DY2 



C 

C*****SINCE THE SUBROUTINES ARE PASSED AS ARGUMENTS IN THE AST ROUTINES 
C THEY MUST BE CALLED IN EXTERNAL STATEMENTS 
C 

C*****CREATE THE VIRTUAL DISPLAY WITH TITLES 
C 

C FIRST, WE SHRINK THE CALLING WINDOW. 

C 

C CALL LIB$SPAWN( ' SET TERMINAL/WIDTH=5 ' ) 

C CALL LIB$SPAWN( ' SET TERMINAL/PAGE=2 ' ) 

C*****THEN CREATE A NEW VIRTUAL DISPLAY 

VD_ID=UIS$CREATE_DISPLAY(-1. ,-l. ,20. ,25.5,40.0,30.0) 

C*****AND ADD COLORS TO IT 

C*****THESE TWO LINES COMMENTED OUT TO ACHIEVE GOOD B&W SCREEN COPIES 
C FOR OVERHEADS. REMOVE COMMENTS IN FINAL VERSION 
C CALL UIS$SET_COLOR(VD_ID,l,. 75,1. ,1. ) 
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C CALL UIS$SET_C0L0R(VD_ID,0,0. ,0. ,.5) 

C*****SET UP SOME CHARACTER TYPES 

CALL U I S $ SET_WRI TING_MODE ( VD_ I D , 0 , 1 , U I S $ C_MODE_RE PLN ) 
CALL U I S $ SET_WRI TI NG_MODE ( VD_ I D , 1 , 2 , U I S $ C_MODE_REPL ) 
CALL UIS$SET_WRITING_M0DE(VD_ID,0,4,UIS$C_M0DE_REPLN) 
CALL UIS$SET_CHAR_SIZE(VD_ID,4,5, 138,. 39) 

CALL UIS$SET_CHAR_SIZE(VD_ID,0,7, 125,. 3) 

c*****initialize 

XO = 0. 0 
XI = 3.0 
Y_LINE = . 7 
FLl = .FALSE. 

FL2 = .FALSE. 

YPOS = 1 

C*****DRAW THE MAIN MENU 

c*****this routine draws the title boxes for the main menu 

DO 10 Y_COOR = 0. 2,5. 8,. 8 
DY = Y_COOR + .4 

CALL UIS$PLOT(VD_ID,0,0. ,Y_COOR,3. ,Y_COOR,3. ,DY,0. , 
& DY,0. ,Y_COOR) 

10 CONTINUE 

C*****sET TITLES FOR MAIN MENU BLOCKS 



BLOCK(l) = ’ 


RESPONSE AREA ’ 


BLOCK(2) = ' 


INSTRUCTIONS ’ 


BLOCK(3) = ' 


OPTIONS MENU ' 


BL0CK(4) = ' 


! 


BLOCK(5) = ' 


SELECTIONS ' 


BLOCK(6) = ' 


f 


BLOCK(7) = ' 


! 


BLOCKC8) = ' 


f 


BLOCK(9) = ' 


! 


BLOCK(IO) = ' 


t 


BLOCK(ll) = ’ 


! 


LONG TITLE(l) 


= ' DATA DISPLAY AREA 


LONG TITLE(2) 


_ t 1 


LONG_TITLE(3) 


_ ! » 


L0NG_TITLE(4) 


_ t t 


LONG TITLE(5) 


_ ! f 


^SET TITLES FOR 


OPTIONS BLOCKS 


OPTION(l) = ’ 


EXIT THE PROGRAM' 


OPTION(2) = ’ 


GENERATE REPORT* 


OPTION(3) = ' 


FILE UTILITIES' 


0PTI0N(4) = ’ 


RESERVE MODULE' 


OPTION(5) = ' 


ENDURANCE ' 


OPTION(6) = ' 


POWER PREDICTION' 


OPTION(7) = ' 


MANEUVERING' 


OPTION(8) = ’ 


STATIC STABILITY' 


OPTION(9) = ’ 


EXIT THE PROGRAM' 


OPTION(IO) = 


' RETURN TO MAIN* 


OPTION(ll) = 


t t 


OPTION(12) = 


t ! 


OPTION(13) = 


' STORE DISPLAYED DATA* 


OPTION(14) = 


' INPUT DATA FROM FILE* 


OPTIONC15) = 


' ENTER DATA FROM KB* 


OPTION(16) = 


' RUN TURNING CIRCLE ' 


OPTION(17) = 


' RUN DYNAMIC PERFORMANCE 
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0PTI0N(18) = 
0PTI0N(19) = 

c*****add options for 

0PTI0N(20) = 
OPTION(21) = 

c*****add options for 

OPTION(22) = 
OPTION(23) = 



LONGITUDINAL STABILITY ' 
INITIAL TRANS. STABILITY' 
VERSIONS 15 & 18 

PLOT GRAPH TO SCREEN ' 
PLOT GRAPH TO DISK ' 
VERSION 17 

CALC ENDURANCE LIMIT’ 
CHANGE HELP LEVEL ON/OFF’ 



C*****NOW FILL THE BLOCKS OF THE MAIN MENU WITH THE STORED TITLES 



DO 20 I = 1,8,1 

CALL UIS$TEXT(VD_ID,0,0PTI0N(I),. 3,Y_LINE) 

Y_LINE = Y_LINE + . 8 
20 CONTINUE 

c*****turn on the main menu viewport 

WD_MAIN=UIS$CREATE_WINDOW(VD_ID, 'SYS$WORKSTATION' , 'MAIN MENU’ , 
& 5,-. 5,10. ,7. ,40. ,30. ) 

C*****sET UP THE AST’S FOR HIGHLIGHTING 



DO 50 Y_COOR = . 2,5. 8,. 8 

DY = Y_COOR + . 2 , 

CALL UIS$SET_P0INTER_AST(VD_ID,WD_MAIN,DARK1, ,X0,Y_C00R,X1, 

& DY,LIGHT1) 

50 CONTINUE 

C*****SET MAIN MENU AST TRAPS FOR MOUSE BUTTONS 

CALL UIS$SET_BUTTON_AST(VD_ID,WD_MAIN,EXIT, ,,X0,. 2,X1,. 6) 

CALL UIS$SET_BUTTON_AST(VD_ID,WD_MAIN, REPORT, , ,XO,l. 0,X1,1. 4) 
CALL UIS$SET_BUTTON_AST(VD_ID,WD_MAIN,UTIL, , ,X0 , 1. 8 ,X1 ,2. 2) 

CALL UIS$SET_BUTTON_AST( VD_ID,WD_MAIN, NOWHERE, , ,X0 , 2. 6 ,X1 , 3. 0) 
CALL UIS$SET_BUTTON_AST(VD_ID,WD_MAIN,ENDURANCE, , ,X0, 3. 4 ,X1 ,3. 8) 
CALL UIS$SET_BUTTON_AST(VD_ID,WD_MAIN, POWER, , ,X0,4. 2,X1,4. 6) 

CALL U I S $ SET_BUTTON_AST( VD_I D , WD_MA I N , MANEUVER ,,,X0,5.0,X1,5.4) 
CALL UIS$SET_BUTTON_AST(VD_ID,WD_MAIN,STATIC, , ,X0,5. 8,X1,6. 6) 
C*****PUT IN THE TITLE AND AUTHORS BLOCK TEXT 
INCLUDE 'HEADER. FOR/LIST’ 

C*****DRAW THE STATIC/MANEUVERING/POWER/ENDURANCE INTERACTIVE DISPLAY 
CALL UIS$PL0T(VD_ID,0,-. 5,10. ,10. ,10. ,10. ,18. 8,-. 5, 

& 18. 8,-. 5,10. ) 

CALL UIS$LINE(VD_ID,0,3. 5, 10. ,3. 5 , 18. 8 , -. 5 , 11. ,3. 5, 11. , 

& -. 5,11. 4,3. 5,11. 4,-. 5,12. 4,3. 5,12. 4,-. 5,12. 8,3. 5,12. 8, 

& -. 5,18. 4,10. ,18. 4) 

DO 30 Y_COOR = 13. , 18.4, .8 

DY = Y_COOR + . 4 
DY2 = DY +. 1 

Z = 9 + INT((Y_COOR - 12. 9)/. 8) 

CALL UIS$TEXT(VD_ID,0,0PTI0N(Z),. 3,DY2) 

C ALL U I S $ PLOT( VD_ I D , 0 , XO , Y_COOR , X 1 , Y_COOR , X 1 , D Y , XO , 

& DY,XO,Y_COOR) 

30 CONTINUE 

C*****FILL THE TITLE BLOCKS OF WORKING DISPLAY 

CALLUIS$TEXT(VD_ID,5,BL0CK(1),-.5,12.85) 

CALL UIS$TEXT(VD_ID,5,BL0CK(2),-. 5,11. 4) 

CALL UIS$TEXT(VD_ID,5,BL0CK(3),-. 5,18. 8) 

CALL UIS$TEXT(VD_ID,5,L0NG_TITLE(1),4. ,18.8) 

C*****THIS LINE CREATES A HARCOPY METAFILE WHEN RUN 
C CALL HCUIS$WRITE_DISPLAY(VD_ID,’HARD. UIS’) 

C*****NOW DRAW UP THE REPORT SCREEN 
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o o 



CALL RPT_GRAPH 

C*****HAVE SYSTEM WAIT FOR USER ACTION 
CALL SYS$HIBER() 

40 END 

C 

C 

C*****eaCH main subroutine is an 'include' FILE TO EASE EDITING AND 
C MODIFICATIONS 

INCLUDE ' REPORT. FOR/LIST' 

INCLUDE 'MANUEVER. FOR/LIST' 

INCLUDE ' STATIC. FOR/LIST' 

INCLUDE ' POWER. FOR/LIST' 

INCLUDE 'ENDURANCE. FOR/ LI ST' 

END 



C 

C 

SUBROUTINE DARKl 



c * 
C THIS SUBROUTINE TURNS THE MENU ITEM IN THE BOX TO REVERSE VIDEO * 
C WHEN THE MOUSE POINTER ENTERS THE REGION DEFINED INSIDE THE * 
C BOX. IT IS USED TO INDICATE TO THE USER THAT THE SELECTION IS * 
C ABLE TO BE SELECTED SINCE THEY ARE IN THE REGION * 
C * 



Q**Vr**Vf**********Vc')V^V**^V*'sV*Vc*Vf***^V***3V')V**Vc*Vf'>V*Vc')V7V*’3V'5V**'5V*'jV****')VVc*^V***'5V’jV'>V'jV* 



c 

INCLUDE 'GENERAL. FOR' 

c***,v*test to see if we have been here before and done the job. 

POSH = UIS$GET_POINTER_POSITION(VD_ID,WD_MAIN,XPOS,YPOS) 

Z = INTCYPOS/. 8)+l 

YO = .8 * INTCYPOS/. 8) + .4 

Y1 = YO + . 3 

IF (FLl) RETURN 

FLl = .TRUE. 

C*****TURN OFF THE LIGHTS 

CALL UIS$TEXT(VD_ID,l,OPTION(Z) ,. 3,Y1) 

RETURN 

END 



SUBROUTINE LIGHTl 



C * 
C THIS SUBROUTINE TURNS THE MENU ITEM IN THE BOX TO REVERSE VIDEO * 
C WHEN THE MOUSE POINTER ENTERS THE REGION DEFINED INSIDE THE * 
C BOX. IT IS USED TO INDICATE TO THE USER THAT THE SELECTION IS * 
C EBLE TO BE SELECTED SINCE THEY ARE IN THE REGION * 
C * 






c 

INCLUDE 'GENERAL. FOR' 

FLl = .FALSE. 

C*****TURN ON THE LIGHTS 

CALL UIS$TEXT(VD_ID,2,OPTION(Z),. 3,Y1) 
RETURN 
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END 

C 

c 

SUBROUTINE DARK2 

Q:iV**^V*'>V*<V*9V'>V*‘5Vx‘>V'>V'5V:^rVcVr*‘>VVcVf'>VVr**->V*9V*9V:V')V*^V:iVVf*Vc*>V9V'3Sr‘>VVr‘>V'>V^V*'>V**’5V'5V:iV*9V*Vc**?V’>V**VoV*')V 



C * 
C THIS SUBROUTINE TURNS THE MENU ITEM IN THE BOX TO REVERSE VIDEO * 
C WHEN THE MOUSE POINTER ENTERS THE REGION DEFINED INSIDE THE * 
C BOX. IT IS USED TO INDICATE TO THE USER THAT THE SELECTION IS * 
C ABLE TO BE SELECTED SINCE THEY ARE IN THE REGION * 
C * 



Q')V***iV'>V**A'iV**')Wc**"A****')V*'>V'A:**'>V**')V*')V>ViV*')V?V*')V**-V')V**iV')V'>V'>V*iV****iV ******* ****** 

C 

INCLUDE 'GENERAL. FOR' 

C TEST TO SEE IF WE HAVE BEEN HERE BEFORE AND DONE THE JOB. 

POSI2 = UIS$GET_POINTER_POSITION(VD_ID,WNDOW,XPOS,YPOS) 

Z2 = INT(YPOS/. 8)-7 

YO = . 8 * INT(YPOS/. 8) + . 3 

Y1 = YO + .4 

IF (FL2) RETURN 

FL2 = .TRUE. 

c*****turn off the lights 

CALL UIS$TEXT(VD_ID,l,OPTION(Z2),. 3,Y1) 

RETURN 

END 

C 

C 



SUBROUTINE LIGHT2 



QVf'5V^V'>Vyr^VVr***Vr*'>VVr'5VVrVrVrVr*Vr*Vc*'5V*Vr*#VVr*^VVr:iVVc'>VVcVr*yr‘A^*Vr*‘5VVcVr9V‘>V*^V')V*VVcVc'3V*'5VVc'>V*5V'5V^V*Vr')V*'>V^V 



c 

c 

c 

c 

c 

c 

c* 

c 



Vc** 



THIS SUBROUTINE TURNS THE MENU ITEM IN THE BOX TO REVERSE VIDEO 
WHEN THE MOUSE POINTER ENTERS THE REGION DEFINED INSIDE THE 
BOX. IT IS USED TO INDICATE TO THE USER THAT THE SELECTION IS 
ABLE TO BE SELECTED SINCE THEY ARE IN THE REGION 



VcVv VcVc Vv Vc Vr Vv Vc V* VcVv Vr VrVcVcVcVcVcVc Vr Vv Vc VcV^Vv Vv ' 



INCLUDE * GENERAL, FOR* 

c*****tell the world we are leaving 

IF (Z2 .LE. ,01) RETURN 
FL2 = .FALSE. 



* 

* 

* 

iV* 



C*****TURN ON THE LIGHTS 

CALL UIS$TEXT(VD_ID,2,OPTION(Z2),. 3,Y1) 

RETURN 

END 

C 

C 

SUBROUTINE EXIT 

Q:VVc^V'5V'5V**Vc^VVr**9V^V9V^V**Vr***VrVr'>V**Vr*Vr*yr^V:V^V'>V*>V'jVV?yr:V**‘5V*Vr**'>VVc'jV'jV‘jV:Vyf'>V‘5VV?*'>V‘5Vyr’jV'>V**yr'A’*'jV 



C * 
C THIS SUBROUTINE IS THE GLOBAL SYSTEM 'EXIT' ROUTINE SELECTED * 
C FROM THE MAIN OR ANY SUB -MENU TO EXIT THE SYSTEM AND RETURN * 
C TO THE CALLING PROCESS * 
C * 



Q*******->V^V*^V*****VoV9VVc‘>V*'>V*‘>V’>V'5VVnVVc')V^V')V'>V'>VVvVc'>V^V^V'3V'5V*^V****'}V'5V*>V'5V*^V***'>V'>V'jV'5V*'5V***‘sV'>V*9V 



c 
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C*****MAKE THE DISPLAY BIG AGAIN SINCE WE SHRUNK IT TO START OFF 
C 

CALL LIB$SPAWN( 'SET TERMINAL/WIDTH=80’ ) 

CALL LIB$SPAWN( ' SET TERMINAL/PAGE=24' ) 

CALL SYS$WAKE(,) 

RETURN 

END 

C 

C 

SUBROUTINE NOWHERE 

■**************•********“*********■*********************************** 



C * 
C THIS SUBROUTINE IS CALLED BY A MOUSE BUTTON SELECTION OF AN * 
C OPTION THAT IS NOT ENABLED YET. IT POPS UP A DIALOG 'HELP' * 
C BOX THAT INSTRUCTS THE USER AND GIVES THEM THE OPTION OF * 
C CONTINUING BACK TO THE PROGRAM AT THE POINT OF INTERRUPTION * 
C OR ENDING THE PROGRAM ENTIRELY. * 
C * 



Q VoV^iWc ****** Vc****:iV ** *Vc** **** *** **** *** ****************** **************** 

c 

IMPLICIT INTEGER(A-Z) 

INCLUDE 'SYS$LIBRARY:UISENTRY' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

CHARACTER*34 HINTl/'*** THIS FEATURE NOT AVAILABLE ***'/ 
CHARACTER*34 HINT2/' TYPE "CONTINUE" TO RESUME '/ 

CHARACTER*34 HINT3/' OR "EXIT" TO END PROGRAM ’/ 

REAL Y_POSN 

C*****CREATE THE DISPLAY WINDOW FOR THE DIALOG BOX 

VD_TEST=UIS$CREATE_DISPLAY(0. 0 , 0. 0 , 10. 0 ,5. 0,20. 0,5. 0) 

C*****CREATE THE WINDOW TO DISPLAY THE HELP TEXT 

WD_TEST=UIS$CREATE_WINDOW(VD_TEST, 'SYS$WORKSTATION' , 'HELP WINDOW' ) 
C*****sET THE SCREEN FONTS WE DESIRE FROM ATTRIBUTE BLOCK 'O' TO BLOCK 
C '20' AND CALL THE FILENAME OF THE FONT 

CALL UIS$SET_F0NT(VD_TEST,0,20, 'DTABER0R03WK00GG0001UZZZZ02A000' ) 

c*****signal the user on process pause for display of help box 

CALL UIS$SOUND_BELL('SYS$WORKSTATION' ,4) 

C*****pRiNT THE TEXT AS WE WANT IT 
Y_P0SN=4. 0 

CALL UIS$TEXT(VD_TEST,20,HINT1,0. 4,Y_P0SN) 

Y_POSN=Y_POSN-l 

CALL UIS$TEXT(VD_TEST,20,HINT2,0. 4,Y_P0SN) 

Y_POSN=Y_POSN-l 

CALL UIS$TEXT(VD_TEST,20,HINT3,0. 4,Y_P0SN) 

PAUSE 

C*****siGNAL THE USER PROCESS RESTARTED 

CALL UIS$SOUND_BELL('SYS$WORKSTATION' ,4) 

c*****return to calling process at point of interrupt 

CALL UIS$DELETE_DISPLAY(VD_TEST) 

RETURN 

END 

C 

C 

SUBROUTINE TO_MAIN 

Q***************'sV*****Vv*******************************************Vc**** 

c * 

C THIS SUBROUTINE IS CALLED TO RETURN TO THE MAIN MENU FROM A * 
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C SUB -MENU. * 

C * 

(]')V***')V***')V********')V****************')Wf**')V*******')V'>V****Vc****'»V********iWf** 

c 

INCLUDE 'GENERAL. FOR' 

C*****DI SABLE THE AST'S 

STA=UIS$SET_POINTER_POSITION(VD_ID,Wndow,9. 7,14. ) 

CALL LIGHT2 

CALL UIS$SET_BUTTON_AST(VD_ID,WNDOW, ,,,X0,13. 0,X1,13.4) 

CALL UIS$SET_BUTTON_AST(VD_ID,WNDOW, , , ,X0,13. 8, XI, 14. 2) 

CALL UIS$SET_BUTTON_AST(VD_ID,WNDOW, , , ,XO , 14. 6 ,X1 , 15. 0) 

CALL U I S $ SET_BUTTON_AST( VD_ ID , WNDOW , , , , XO , 1 5 . 4 , X 1 , 1 5 . 8 ) 

CALL UIS$SET_BUTTON_AST(VD_ID, WNDOW, , , ,X0 , 16. 2 ,X1 , 16. 6) 

CALL UIS$SET_BUTTON_AST(VD_ID, WNDOW, , , ,X0, 17. 0,X1, 17. 4) 

CALL UIS$SET_BUTTON_AST(VD_ID, WNDOW, , , ,X0,17. 8, XI, 18. 2) 
C*****NOW CALL UP MAIN MENU 

CALL UIS$DELETE_WINDOW(WNDOW) 

RETURN 

END 



C 

C 

SUBROUTINE SHOW_SAVE 

C * 

C THIS SUBROUTINE POPS UP A DIALOG BOX TO INSTRUCT THE USER NOT * 

C TO USE THE SAME DATA FILE NAME AS USED IN A PREVIOUS PROGRAM * 

C SECTION [SINCE IT WILL OVERWRITE THE OLD ONE CAUSING ERRORS * 

C OR LOST DATA] . II IT CALLED BY THE SUBROUTINE THAT SAVES THE * 

C DATA FILES 
C 



Q Vf ir * * * * Vr * Vc Vc * * Vc * • 



• V? -sV 1 Wr * Vf -jW? 'sWf * * Vr * -s Wc * * ■sWc'jWc - j V * * Vf ?Wr 'sWc Vc * VcVc -jVVrjV ■? Wc Vf Vf 



•^Wf^V“Wc’^Wr** 



c 

IMPLICIT INTEGER(A-Z) 

INCLUDE 'SYS$LIBRARY:UISENTRY' 

INCLUDE 'SYS$LIBRARY:UISUSRDEF' 

CHARACTER*34 SAVE/'*** SAVING DEFINED DATA FILE *** '/ 
CHARACTER*34 SAVE2/' DO NOT USE SAME NAME & EXTENSION '/ 
CHARACTER*34 SAVE3/'AS DATA FILE FOR OTHER SUB-SECTION'/ 

CHARACTER*34 SAVE4/ ' PLEASE WAIT '/ 

CHARACTER*34 SAVES/' PROGRAM WILL RESUME'/ 

REAL Y_POSN 

C*****CREATE THE DISPLAY FOR THE DIALOG BOX 

VD_SAVE=UIS$CREATE_DISPLAY(0. 0,0. 0,10. 0,5. 0,20. 0,6. 0) 

C*****CREATE THE WINDOW TO DISPLAY THE TEXT IN 

WD_SAVE=UIS$CREATE_WINDOW(VD_SAVE, ' SYS$WORKSTATION' , 'HELP WINDOW' ) 
C*****COPY ATTRIBUTE BLOCK 'O' AS BLOCK '20' AND CHANGE THE FONT SIZE 

CALL UIS$SET_FONT(VD_SAVE,0,20, 'DTABER0R03WK00GG0001UZZZZ02A000' ) 
C*****COPY ATTRIBUTE BLOCK '20' AS BLOCK '2l' AND CHANGE TO BOLD FONT 

CALL UIS$SET_FONT(VD_SAVE,20,21, 'DTABER0R03WK00PG0001UZZZZ02A000' ) 
C*****SIGNAL THE USER TO GET THEIR ATTENTION 

CALL UIS$SOUND_BELL('SYS$WORKSTATION' ,4) 

C*****WRITE THE TE.XT INTO THE WINDOW AND SPAWN PROCESS TO 'WAIT' 15 SEC 
Y_P0SN=4. 0 

CALL UIS$TEXT(VD_SAVE,21,SAVE,0. 2,Y_POSN) 

Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE,20,SAVE2,0. 2,Y_POSN) 
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Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE,20,SAVE3,0. 2,Y_P0SN) 

Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE,21,SAVE4,0. 2,Y_P0SN) 

Y_POSN=Y_POSN-. 8 

CALL UIS$TEXT(VD_SAVE, 21, SAVES, 0. 2,Y_POSN) 

CALL LIB$SPAWN( 'WAIT 00:00:15') 

CALL UIS$SOUND_BELL( 'SYS$WORKSTATION' ,4) 

CALL UIS$DELETE_DISPLAY(VD_SAVE) 

RETURN 

END 

C 

C 

SUBROUTINE UTIL 

*******Vf********Vf*Vc*Vf***Vr*Vf*'}VVc:IV*Vf*Vc*VfVc’}VVc**VcVoVVc*’}V:iV**'5VVe**Vf’)VVf***’}VVc***Vc 



* 

* UTILITY IS A SUBROUTINE WHICH ALLOWS THE USER * 

* TO SHELL BACK TO THE CALLING WINDOW WITHOUT LEAVING THE * 

* PROGRAM ENVIRONMENT. THIS ALLOWS A USER TO RUN A DCL * 

* COMMAND ( SAY A 'DIR' ) WITHOUT LEAVING THE PROGRAM. * 

* WHEN FINISHED, THE USER ENTERS THE COMMAND 'EXIT' OR 

* A RETURN ON A BLANK LINE AND CONTROL RETURNS TO THE TOOL * 

* BOX MENU. * 

* Vf 



^‘V*v^VVf*VfVcVcV?**^V‘5V^V’}V**Vr*}VVc***Vr^V*'jVV?*Vc:iV***VcVc'jVVcVcVc*Vc'5V**jV^V’}VVf*'}V’}V'}V‘VVc7V'jVVc')V’}VVrVf'5VVc'jV‘;V 

INCLUDE 'GENERAL. FOR' 

CHARACTER*64 CMMD 

STA=UIS$SET_POINTER_POSITION(VD_ID,WD_MAIN,9. 9,5. ) 
C*****RETURN THE SCREEN TO DEFAULT SIZE 

CALL LIB$SPAWN( 'SET TERMINAL/WIDTH=80 ' ) 

CALL LIB$SPAWN( 'SET TERMINAL/PAGE=24 ' ) 

C*****PUSH THAT SCREEN TO THE FOREFRONT 
CALL UIS$PUSH_VIEWPORT(WD_MAIN) 

1 TYPE 'ENTER THE DCL COMMAND ...OR [RETURN] TO EXIT’ 

PRINT *, ' ' 

READ (*,FMT=’(BN,A)') CMMD 
PRINT *, CMMD 

C*****IF COMMAND IS A 'RETURN', XFER BACK TO 'TOOL BOX’ 

IF (CMMD .EQ. ' ') THEN 

CALL UIS$POP_VIEWPORT(WD_MAIN) 

RETURN 

ENDIF 

C*****sPAWN A NEW PROCESS FOR THE ENTERED COMMAND 
CALL LIB$SPAWN(CMMD) 

PRINT ’ ’ 

GOTO 1 
END 
C 
C 
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B. TOOL BOX MAIN PROGRAM VARIABLE DECLARATIONS 

This file GENERAL.FOR is INCLUDE'ed in each subroutine of the main core 
code to provide commonality of variable names and types. 



q*****gENERAL. FOR**’''^***’^**’^**’'^*********’^****’’^******’^**********’^*******’^* 
C general’ HEADER INFORMATION 

IMPLICIT INTEGER(A-Z) 

INCLUDE ’SYS$LIBRARY;UISENTRY' 

INCLUDE 'SYS$LIBRARY; UISUSRDEF' 

COMMON /KB/ KB_ID 

COMMON /DISPLAY_1/ VD_ID,VvT)_MAIN,WD_STATIC ,WD_MANU,WNDOW,VD_TEST 

COMMON /DISPLAY_2/ WD_RPT,WD_POVER,WD_FAIL,VD_FAIL,VD_SCR 

COMMON /DISPLAY_3/ WD_NDUR ,WD_TM,WD_MTM,WD_SCR 

COMMON /DISPLAY_4/ WD_SH,VD_SH 

REAL*4 XO , YO , XI , Y1 , XPOS , YPOS , Y_COOR , Y_LINE , DY 

COMMON /REAL_N/ X0,X1,Y0,Y1,Y_C00R,Y_LINE,DY,YP0S 

COMMON /INTEG_R/ Z,Z2 

LOGICAL*! FL1,FL2,KEYBUF(4) 

COMMON /LOGI/FLl,FL2 

CHARACTER*26 BL0CK*29, OPTION, LONG_TITLE*40 
COMMON /CHAR/ BLOCK( 11) ,OPTION( 23) , LONG_TITLE( 5) 



C. TOOL BOX MAIN SUB-PROGRAxM VARIABLE DECLARATIONS 

This file TOP. FOR is INCLUDE'ed in each of the subroutines throughout the 
program that require commonality of only a few of the main program window and de- 
vice variable declarations. 



FOP. FOR '>VVr***>V‘>WciV:»ViV^V^V7V'A’Vc’yc*^Wr‘jV?V'jWoWr*:WcVc‘jV:V^V;WrVv-5Wf «Wf^V*>V:;Sr 
IMPLICIT INTEGER (A-Z) 

INCLUDE ’ SYS$LIBRARY: UISENTRY* 

INCLUDE ' SYS $ LIBRARY: UISUSRDEF' 

COMMON /KB/ KB_ID 

COMMON /DISPLAY_1/ VD_ID , WD_MAIN ,WD_STATIC ,WD_MANU, WNDOW , VD_TEST 
COMMON /DISPLAY_2/ WD_RPT,WD_POWER,WD_FAIL, VD_FAIL, VD_SCR 
COMMON /DISPLAY_3/ WD_NDUR,WD_TM,WD_MTM,WD_SCR 
COMMON /DISPLAY_4/ WD_SH,VD_SH 
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